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ABSTRACT 


This  work  is  concerned  with  the  development  and  application  of  an 
analytical  model  for  "impulsive"  phenomena.  Specifically,  a  new  model 
is  developed  for  atmospheric  noise,  which  is  radio  noise  that  originates 
in  lightning  discharges.  This  model  is  applied  to  the  detection  of  known 
signals  in  additive  atmospheric  noise. 

The  modeling  approach  used  here  is  based  on  the  observation  that  a 
gaussian  model  is  inappropriate  for  received  atmospheric  noise,  primarily 
because  gaussian  noise  does  not  have  the  large  dynamic  range  exhibited 
by  the  received  noise  power.  The  generalized  "t "  model  resulting  from 
the  approach  used  in  this  work  is  in  good  agreement  with  measured  data, 
and  describes  the  received  atmospheric  noise  y(t)  as 

y (t )  =  a(t)  n(t) 

where  n(t)  is  a  narrowband  gaussian  process  and  a(t)  is  a  slowly 
varying  process,  independent  of  n(t),  that  modulates  n(t). 

The  detection  of  known  signals  in  additive  generalized  "t"  noise — 
i.e.,  in  the  presence  of  additive  noise  described  by  the  generalized  "t" 
model,  is  considered,  and  the  receiver  that  minimizes  the  probability  of 
error  is  found  to  be  a  "logarithmic-correlator"  receiver  that  implements 
the  rule: 

Decide  that  the  known  signal  is  present  iff 


where  m(t)  is  the  known  signal,  x(t)  is  the  received  signal,  | x ( t ) ( 
is  the  envelope  of  x(t),  and  |x(t)  -  p(t)|  is  the  envelope  of 
[x(t)  -  m(t)].  An  "exponentially  correct"  estimate  of  the  probability  of 
error  achieved  by  this  receiver  is  given  by 


P  S3  K  exp 
e 


1.35  E 
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2  , 

where  E  is  the  energy  of  the  known  signal,  Y  =  2y  /B  is  proportional 

™  -1  /2 

to  the  average  noise  energy,  and  K  is  a  polynomial  in  (2TB ) 

The  performance  of  the  log-corrclator  receiver  is  compared  with  that 
achieved  by  other  receiver  forms  in  the  presence  of  additive  generalized 
"t"  noise.  Analytical  results  indicate  that  the  log-correlator  receiver 
performs  significantly  better  than  any  linear  receiver  in  the  presence  of 
additive  atmospheric  noise.  Available  experimental  data  support  this 
conclusion . 
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I .  INTRODUCTION 


A.  PURPOSE 

The  purpose  of  this  work  is  two-fold:  1 )  to  develop  an  analytical 
model  for  "impulsive"  phenomena,  i.e.,  phenomena  having  impulse  charac¬ 
teristics,  and  2)  to  apply  this  model  to  the  study  of  communication 
channels  having  additive  impulse  noise.  In  particular,  this  work  focuses 
on  a  specific  impulsive  phenomenon,  atmospheric  radio  noise.  Atmospheric 
radio  noise  is  selected  as  a  case  of  special  interest  because:  l)  there 
is  a  reasonable  amount  of  experimental  data  available  against  which  pro¬ 
posed  models  can  be  checked,  and  2)  additive-atmospheric-noise  channels 
are  important  both  for  communication  purposes  and  as  a  source  of  infor¬ 
mation  on  naturally  occurring  signals  of  scientific  interest.  (An  example 
of  such  signals  is  the  whistler-mode  signal  discussed  in  Chapter  IV.) 

B .  BACKGROUND 

Atmospheric  noise,  which  is  described  in  more  detail  in  Chapter  II, 
is  radio  noise  that  originates  in  lightning  discharges.  At  a  represen¬ 
tative  receiving  site,  lightning  discharges  are  typically  observed  to 
occur  at  a  rate  of  the  order  of  10  discharges  per  second,  although  it  is 
emphasized  that  this  typical  rate  can  vary  from  about  one  per  second  to 
nearly  100  per  second  dependent  on  several  factors  men  tiered  in  C,.  ipte-r 
II.  Since  the  principal  noise  pulse  produced  by  each  discharge  has  a 
duration  of  about  100  ^.sec,  it  follows  that  the  obs-.rvtd  noise  has  a  dis¬ 
tinctive  impulsive  character.  Various  analytical  models  have  been  pro¬ 
posed  for  received  atmospheric  noise,  and  these  usually  are  derived  from 
one  of  two  points  of  view,  which  will  be  discussed  in  more  detail  in 
Chapter  II.  Briefly,  the  most  interesting  approach  from  a  physical  point 
of  view  takes  the  received  noise  to  be  the  weighted  sum  oi  contributions 
from  individual  lightning  discharges.  Although  this  approach  is  we'  1 
motivated  physically,  it  has  the  disadvantage  that  the  resulting  models 
are  not  analytically  tractable.  The  alternative  approach  is  an  empirical 
method  that  yields  analytical  models  chosen  to  fit  the  measured  statistics 
of  the  noise.  However,  in  addition  to  the  lack  of  direct  physicnl  support 
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for  this  procedure,  the  existing  empirical  models  suffer  from  the  fact 
that  they  consider  only  the  first-order  statistics  of  the  noise,  while 
neglecting  the  higher-order  statistics. 

Although  several  models  for  received  atmospheric  noise  have  been 
proposed  in  the  literature,  very  few  analytical  results  are  available  on 
problems  of  signal  detection  or  estimation  m  the  presence  of  additive 
atmospheric  noise.  The  published  results,  which  are  discussed  in  detail 
in  Chapter  III,  fall  into  two  categories: 

1.  A  signal-to-noise  ratio  (SNR)  criterion  has  been  used  to  compare 
specific  receiving  techniques. 

2.  Several  workers  have  employed  empirical  models  to  calculate  the 
probability  of  error  resulting  from  the  use  of  linear  receivers 
in  the  presence  of  additive  atmospheric  noise. 

C.  CONTRIBUTIONS  OF  THIS  WORK 

Using  an  approach  which  differs  from  both  of  those  mentioned  above, 
a  new  analytical  model  for  received  atmospheric  noise  is  developed  in 
Chapter  II.  This  model,  called  the  generalized  "t"  model,  is  attractive 
from  an  analytical  point  of  view  and  also  provides  a  good  fit  to  the  mea¬ 
sured  statistics  of  received  atmospheric  noise,  including  statistics  of 
higher-order  as  well  as  those  of  first-order.  This  model  is  applied  in 
Chapter  III  to  the  detection  of  known  signals  in  additive  atmospheric 
noise.  A  receiver  called  the  "logarithmic-correlator"  receiver  is  found 
to  minimize  the  probability  of  error  in  the  presence  of  additive  general¬ 
ized  "t"  noise.  The  probability  of  error  achievable  using  this  receiver 
is  calculated  and  provides  the  only  theoretical  estimate  known  to  the  au¬ 
thor  of  the  error  performance  achievable  in  the  presence  of  additive  at¬ 
mospheric  noise.  In  Chapter  IV,  the  performance  of  the  log-correlator 
receiver  is  compared  with  that  of  receivers  commonly  used  in  practice  in 
the  presence  of  additive  atmospheric  noise. 


II.  THE  GENERALIZED  "t11  MODEL 


A.  INTRODUCTION 

The  purpose  of  the  work  described  in  this  chapter  is  to  obtain  a 
class  of  nongaussian  random  processes  that  can  serve  as  a  model  for  cer¬ 
tain  impulsive  phenomena  observed  in  nature.  Specifically,  it  is  desired 
to  develop  a  model  for  certain  impulsive  noises  observed  in  communications 
engineering,  with  particular  emphasis  on  the  additive  impulse  noise,  com¬ 
monly  termed  atmospheric  noise,  observed  or  radio  circuits  operating  at 
frequencies  below  approximately  100  Me.  The  work  in  this  chapter  is  con¬ 
cerned  with  the  formulation  of  a  model  for  impulsive  phenomena,  and  the 
verification  of  the  applicability  of  this  model,  when  appropriately  spe¬ 
cialized,  as  a  model  for  received  atmospheric  noise. 

1 .  Physical  Description  of  Atmospheric  Noise 

Before  proceeding  with  the  modeling  problem,  a  brief  physical 
description  of  atmospheric  noise  is  presented  as  follows:  Atmospheric 
noise  is  taken  hire  to  be  radio  noise  that  has  its  source  in  lightning 
discharges.  A  descriptive  discussion  of  both  this  noise  and  its  lightning 
sources  is  presented  by  Watt  and  Maxwell  [Ref.  l],  who  point  out  that  the 
lightning  discharge  consists  essentially  of  a  slowly  developing  leader 
stroke  (predischarge)  of  about  one  msec  duration  followed  by  a  return 
stroke  (main  discharge)  of  about  100  psec  duration.  Although  the  details 
of  the  discharge  are  complicated,  the  observed  noise  can  be  explained  by 
noting  that  the  leader  stroke  is  actually  made  up  of  a  series  of  discrete 
leaders  occurring  at  a  rate  of  one  every  25  to  100  j-sec.  Each  of  these 
leaders  travels  over  a  successive  30  ft  to  200  ft  portion  of  the  cloud-to- 
ground  path,  producing  a  current  pulse  of  about  one  Msec  duration  with  an 
intensity  of  the  order  of  300  amp.  The  result  of  this  recurrent  process 
is  that  the  energy  radiated  by  the  predischarge  has  a  3-db  bandwidth  of 
the  order  of  40  kc  centered  at  a  frequency  of  approximately  30  kc.  Fol¬ 
lowing  this  predischarge,  the  single  return  stroke  takes  place,  producing 
a  current  pulse  of  about  100  ^sec  duration  having  an  average  peak  intensity 
of  20  kiloamperes.  The  energy  radiated  by  this  return  stroke,  which  ac¬ 
counts  for  about  95  percent  of  the  total  energy  radiated  by  the  lightning 
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discharge,  has  a  3-db  bandwidth  of  the  order  of  10  kc  centered  at  a  fre¬ 
quency  of  approximately  10  kc.  This  behavior  of  the  lightning  discharge 
as  a  function  of  time  and  frequency  is  illustrated  in  Fig.  1  (taken  from 
Watt  and  Maxwell)  which  shows  both  the  waveform  and  the  spectrum  of  the 
radiated  field  (see  also  Refs.  2  and  3).  The  total  energy  radiated  by  a 
single  lightning  discharge  is  of  the  order  of  200,000  joules  (Ref.  4 j ; 
therefore  it  is  not  surprising  that  the  resulting  atmospheric  noise  is 
the  predominant  radio  noise  at  vlf  and  If.  However,  because  the  spectrum 
of  the  lightning  discharge  is  approximately  proportional  to  l/f  at  fre¬ 
quencies  above  the  upper  "3-db-down  point,"  as  indicated  in  Fig.  1,  it 
is  often  true  that  atmospheric  noise  is  the  predominant  radio  noise  even 
at  hf .  It  should  be  noted,  however,  that  hf  atmospheric  noise  differs 
from  vlf  atmospheric  noise  in  that  the  former  is  primarily  the  result  of 
the  predischarge  with  its  discrete  fine  structure,  whereas  the  latter  is 
predominantly  the  result  of  the  return  stroke.  With  the  above  description 
of  the  individual  lightning  discharge  in  mind,  and  noting  that  the  average 
number  of  lightning  discharges  observed  per  unit  time  at  a  receiving  site 
typically  lies  in  the  range  from  one  per  second  to  100  per  second,  it  is 
clear,  in  agreement  with  experience,  that  the  received  atmospheric  noise 
will  be  impul.'.ive  in  nature.  It  is  pointed  out  that  the  large  range  ex¬ 
hibited  by  the  average  rate  of  iserved  lightning  discharges  is  the  result 
of  several  factors  detailed  in  the  next  section.  Of  particular  importance 
among  these  factors  is  the  geographical  position  of  the  receiver.  For 
example,  the  average  rate  can  vary  from  about  one  distinguishable  dis¬ 
charge  per  second  at  temperate  latitudes  to  nearly  100  per  second  in  the 
tropics . 

2 .  Summary  of  Existing  Models  for  Atmospheric  Noise 

Turning  attention  now  to  the  problem  of  developing  a  model  for 
received  atmospheric  noise,  it  is  seen  that  this  noise  is  always  observed 
through  the  passband  of  some  receiver  filter.  Thus  the  noise  observed  at 
any  instant  of  time  is  the  resultant  of  several  lightning  discharges  spaced 
in  time  and  location  of  occurrence,  with  the  observed  number  and  intensity 
of  the  lightning  discharges  dependent  upon  a  number  of  factors.  Notable 
among  these  factors  are  receiver  bandwidth  and  band  center  frequency, 
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a.  Waveform 


b.  Spectrum 

FIG.  1.  WAVEFORM  AND  SPECTRUM  OF  THE  RADIATED  FIELD  OF  A 
LIGHTNING  DISCHARGE.  (Referred  to  a  distance  of  one  mile.) 
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receiver  location,  thunderstorm  activity,  and  propagation  conditions. 
These  factors  must  be  taken  into  consideration  in  the  model. 

To  provide  motivation  for  the  development  of  a  new  model  for  at¬ 
mospheric  noise,  a  brief  summary  of  available  models  is  now  presented. 

It  is  reasonable  to  categorize  the  available  models  for  observed  atmo¬ 
spheric  noise  into  two  groups,  to  be  called  here  the  f iltered-impulse 
models  and  the  empirical  models,  respectively. 

a.  The  Filtered-Impulse  Models 

The  f iltered-impulse  models  [Refs.  5  6,  7,  8]  are  well  mo¬ 

tivated  physically  and  are  similar  to  one  another  in  that  they  all  take 
the  received  noise  n(t)  to  have  essentially  the  form 


N 

n(t)  =  ^  >  (2.l) 

i=l 


where  the  are  independent,  identically  distributed  random  variables 

whose  distribution  is  deduced  primarily  from  lightning  discharge  statis¬ 
tics  as  affected  by  receiver  location  and  propagation  considerations;  p(t 
gives  the  form  of  the  noise  pulse  resulting  from  an  individual  lightning 
stroke  as  shaped  by  an  assumed  filter  in  the  front  end  of  the  receiver; 
and  the  t  are  the  occurrence  times  of  the  N  individual  lightning 
strokes.  With  this  formulation  for  the  received  noise,  and  with  the  as¬ 
sumption  that  the  contributions  to  the  received  noise  by  the  individual 
lightning  strokes  are  statistically  independent,  various  results  [Refs.  6, 
7,  8]  have  been  obtained  on  the  first-order  statistics  of  the  received 
noise  using  the  method  of  characteristic  functions.  In  particular,  the 
workers  referenced  above  have  focused  on  the  problem  of  determining  the 
probability  distribution  of  the  envelope  of  the  received  noise,  since 
this  distribution  has  been  measured  quite  extensively  [Refs.  6,  9,  10,  11, 
12],  Although  the  analytical  results  vary  in  their  agreement  with  the 
measured  data  as  a  function  of  both  tne  choice  of  p(t)  and  the  choice 
of  the  statistics  of  N,  of  the  t^,  and  of  the  a^,  the  following  con¬ 
clusions  on  the  f i ltered-impulse  models  appear  to  be  warranted  in  general: 
The  strong  point  of  the  filtered  impulse  approach  is  of  course  that  the 
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model  is  well  motivated  physically.  Thus  it  is  possible  to  specify  the 
model  and  to  determine  the  degree  to  which  the  received  noise  is  expected 
to  be  nongaussian  in  terms  of  physically  meaningful  parameters.  On  the 
other  hand,  there  are  several  disadvantages  that  arise  as  a  consequence 
of  the  f iltered-impulse  approach.  These  can  be  summarized  as  follows: 

1.  The  resulting  probability  distributions  are  quite  complicated,  and 
one  is  not  able  to  put  them  in  closed  form  except  in  certain  limiting 
situations.  As  a  result  of  the  complexity  of  these  distributions, 
they  are  not  of  much  use  in  the  consideration  o.  the  statistical  de¬ 
tection  problem  at  the  receiver.  The  available  results  [Ref.  7]  are 
achieved  by  resorting  to  a  signal-to-noise  ratio  SNR  argument  as  a 
criterion  for  comparing  specific  receiving  techniques;  it  is  well 
known  that  SNR  is  not  necessarily  a  good  measure  of  the  ability  of 

a  receiver  to  process  digital  information. 

2.  The  assumption  that  the  contributions  of  the  individual  lightning 
strokes  to  the  received  noise  are  statistically  independent  of  one 
another,  which  is  crucial  to  the  solution  for  the  first-order  sta¬ 
tistics  of  n(t),  does  not  appear  to  be  true,  at  least  for  vlf  at¬ 
mospheric  noise.  This  conclusion  follows  from  the  experimental 
results  obtained  by  Watt  and  Maxwell  [Ref.  9j  on  the  distribution 
of  the  time  interval  between  crossings  of  a  fixed  level  b,>  the  en¬ 
velope  of  vlf  atmospheric  noise  in  which  it  has  been  found  that  re¬ 
ceived  noise  pulses  are  usually  statistically  dependent  on  preceding 
pulses . 

The  second  disadvantage  is  particularly  disturbing,  because 
the  strongest  argument  in  support  of  the  f iltered-impulse  models  has  been 
the  closeness  of  their  approximation  to  the  physics  of  the  noise.  In  this 
regard  it  is  noted  that  the  f il tered-impulse  models  in  the  literature  have 
typically  been  based  on  the  assumption,  either  that  the  number  of  noise 
pulses  influencing  the  received  noise  is  known  .Ref.  81,  or  that  the  noise 
pulses  occur  in  a  Poisson  fashion  [Refs,  5,  6,  7j.  The  latter  assumption, 
which  enables  one  to  think  of  the  received  noise  as  «  filtered  Poisson 
process,  seems  quite  reasonable,  but  does  not  admit  the  particular  type  of 
statistical  dependence  between  adjacent  noise  pulses  that  is  usually  ob¬ 
served  in  practice. 

One  approach  that  may  be  useful  in  introducing  this  dependence 
into  the  f iltered-impulse  model  is  that  taken  by  Furutsu  and  Ishida  :kof. 
6],  who  have  considered  the  case  of  Poisson-Poisson  noise.  In  this  case 
the  noise  is  assumed  to  consist  of  clusters  of  noise  pulses,  where  the 
pulses  within  the  cluster  occur  in  a  Poisson  fashion  with  mean  rate  v, 
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while  the  clusters  themselves  occur  in  a  Poisson  fashion  with  mean  rate 
v',  where  v'  <  v.  It  is  to  be  noted  that  Furutsu  and  Ishida  were  pri¬ 
marily  studying  hf  atmospheric  noise,  and  introduced  the  Poisson-Poisson 
noise  model  as  a  means  of  representing  the  fine  structure  of  the  leader 
stroke.  It  has  not  yet  beer,  demonstrated  that  the  Poisson-Poisson  model 
yields  the  distribution  of  inter-level-crossing  intervals  observed  by 
Watt  and  Maxwell  at  vlf;  but  it  seems  reasonable  as  a  model  which  imitates 

the  tendency  of  received  atmospheric  noise  to  consist  of  clusters  of  noise 

« 

pulses . 

Finally,  it  should  be  mentioned  as  noted  by  Galejs  [Ref.  8] 
that  the  f iltered-impulse  models  may  be  most  useful  as  a  method  for  study¬ 
ing  the  problem  (the  inverse  of  that  being  studied  here)  of  using  the  rela¬ 
tively  large  amount  of  data  on  received  atmospheric  noise  to  determine  the 
statistics  of  the  lightning  source  itself. 

b.  The  Empirical  Models 

The  empirical  models  differ  fundamentally  in  their  concept 
from  the  f iltered-impulse  models  in  that  they  result  from  an  attempt  to 
construct  a  mathematical  expression  that  fits  the  observed  data  without 
regard  for  the  physics  of  the  noise  source.  In  every  case,  these  empir¬ 
ical  models  consist  of  mathematical  expressions  constructed  in  an  attempt 
to  fit  the  measured  data  on  the  first-order  statistics  of  the  envelope  of 
the  received  noise.  A  recent  summary  of  the  various  empirical  models  that 
have  been  proposed  i.-  presented  by  Ibukun  [Ref.  12],  Despite  the  diver¬ 
sity  of  these  models,  the  following  conclusions  appear  to  be  warranted  in 
general:  The  principal  advantage  that  accrues  from  use  of  the  empirical 
modeling  procedure  is  that  the  resulting  model  for  the  first-order  statis¬ 
tics  of  the  receiv'd  noi:  e  is  much  simpler  than  that  obtained  from  any  of 
the  f il fcered-impulse  models.  Hopefully,  this  simplicity  will  make  the 
empirical  model  useful  in  the  consideration  of  the  signal-detection  prob¬ 
lem  at  the  receiver.  On  the  other  hand,  disadvantages  that  result  from 
use  of  the  empirical  modeling  procedure  are: 

1.  It  may  not  be  possible  to  justify  the  empirical  model  on  physical 
grounds. 

2.  By  the  nature  of  their  construction,  the  empirical  models  consider 
only  the  first-order  statistics  of  the  noise.  In  order  to  extend 
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those  models  to  the  point  that  one  has  a  model  for  the  random  pro¬ 
cess,  it  is  necessary  to  make  further  assumptions  concerning  the 
higner-order  statistics  of  the  noise.  Because  of  the  distinctly 
nongaussian  nature  of  the  noise,  it  appears  that  the  assumptions 
required  to  produce  an  analytically  tractable  model  must  include 
the  assumption  that  the  contributions  of  individual  lightning 
strokes  to  the  received  noise  are  statistically  independent.  As 
discussed  earlier,  this  assumption  is  not  in  agreement  with  measured 
results . 

Although  (l)  above  is  a  fundamental  disadvantage  of  the  empir¬ 
ical  approach,  it  is  noted  that  several  workers  [Refs.  13,  14,  15]  have  had 
some  success  in  finding  physical  justification  for  empirical  models.  Per¬ 
haps  the  most  interesting  of  these,  because  of  its  simplicity  and  its 
closeness  of  fit  to  measured  data,  is  the  model  that  takes  the  envelope 
of  the  received  atmospheric  noise  to  be  Rayleigh-distributed  at  low  values 
of  the  envelope  and  log-normally  distributed  at  high  values.  Beckmann 
[Ref.  14]  has  given  a  physical  argument  which  supports  this  model,  partic¬ 
ularly  in  the  situation  wl  ere  there  is  little  local  thunderstorm  activity. 
It  is  noted  that  several  workers  [Refs.  12,  14,  15,  16]  have  proposed  mod¬ 
els  similar  to  that  considered  by  Beckmann,  although  they  differ  in  regard 
to  how  the  two  distributions  should  be  combined  to  give  the  best  resultant 
model . 

With  this  brief  summary  of  available  models  in  mind,  the  de¬ 
velopment  of  a  new  model  called  the  generalized  "t"  model  will  proceed 
in  the  next  section  from  a  point  of  view  different  than  those  leading  to 
either  the  f iltered-impulse  models  or  the  empirical  models  discussed  above. 

It  is  noted  at  this  point,  however,  that  the  generalized  "t" 
model  will  be  seen  to  give  some  physical  support  to  perhaps  the  simplest 
of  the  empirical  models.  This  empirical  model  is  the  one  [Refs.  12,  17] 
that  takes  the  probability  P  (V )  that  the  noise  envelope  exceeds  the 
value  V  to  be 


po(v>  - 


(2.2) 


where  V  is  the  average  value  of  the  envelope  and  0:  and  r  are  two 
parameters  to  be  chosen.  This  empirical  model  is  also  essentially  identical 
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to  the  model  studied  in  detail  by  Mertz  [Refs.  18,  19],  In  addition  to 
the  fact  that  the  first-order  statistics  of  the  generalized  "t"  model 
agree  closely  with  this  class  of  empirical  models,  it  will  be  shown  that 
the  generalized  "t"  model  has  the  advantage  over  the  empirical  models 
in  general  that  it  can  be  specified  to  also  give  a  good  fit  to  the  higher- 
order  statistics  of  the  noise. 

Before  proceeding  with  the  development  of  a  new  model  for  impul¬ 
sive  phenomena,  it  is  acknowledged  that  the  conceptual  motivation  for  this 
development  derives  from  the  work  of  Mandelbrot  [Ref.  20 ]  on  the  applica¬ 
bility  of  a  class  of  "self-similar"  random  processes  as  a  model  for  cer¬ 
tain  intermittent  phenomena.  Mandelbrot  introduces  the  concept  of  a  random 
process  that  is  controlled  by  one  "regime"  for  the  duration  of  observation, 
where  this  regime  is  itself  a  random  variable.  Although  the  class  of 
"self-similar"  random  processes  considered  by  Mandelbrot  is  not  of  inter¬ 
est  as  a  model  for  atmospheric  noise,  the  regime  concept  will  be  useful 
in  suggesting  a  model  of  the  form  of  the  generalized  "t"  model. 


B.  DEVELOPMENT  OF  THE  NEW  MODEL 

Before  mathematically  formulating  a  new  model  for  impulsive  phenomena, 
a  physical  discussion  focusing  on  the  particular  case  of  received  atmo¬ 
spheric  noise  is  presented  in  order  to  lend  physical  significance  to  the 
mathematical  formulation.  Although  this  discussion  is  based  on  experience 
with  received  atmospheric  noise,  it  is  seen  that  similar  models  may  be 
applicable  to  other  impulsive  phenomena,  e.g.,  to  the  impulse  noise  ob¬ 
served  on  telephone  lines. 

As  mentioned  in  Section  A2,  the  atmospheric  noise  observed  at  an  in¬ 
stant  in  time  through  the  passband  of  a  linear  receiver  is  dependent  on  a 
number  of  factors;  notably  receiver  bandwidth  and  band  center  frequency, 
receiver  location,  thunderstorm  activity,  and  propagation  conditions.  Now, 
if  the  receiver  is  sufficiently  narrowband,  it  is  reasonable  to  assume  that 
the  noise  at  the  receiver  output  is  modeled  well  as  a  gaussian  process. 

This  follows  because  of  the  fact  that  the  filtered  noise  is  the  sum  of  con¬ 
tributions  from  a  large  number  of  independent  lightning  discharges,  none 
of  which  is  dominant  at  the  filter  output.  However,  experimental  data 
[Ref,  9]  indicate  that  the  bandwidth  required  to  achieve  this  condition  at 
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vlf  is  substantially  less  than  one  cps,  so  that  a  gaussian  assumption  is 
not  physically  viable  at  vlf.  The  goal  of  the  model  development  here  is 
the  formulation  of  an  analytical  model  that  is  both  an  accurate  descrip¬ 
tion  of  the  received  noise  and  suitable  for  application  to  the  statistical 
detection  problem  at  the  receiver.  As  far  as  the  detection  problem  is 
concerned,  it  appears  necessary  to  model  the  atmospheric  noise  prior  to 
the  performance  of  any  receiver  operations,  so  that  a  gaussian  assumption 
is  certainly  not  justified. 

Although  the  gaussian  assumption  cannot  be  justified,  the  modeling 
problem  can  nevertheless  be  simplified  by  noting  that  in  practice  atten¬ 
tion  can  be  restricted  to  receiver  bandwidths  substantial. y  smaller  than 
the  band  center  frequency.  This  result,  which  enables  the  received  at¬ 
mospheric  noise  to  be  thought  of  as  a  narrowband  random  process,  follows 
from  the  fact  of  limited  spectrum  availability  at  the  frequencies  where 
atmospheric  noise  is  predominant.  This  limited  spectrum  availability  means 
that  the  receiver  bandwidth  must  be  sufficiently  restricted  to  validate 
the  assumption  that  interference  from  strong  adjacent  channel  signals  is 
negligible.  Further,  it  is  noted  that  this  "narrowband"  assumption  is  pre¬ 
cisely  what  is  needed  to  enable  verification  of  the  accuracy  of  the  analyt¬ 
ical  model,  since  almost  all  avaixable  experimental  data  [Refs.  6,  9,  10, 
11,  12]  have  been  obtained  in  narrowband  situations. 

1 .  Specification  of  the  Generalized  "t"  Model 

With  the  above  general  features  of  the  atmospheric  noise  in  mind, 
the  formulation  of  tne  generalized  "t"  model  proceeds  as  follows:  Al¬ 
though  the  gaussian  assumption  is  not  justified  for  the  atmospheric  noise 
of  interest,  it  is  conjectured,  following  loosely  the  intuitive  notion  of 
regime  suggested  by  Mandelbrot  [Refs.  20,  2 1 J ,  that  this  may  be  closely 
related  to  the  fact  that  the  regime  of  sources  controlling  the  received 
noise  varies  with  time  over  a  large  dynamic  range.  That  is,  the  gaussian 
model  may  be  inappropriate  principally  because  the  received  noise  at  an 
instant  in  time  is  controlled  by  a  few  lightning  discharges  that  vary  in 
intensity  and  location  as  a  function  of  time.  Continuing  with  this  line 
of  reasoning,  it  is  proposed  that  a  model  for  received  atmospheric  noise 
worth  consideration  is  one  that  takes  the  received  noise  to  be  a  narrow- 
band  gaussian  noise  multiplied  by  a  weighting  factor  that  varies  with  time 
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in  accordance  with  the  regime  of  sources  in  control  of  the  noise.  Thus 
it  is  proposed  that  the  narrowband  received  atmospheric  noise  y(t)  be 
considered  to  have  the  form 


y(t)  =  a(t)  n(t)  ,  (2.3) 

where  n(t)  is  a  zero-mean  narrowband  gaussian  process  with  covariance 
function  R^x);  and  a(t),  the  "regime"  process,  is  a  stationary  ran¬ 
dom  process,  independent  of  n(t),  whose  statistics  are  to  be  chosen  so 
that  the  product  a(t)  n(t)  is  an  accurate  description  of  the  received 
atmospheric  noise.  It  is  noted  that  the  above  discussion  leading  to  this 
model  formulation  suggests  that  a(t)  will  be  a  slowly  varying  random 
process  as  compared  to  the  narrowband  process  n(t).  It  will  be  shown 
below  that  this  is  indeed  the  case  in  the  sense  that  for  a  good  fit  to 
measured  data  the  power  spectrum  of  the  lowpass  function  a(t)  needs  to 
have  negligible  overlap  with  the  power  spectrum  of  the  band-centered  gaus¬ 
sian  process  n(t). 

Before  attempting  to  specify  the  statistics  of  the  regime  process 
a(t),  an  alternate  argument  suggesting  the  existence  of  such  a  regime 
process  is  presented  in  support  of  the  model  proposed  in  Eq.  (2.3).  Set¬ 
ting  aside  the  regime  notion  and  observing  the  narrowband  atmospheric  noise 
that  it  is  desired  to  model,  it  is  interesting  to  note  that  a  gaussian 
model  is  inappropriate  primarily  for  the  following  two  reasons  [Refs.  9, 
17]: 

1.  The  atmospheric  noise  process  does  not  tend  to  deliver  energy  at  a 
constant  rate  as  would  a  gaussian  process. 

2.  The  contribution  of  a  single  source  component  may  not  be  small  com¬ 
pared  to  the  total  received  energy. 

Further,  it  is  known  [Ref.  17]  that  in  almost  every  case  where  the  received 
noise  satisfies  requirement  1,  so  that  excursions  far  above  the  rms  value 
are  extremely  unlikely,  the  process  can  be  well  approximated  as  a  gaussian 
process  despite  the  fact  that  it  may  exhibit  a  preferred  waveform  in  vio¬ 
lation  of  requirement  2.  This  is  particularly  true  in  narrowband  situa¬ 
tions  where  the  signal  shape  is  essentially  a  sinusoid  at  the  band  center 
frequency,  and  in  fact  an  example  of  this  result  is  that  in  most  calcula¬ 
tions  a  sinusoid  is  well  approximated  by  a  narrowband  gaussian  process. 
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Thus  it  is  concluded  that  it  is  primarily  the  large  dynamic  range  of  at¬ 
mospheric  noise  that  makes  a  gaussian  model  inappropriate.  Now,  noting 
that  it  is  desirable  from  an  analytical  point  of  view  to  model  the  re¬ 
ceived  atmospheric  noise  in  terms  of  gaussian  processes,  it  seems  naiu. 
to  ask  how  we  can  operate  on  a  gaussian  model  to  introduce  the  large  dy¬ 
namic  range  observed  in  atmospheric  noise.  One  possibility,  which  in  fact 
appears  to  be  the  simplest  feasible  possibility,  is  to  multiply  (modulate) 
the  narrowband  gaussian  process  by  a  process  having  the  desired  dynamic 
range.  Noting  that  this  modulating  process  is  expected  to  vary  as  the 
envelope  of  the  narrowband  atmospheric  noise,  we  are  again  led  to  a  model 
of  the  form  given  by  Eq.  (2.3),  where  a(t)  is  now  thought  of  as  a  slowly 
varying  modulating  function. 

The  above  discussion  based  on  the  physical  characteristics  of  re¬ 
ceived  atmospheric  noise  has  been  presented  to  give  physical  support  to 
the  model  proposed  in  Eq.  (2.3)  and  to  indicate  the  existence  of  the  mod¬ 
ulating  process  a(t).  It  remains  to  specify  the  statistics  of  a(t), 
and  this  crucial  step  will  now  be  considered  from  an  empirical  point  of 
view. 

Recalling  that  a(t)  is  expected  to  be  a  slowly  varying  random 
process,  whereas  n(t)  is  a  narrowband  gaussian  process,  it  is  assumed 
that  the  envelope  v(t)  of  the  received  atmospheric  noise  y(t)  is  given 


V(t)  =  | a ( t ) |  E(t)  , 


(2.4) 


where  E(t)  is  the  envelope  of  n(t).  It  is  now  observed  that  the  first- 
order  statistics  of  the  product  | a ( t ) |  E(t)  have  the  following  "repro¬ 
ducing"  property:  If  the  probability  density  function  of  | a|  has  she 


Hm  p.  ,  (x)  oc  -i- 

X-cc  M  / 


'2.5) 


where  p  >  1,  then  the  probability  density  function  of  V  =  j a|  E  has 
the  form 


lim  py(x)  oc 

X  -*oo  xP 
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Now,  it  is  seen  thnt  the  hyperbolic  distribution  specific  Eq.  (2.6) 
is  asymptotically  identical  in  form  to  perhaps  the  simples  if  the  empir¬ 
ical  models  [Refs.  12,  17,  18,  19]  proposed  from  observation  of  measured 
data  on  the  envelope  of  received  atmospheric  noise  [see  Eq.  (2.2)].  Thus 
it  is  concluded  that  Eq.  (2.5)  gives  a  reasonable  specification  of  the 
asymptotic  behavior  of  the  first-order  statistics  of  a(t). 

It  is  next  observed  that  if  a  random  variable  b  is  defined  to 
be  distributed  according  to  the  "tw<  sided"  chi  distribution 


Pb(b) 


(if 


/2 


m  r/m\ 
a  r\2l 


,  i  m-1 

b|  exp 


m  .2 

I7b. 


-oo  <  b  <  oo  (2.7) 


designated  as  x  (m,cr),^  then  the  random  variable  a  =  l/b  is  distrib- 
uted  according  to 


Pa(a) 


■  4  »>({) 


m  , 


>/2 


'(S) 


i  m+1 


exp 


r  _hl  i_' 

I  o  2  2 

L  2cr  a  _ 


<  a  <  oo  .  (2.8) 


Thus  it  is  seen  that  p  (a)  given  by  Eq .  (2.8)  has  precisely  the  asymp- 

a 

totic  behavior  specified  by  Eq.  (2.5).  Now  it  is  certainly  true  that 
other  distributions  could  be  formulated  having  this  same  asymptotic  be¬ 
havior,  but  that  given  by  Eq.  (2.8)  is  preferred  in  the  model  development 
here  for  the  following  reasons: 


1.  a(t)  distributed  according  to  Eq.  (?-.8)  provides  a  distribution 
for  a(t)  n(t)  that  agrees  very  well  with  available  measured 
data  on  the  first-order  statistics  of  received  atmospheric  noise. 
This  is  demonstrated  in  detail  in  Section  Cl. 

2.  The  model  development  to  this  point  has  considered  only  the  first- 
order  statistics  of  the  noise.  However,  the  fact  that  x2(l,cr) 

is  identical  to  N(o,o2)  means  that  b(t)  =  l/a(t)  is  a  gaussian 
process  in  this  special  case.  Thus  the  model  becomes  the  quotient 


This  designation  follows  from  the  fact  that  | b|  is  distributed  according 
to  the  well-known  chi  distribution  [Ref.  22], 
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of  two  gaussian  processes,  and  it  is  anticipated  that  this  fact  will 
be  useful  in  fitting  the  model  to  the  observed  higher-order  statis¬ 
tics  of  the  noise. 

Now,  with  the  above  specification  of  the  first-order  statistics  of 
a(t)  [or  equivalently  b ( t ) ] ,  a  direct  solution  can  be  obtained  for  the 
first-order  statistics  of  y ( t )  =  a(t)  n(t)  as  follows:  Since  n(t)  and 
b(t)  are  assumed  to  be  statistically  independent,  it  follows  that  [Ref.  22]: 


where 


py(y)  =  J  dx  lxi  pn,b^yX,X^’ 


=  j  dx  |xj  Pn(yx)  Pb(x)  ,  -00  <  y  <  oo  ,  (2.9) 


P  (n)  =  N[0,R  (0)J  ,  R  (0)  =  a\  , 


so  that 


P_(n)  =  -7==^  exp 


2a, 


-co  <  n  <  00 


(2.10) 


and  pb(b)  is  given  by  Eq.  (2.7).  Thus,  using  Dwight  [Ref.  23],  item 
860.17,  this  gives 


pyM 


r/m \  l/2 


'(f) 


l(m+l)/2 


2  U1 

y  +  m  ~o 


m  >  0  , 


(2.11) 
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which  is  conveniently  written  as 


* 


-00  <  y  <  00  (2.12) 


where  y  =  mA '  cr^/a  and  9  =  m+1  >  1.  Now,  noting  that  in  the  special 
case  cr^  =  c,  p^(y)  given  by  Eq.  (2.1l)  is  Student's  "t"  distribution 
(Ref.  22,  p.  180]  with  parameter  m,  it  will  henceforth  be  said  that  the 
first-order  statistics  of  y(t)  are  given  by  the  generalized  "t"  dis¬ 
tribution  defined  by  Eq.  (2.12)  with  parameters  0  and  y. 

Before  proceeding  to  check  the  valiJ<ty  of  the  generalized  "t" 
model  as  a  model  for  atmospheric  noise  by  comparing  it  with  measured  sta¬ 
tistical  data,  there  are  several  points  of  interest  to  note  concerning  the 
model.  In  conjunction  with  the  modulating  function  b(t)  whose  first- 
order  statistics  are  given  by  the  xg  distribution  specified  in  Eq.  (2.7), 

it  is  interesting  to  note  (Ref.  22j  that  if  x.  ,  x  ,  ...,  x  are  indepen- 

X  "2  ^ 

dent  gaussian  variables,  each  distributed  as  N(o,cx  ),  then 


l  k=l 


is  distributed  as  x(n,a).  Thus  it  is  seen  that  |b(t)|  has  a  relation¬ 
ship  to  a  combination  of  gaussian  variables  that  may  be  useful  in  later 
calculations.  It  will  be  shown  below  that  9  in  the  range  2  <  9  <  4  is 
appropriate  to  fit  measured  data  on  atmospheric  .oise,  with  0=3  being 
appropriate  to  fit  a  large  body  of  data  at  vlf  and  If.  Now 


-00  <  y  <  00  , 


(2.14) 
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so  that  E [y ]  j_  =  0  and  var(y)  i  =  Thus  it  is  seen  that  8=3 
results  in  an  infinite  variance,  so  that  this  is  not  a  model  for  a  phys¬ 
ical  noise,  although  it  will  be  found  to  fit  the  data  very  closely.  This 
problem  of  infinite  variance  will  be  considered  in  detail  in  Section  C, 
and  it  will  be  shown  that  there  are  two  possible  solutions  to  the  diffi¬ 
culty: 

1.  The  measured  data  actually  requires  9  >  3  for  an  optimum  fit  at 
large  values  of  y,  so  that  var(y)  is  bounded;  or 

2.  The  generalized  "t"  model  as  formulated  in  this  development  re¬ 
quires  modification  in  that  the  measured  data  depart  from  the  model 
for  cumulative  envelope  probabilities  greater  than  about  one  minus 
10-5.  xt  will  be  shown  that  the  measured  data  in  this  region  of 
envelope  probabilities,  which  is  somewhat  suspect  because  of  the 
tremendous  dynamic  range  required  for  its  observation,  can  be  fitted 
well  by  truncating  pa(a)  given  by  Eq.  (2.8)  at  some  point.  This 
truncation  will  have  negligible  effect  on  py(y)  as  determined 
above  in  Eq.  (2.12)  for  cumulative  probabilities  smaller  than  about 
1-10”5,  so  that  the  truncation  is  not  of  practical  significance  in 
most  problems  of  interest. 

It  should  be  mentioned  that  the  above  discussion,  which  indicates 
that  the  atmospheric  noise  has  very  large  average  power,  agrees  in  princi¬ 
ple  with  the  results  of  Mandelbrot  [Ref.  21 ]  in  which  he  notes  that  inter¬ 
mittent  phenomena  often  appear  to  have  barely  convergent,  or  even  divergent, 
second  moments  dependent  strongly  on  sample  size.  Mandelbrot  [Ref.  21 ] 
also  gives  an  excellent  discussion  of  how  the  mathematical  result  of  in¬ 
finite  variance  is  to  be  interpreted  in  a  physical  application. 

Finally,  it  is  noted  that  p^(y)  as  given  by  Eq.  (2,12)  is  a 
function  only  of  the  two  parameters  8  and  y .  This  result  is  in  agree¬ 
ment  with  measured  data  [Refs.  9,  10,  11 ]  on  the  first-order  statistics 
of  atmospheric  noise,  in  which  it  is  found  that  these  statistics  are  ade¬ 
quately  described  by  the  rms  value  of  the  noise  together  with  the  ratio 
of  this  rms  value  to  the  average  value  of  the  noise  envelope. 


C.  VERIFICATION  OF  APPLICABILITY  OF  THE  GENERALIZED  "t"  MODEL 


Although  a  heuristic  argument  indicating  the  plausibility  of  choosing 
a  model  of  the  type  of  the  generalized  "t"  model  accompanied  the  develop¬ 
ment  of  this  model,  it  remains  to  be  shown  that  a  specific  form  of  the  model 
is  in  fact  of  any  practical  use  as  a  model  for  received  atmospheric  noise. 
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Since  we  are  attempting  to  develop  a  model  for  a  random  process,  the  mea¬ 
surements  required  to  check  the  validity  of  the  model  fall  into  two  cate¬ 
gories,  as  follows: 

The  first  of  these  categories,  and  the  one  for  which  the  greatest 
amount  of  experimental  data  are  available,  is  concerned  with  the  first- 
order  statistics  of  the  random  process.  The  particular  measurements  in 
this  category  that  have  been  reported  in  the  literature  are  measurements 
of  the  probability  distribution  of  the  envelope  of  the  received  noise 
[Refs.  6,  9,  10,  11,  12 ] ,  and  measurements  of  the  average  number  of  level 
crossings  per  unit  time  of  a  fixed  level  by  the  envelope  of  the  received 
noise  [Refs.  9,  10],  Relatively  extensive  measurements  have  been  made  of 
the  probability  distribution  function  of  the  envelope,  and  these  measure¬ 
ments  have  been  used  exclusively  in  the  verification  of  the  empirical  mod¬ 
els  discussed  in  Section  A2.  Thus,  as  the  first  step  in  the  verification 
of  the  applicability  of  the  generalized  "t"  model,  the  probability  dis¬ 
tribution  function  of  the  envelope  will  be  calculated.  Particular  interest 
will  center  on  the  fit  of  the  model  as  a  function  of  band  center  frequency 
and  receiver  bandwidth.  Because  of  the  fact  that  we  are  interested  in 
modeling  the  atmospheric  noise  in  as  wide  a  bandwidth  as  possible  subject 
to  the  interfering  signal  constraint  discussed  in  the  model  development, 
the  model  will  be  checked  using  the  widest  bandwidth  data  available.  It 
is  noted,  however,  that  all  reported  data  have  been  obtained  in  bandwidths 
justifying  a  narrowband  assumption  for  the  received  noise.  In  this  narrow- 
band  case,  the  envelope  V(t)  and  phase  cp(t)  of  the  received  noise 
waveform  have  an  unambiguous  and  physical  (operationally  meaningful)  signif¬ 
icance.  Further,  the  available  data  on  the  first-order  statistics  of  v(t), 
plus  the  intuitive  notion  that  qp(t)  must  be  uniformly  distributed  in  the 
interval  [o,2jt],  provide  the  opportunity  to  check  conclusively  the  appli¬ 
cability  of  the  generalized  "t"  model  as  far  as  the  first-order  statis¬ 
tics  are  concerned. 

As  a  further  check  on  the  first-order  statistics  of  the  model,  the 
average  number  of  level  crossings  by  the  envelope  per  unit  time  will  be 
calculated.  This  calculation  is  of  interest,  because  it  bears  out  the 
experimental  fact  that  the  average  number  of  level  crossings  per  unit  time 
for  the  case  of  atmospheric  noise  is  not  given  by  the  product  of  the 
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probability  density  Junction  ol  the  envelope  with  .a  suitable  bandwidth 
factor,  as  it  would  be  if  the  noise  were  a  gaussian  process  [kef .  24], 

It  is  noted  that  this  calculation  has  not  been  carried  out  for  any  of  the 
empirical  models  oiscussed  in  Section  A2,  although  Nakni  [Ref.  2,1  ]  has 
obtained  numerical  results  in  agreement  with  the  experimental  data  of  Watt 
and  Maxwell  for  a  f il tereo-impulse  model  in  which  the  noise  pulses  occur 
in  a  Poisson  fashion. 

The  second  category  of  measured  data,  and  the  one  where  much  less  data 
are  available,  is  concerned  with  the  second-  and  higher-order  statistics 
of  the  random  process.  The  particular  measurements  that  have  been  reported 
in  this  category  are  measurements  of  the  probability  distribution  of  the 
time  interval  between  crossings  of  a  specified  level  by  the  envelope  of 
the  noise  [Ref.  9].  While  availaole  experimental  data  are  spa’  >e,  being 
restricted  to  a  few  measurements  at  vlf,  it  is  also  true  that  the  analyt¬ 
ical  derivation  of  these  statistics  is  complicated,  requiring  machine  com¬ 
putation  in  the  general  case.  Nevertheless,  this  derivation  will  be 
considered  in  detail  m  Section  C3  and  limiting  cases  will  be  presented  to 
demonstrate  that  the  higher-order  statistics  of  the  generalized  "t"  model 
can  be  selected  to  fit  the  measured  data.  It  is  noted  that  the  interpre¬ 
tation  that  the  generalized  "t"  model  is  made  up  of  a  narrowband  gaussian 
noise  modulated  by  a  slowly  varying  random  process  will  be  useful  in  the 
specification  of  these  higher-order  statistics. 

1 .  First-Order  Distributions  of  Envelope  and  Phase 

As  discussed  above,  the  significance  of  the  representation 

y(t)  =  V(t)  cos^t  +  9(t)j  (2.15) 

for  narrowband  received  noise  is  clear.  In  order  to  check  the  applicabil¬ 
ity  of  the  generalized  "t"  model  as  far  as  the  first-order  statistics 
are  concerned,  it  is  sufficient  to  calculate  the  joint  probability  density 
Pv  -  (v>9)  °f  envelope  and  phase  which  can  be  accomplished  as  follows:  It 
is  well  known  that  [Ref.  26] 

Pv  ~(vi9)  =  v  Pv  ~(v  cos  '  sin  9)  ,  (2. 16) 

v ,  y ,  y 
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where  V  =  (y2  +  y2)1/2,  cp  =  tan  l(y/y),  and  y(t)  is  the  quadrature 
component  corresponding  to  y(t)  which  can  be  found  by  taking  the  Hilbert 
transform^  of  y(t).  Perhaps  the  easiest  way  of  finding  y(t)  here  is 
to  pass  y(t)  through  a  Hilbert  transforming  filter  as  shown  in  Fig.  2. 
Now,  recalling  that 

y ( t )  =  a(t)  n (t )  ,  (2.3) 

where  n(t)  is  a  zero-mean  stationary  narrowband  gaussian  pr^.-tss  and 
a(t)  is  assumed  to  be  a  slowly  varying  stationary  random  process,  it  is 
useful  to  make  the  specific  assumption,  depicted  in  Fig.  3,  that  in  fact 


FIG.  2.  HILBERT  TRANSFORMING  FILTER. 


FIG.  3.  TYPICAL  FOURIER  TRANSFORMS,  A(l )  AND  N(f),  OF  a(t)  AND 
n(t),  (Magnitude  only) 

^See  "Conventions,"  p.  xi. 
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the  Fourier  transform"^  A(f)  of  a(t)  has,  with  high  probability,  neg¬ 
ligible  overlap  in  the  frequency  domain  with  the  Fourier  transform  N(f) 
of  n ( t ) . 

Define  Y(f)  as  the  Fourier  transform  of  y(t);  then  with  the 
above  assumption,  it  follows  directly  that  Y(f),  the  Fourier  transform 
of  y(t) ,  is  given  by 

Y(f)  -  (-i  sgn  f)  Y(f) 


i  J  N(a)  A(f-a)  da  -  i  J  s(a )  A(f-a) 


«  A(f)  *  N(f)  ,  i  =  yfl  , 


(2.17) 


where  *  signifies  the  convolution  of  A(f)  with  N(f).^  Thus  it  is 
found  that,  with  the  assumption  of  negligible  overlap  of  A(f)  and  N(f) 
(which  is  intuitively  reasonable  for  the  noise  of  interest), 


y(t)  »  a (t )  n(t)  . 


(2.18) 


Now,  n(t)  being  a  gaussian  process  implies  that  n(t)  is  a  gaussian 
process  where  it  is  easily  shown  that  n(t)  and  n(t)  at  the  same  time 
instant  t  are  independent  and  identically  distributed  random  variables. 


Thus  it  follows  that 


py,j(y.y)  -  pnn,a;(y’y) 


00 

=  [dxip  ~(x, 

J  2  a»n,n\  x  x/ 


00 

=  J  to  i  pa(x)  P„(j)p^§)  ,  (2.19) 


See  "Conventions,"  p.  xi . 
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v/hero 


p(x)  =  P~(x)  = 
n  n 


i 


and  where  pfl(x)  is  given  by  Eq.  (2.8).  Making  a  change  of  variable  and 
using  Dwight,  item  860.17,  this  yields 


p  ~(v  y)  -  -fa)  7 

pyty  y,y  2jt 


e-i 


~ 2  -v2  2]^®+1^/2 

y  +  y  +  y  J 


(2.20) 


where  y  and  0  are  defined  in  Eq.  (2.12),  Now,  substituting  £q.  (2.20) 
into  Eq.  (2.16): 


=  (e-i) 


9-1 


2jt 


P  *  /] 


(e+i)/2 


t 


=  PV(V)  p  (<p)  0  <  V  <  °°  ,  0<cp<2rt.  (2.2l) 


Thus,  finally: 


p7(v)  =  (e-i)  79”1  — - b;e-i-j/2  .  0<v<»  (2.22) 

Lv  *  y  \ 


and 


P<p(cP)  = 


_1_ 

2jt 


0  (p  <  2jt 


(2.23) 
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a.  Discussion  of  Calculated  Results 


i  m 


The  results  calculated  above  lead  to  the  following  comments 


and  conclusions: 


As  anticipated,  <p( t )  is  found  to  be  uniformly  distributed  in  the 
interval  [0,2n]  in  accord  with  intuition. 


It  is  noted  that  the  results  on  Py(v)  and  Pm(cp)  straightforwardly 
calculated  above  could  have  been  deduced  directly  using  the  formula¬ 
tion  assumed  in  Eq.  (2.4)  along  with  the  reasonable  assumption  that 
the  phase  is  uniformly  distributed  in  [0,2«J. 


3.  Finally,  it  is  important  to  note,  in  support  of  the  generalized  1 
model  as  a  model  for  received  atmospheric  noise,  the  following  as¬ 
ymptotic  forms  of  Py(v): 


lim  p 

V-»  00  V 


(v)  =  (e-i)  79"1  ~ 


(2.24) 


pv(v) 


/  \  V  V 

(0-1 )  —  =  lim  — 2  exp 

y  V+0  o’- 


(2.25) 


where 


a2  0-1 


Thus  it  is  seen  that  for  large  v/y  the  generalized  "t"  model  has 
precisely  the  desired  behavior  specified  by  Eq.  (2.6),  while  for 
small  v/y  the  envelope  of  the  generalized  **t **  model  behaves  as 
if  it  were  Rayleigh-distributed.  This  is  of  course  in  agreement  with 
both  experimental  results  and  intuition,  since  the  noise  at  low  lev¬ 
els  is  expected  to  be  the  result  of  contributions  from  a  large  num¬ 
ber  of  independent  noise  sources. 

Now,  to  proceed  with  the  comparison  of  the  distribution  of  the 
envelope  given  by  the  generalized  "t,!  model  with  available  measured  data, 
it  is  necessary  to  compute  the  complement  of  the  envelope  probability  dis¬ 
tribution  function,  i.e.,  go  find  PQ(v),  where 
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P0(v)  =  1  -  p(v)  =  J  Pv(x)  dx  , 


(2.26) 


in  which  Py(x)  is  given  by  Eq.  (2.22).  Evaluation  of  this  integral 
gives 


po<v>  - 


(v2  *  r2) 


le-D/2  • 


(2.27) 


P  (V)  calculated  as  above  is  plotted  in  Figs  4  through  13 
along  with  measured  data  from  Watt  and  Maxwell  [Ref.  9]  (Figs.  4-9)  and 
Clarke  et  al  [Ref.  10 ]  (Figs.  10-13).  (Note  that  in  the  figures, 
log| log  PQ(v)|  is  plotted  vs  log  V.  This  choice  of  scales  has  the 
interesting  property  that  the  Rayleigh-distributed  envelope  plots  as  a 
straight  line.)  In  addition,  Fig.  14  using  data  from  Furutsu  and  Ishida 
[Ref.  6]  is  presented  to  demonstrate  the  variation  observed  in  the  sta¬ 
tistics  of  the  received  noise  at  hf  where  the  choice  of  an  appropriate 
noise  model  is  strongly  a  function  of  local  thunderstorm  activity. 

In  order  to  produce  these  plots,  it  was  necessary  to  specify 

values  to  be  taken  on  by  the  two  parameters  of  the  model  y  and  9.  The 

plots  in  these  figures  were  produced  by  letting  9  take  on  integer  values 

as  labeled,  while  y  was  chosen  in  the  following  manner:  For  the  case 

9  =  3,  y  was  chosen  to  make  E[v)  =  V  ,  where  E[v]  is  the  expected 

value  of  the  envelope  as  given  by  the  model,  and  V  is  the  time  aver- 

a  vg  ^ 

age  of  the  envelope  of  the  received  noise  obtained  by  direct  measurement. 


v 

E[V]  |  =  f  V  p  (V)  |  dV  , 

9=3  J  9=3 


(2.28) 


^For  Figs,  10  and  11,  Vavg  was  obtained  by  using  the  ratio  ^  g/Va 
presented  by  Clarke,  e_t  a!  in  their  Section  3.1.  ’  “'b 
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Figures  4  through  14 
appear  on  the  following  pages. 
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ENVELOPE  INTENSITY  ,  V  (/x-v/m) 
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Rayleigh-distributed  envelope 
E[vJ  =  5.0  pv/m 
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ENVELOPE  INTENSITY  ,  V  (/tv/m  ) 
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ENVELOPE  INTENSITY 


FIG.  10.  COMPOSITE  PROBABILITY  DISTRIBUTION  FUNCTION  OF  THE  ENVELOPE 
OF  RECEIVED  VLF  ATMOSPHERIC  NOISE:  COMPARISON  OF  MODEL  RESULTS  WITH 
DATA  MEASURED  AT  SLOUGH,  ENGLAND  AND  AT  SINGAPORE. 


P  (v)  =  Probability  that  envelope  intensity  V 
exceeds  abscissa  value. 

[Note:  Vertical  scale  is  log| log  P  (v ) j  . 3 


LEGEND: 


■ ■  Measured  data,  Slough,  England 

Band  center  frequency,  f  =  24  kc 
Receiver  power  bandwidth  =  425  cps 
(Clarke,  ct  al ,  [Ref.  IOj) 

2 

— — —  Measured  data,  Singapore 

(Same  conditions  as  above) 

Generalized  "t"  model 
9=3 

y  =  0.636  V,  V  =  E[v]  «  13  db  below  V 

Modified  generalized  "t"  modol 
9=3 

y  =  0.636  V 
S2o-2  =  2.5  X  106 
5 

♦  ■"  ■—  •  Rayleigh-distributed  envelope 

e[v2J  =  v2 

rms 
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FIG.  11.  COMPOSITE  PROBABILITY  DISTRIBUTION  FUNCTION  OF  THE  ENVELOPE 
OF  RECEIVED  LF  ATMOSPHERIC  NOISE:  COMPARISON  OF  MODEL  RESUl-S  WITH 
DATA  MEASURED  AT  SLOUGH,  ENGLAND  AND  AT  SINGAPORE. 


PQ(V)  =  Probability  that  envelope  intensity  V 
exceeds  abscissa  value. 

[Note:  Vertical  scale  is  log|  log  Pq(v)| .  ] 


LEGEND: 


-  Measured  data,  Slough,  England 

Band  center  frequency,  f  =  135  kc 
Receiver  power  bandwidth  =  425  cps 
(Clarke,  et  £l,  [Ref.  10  j ) 


2 

• . .  Measured  data,  Singapore 

Band  center  frequency,  fQ  =  145  kc 
Receiver  power  bandwidth  =  425  cps 
(Clarke,  et  £l^  [Ref.  10  ] ) 


3 


Generalized  "t"  model 
0  =  3 

y  =  0.636  V,  V  =  E [V ] 


6  db  below  V 

rms 


Generalized  "t"  model 
0  =  2 

y  =  33  db  below  V 
'  rms 


5 


Rayleigh-distributed  envelope 


E[V2]  =  V 


2 

rms 
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FIG.  12.  COMPOSITE  PROBABILITY  DISTRIBUTION  FUNCTION  OF  THE  ENVELOPE 
OF  RECEIVED  HF  ATMOSPHERIC  NOISE:  COMPARISON  OF  MODEL  RESULTS  WITH 
DATA  MEASURED  AT  SLOUGH,  ENGLAND  AND  AT  SINGAPORE. 


P  (v)  =  Probability  that  envelope  intensity  V 
exceeds  abscissa  value. 

[Note:  Vertical  scale  is  log| log  PQ (v ) j  .  ] 


LEGEND: 


Measured  data,  Slough,  England 

Band  center  frequency,  f  =  1]  Me 
Receiver  power  bandwidth  =  425  cps 
(Clarke,  el  jal,  [Ref.  10  ] ) 

Measured  data,  Singapore 

(Same  conditions  as  above) 

Generalized  "t"  model 
e  =  4 

y  =  0.707  V  ,  V*  =  Efv2]  ^  10  db  above  V2 
'  rms  r.  is  rms 

(see  footnote,  page  48) 


Raylei;  a- jistributed  envelope 

E[v2]  =  V2 

rms 
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ENVELOPE  INTENSITY,  V  (db  above  Vrms) 
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FIG.  13.  COMPOSITE  PROBABILITY  DISTRIBUTION  FUNCTION  OF  THE  ENVELOPE 
OF  RECEIVED  HF  ATMOSPHERIC  NOISE:  COMPARISON  OF  MODEL  RESULTS  WITH 
DATA  MEASURED  AT  SLOUGH,  ENGLAND  AND  AT  SINGAPORE. 

» 

PQ(V)  =  Probability  that  envelope  intensity  V 
exceeds  abscissa  value. 

[Note:  Vertical  scale  is  log j log  Pq(v)|.J 


LEGEND: 


Measured  data,  Slough,  England 

Band  center  frequency,  fQ  =  20  Me 
Receiver  power  bandwidth  =  425  cps 
(Clarke,  et  a_l,  [Ref.  10]) 

Measured  data,  Singapore 

(Same  conditions  as  above) 
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Generalized  "t"  model 
0  =  4 

y  =  0.07  V  V2  f  E[v2]  =  10  db  below 

'  rms  rms 

(see  footnote,  page  48) 
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Rayleigh-distributed  envelope 


E[V2]  =  V 


2 

rms 
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rms 
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PROBABILITY  ,  PA  (V)  =  I-  P(V) 


< 


FIG.  14.  PROBABILITY  DISTRIBUTION  FUNCTION  OF  THE  ENVELOPE  OF  RECEIVED 
HF  ATMOSPHERIC  NOISE:  DEMONSTRATION  OF  EFFECT  OF  A  PREDOMINANT  LOCAL 
NOISE  SOURCE. 

P  (V)  =  Probability  that  envelope  intensity  V 
u  exceeds  abscissa  value 

[Note:  Vertical  scale  is  log[  log  P  (v)|.j 

LEGEND: 

Measured  data  (Furutsu  and  Ishida  [Ref.  6]) 

Ohira,  Japan,  0218,  Feb.  19,  1957 
(Predominant  local  noise  source) 

Band  center  frequency,  f^  =  3,5  Me 
Receiver  bandwidth  =  1.2  kc 
a  =  5  db  below  cr  [see  [Ref.  6],  Fig.  8(a)} 

Measured  data  (Furutsu  and  Ishida  [Ref.  6]) 

Ohira,  Japan,  0605,  Feb.  19,  1957 
(No  predominant  local  noise  source) 

Band  center  frequency,  f  =  3.5  Me 
Receiver  bandwidth  =  1.2  kc 
0—10  db  above  aQ  [see  [Ref.  6],  Fig.  9(a)} 

Rayleigh-distributed  envelope 
a  =  V 

rms 


t 
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whore 


p  (V)  i  =  2-/  - i - = 

(v2  ♦  /) 


follows  from  Eq.  (2.22).  Thus  it  is  found  that  E [V ]  I  -  %  y,  so  that, 

0—3  2 

lor  the  case  0  =  3,  y  was  ch<  j  to  satisfy 


7  =  2  V 


n  avg 


(2.29) 


For  the  case  9=4,  because  of  the  fact  that  V  and  not  V  was 

rms  avg 

presented  in  the  measured  data  where  0=4  was  appropriate,  y  was 

2  1  /2  “t 

chosen  to  make  {E[v  ])  '  =  V  .  (it  is  noted  that,  if  available, 

rms 

V  is  preferred  to  V  for  calibration,  since  it  is  less  dependent 
avg  rms  ^ 

on  large  values  of  V,  at.  which  values  the  mensured  data  are  least  reli¬ 
able  because  of  the  tremendous  dynamic  range  and  long  observation  time 
required  to  obtain  valid  results.)  Now 


E[V2] 


1  I  -  /  Vs 

0=4  J 


PV(V)  |  dV 
9=4 


■/  *s7 


3  V~  2 

■-  ,  ,&*  dV  = 

0  (v2  +  /) 


(2.30) 


so  that,  for  the  case  9=4 ,  y  was  chosen  to  satisfy 


,  =  a'1/2  V 


rms 


(2.31) 


» 


« 


1 


For  Figs.  12  and  13,  V’rms  d  Vrms  was  chosen  to  give  the  best  fit,  since 
values  of  V^g  presented  by  Clarke,  et  a_l  in  their  Fig,  6,  for  the  cases 
at  11  Me  and  20  Me,  are  not  compatible  with  each  other. 


SEL-66-052 


-  48  - 


« 


4 


t 


-t 


Finally,  for  the  case  0  =  2,  y  was  chosen  empirically  to  give  the  test 
fit  to  the  measured  data,  since  none  of  the  moments  of  the  envelope  of 
the  generalized  "t"  model  are  finite  in  this  case.  (The  nonexistence 
of  these  moments  is  discussed  further  in  the  next  section.) 

b.  Discussion  of  Plotted  Results 


The  plotted  results  are  reasonably  self-explanatory,  but  sev¬ 
eral  comments  and  conclusions  of  interest  are  warranted.  Because  of  the 
extremely  impulsive  nature  of  the  received  atmospheric  noise  at  vlf ,  it 
is  anticipated  that  the  generalized  "t"  model  for  the  received  noise 
will  be  most  useful  in  this  portion  of  the  spectrum.  Thus  the  plotted 
results  emphasize  the  vlf  region,  the  applicable  figures  being  Figs.  4  to 
10.  Further,  because  of  the  fact  mentioned  in  the  introduction  that  the 
spectral  shape  of  the  received  noise  is  determined  by  the  receiver  in  all 
cases  of  available  measured  data,  and  since  it  is  desired  to  model  the 
noise  as  little  affected  by  the  receiver  as  possible,  the  plotted  data  in 
Figs.  4  to  9  make  use  of  the  widest-band  measured  data  known  to  the  author. 
It  is  noted,  however,  that  a  closely  similar  form  is  obtained  for  the  dis¬ 
tribution  of  the  envelope  for  bandwidths  in  the  range  from  about  100  cps 
up  to  at  least  2  kc,  which  is  the  largest  bandwidth  for  which  any  experi¬ 
mental  results  have  been  reported.  Verification  of  this  statement  is  con¬ 
tained  in  the  experimental  results  of  Watt  and  Maxwell  [Ref.  9]  and  Harwood 
[Ref.  27],  and  it  is  supported  by  Fig.  10,  where  the  data  were  obtained  in 
a  425-cps  bandwidth  rather  than  in  the  1000-to  1300-cps  bandwidths  used  to 
obtain  Figs,  4  to  9.  With  this  explanation  of  the  choice  of  vlf  experi¬ 
mental  data,  the  conclusions  to  note  are: 


At  vlf  the  generalized  "t"  model  is  obviously  a  much  better  fit  to 
the  received  noise  than  is  a  gaussian  model  (Rayleigh  distributed 
envelope).  Quantitatively,  it  is  seen  that  in  most  cases  the  gen¬ 
eralized  "t!'  model  gives  a  good  fit  for  cumulative  envelope  prob¬ 
abilities  smaller  than  about  1-10-5  or  more,  while  the  gaussian  model 
gives  a  good  fit  only  for  cumulative  envelope  probabilities  smaller 
than  about  1-10--. 


2.  It  is  seen  that  a  good  fit  tc  the  measured  vlf  data  is  attained  by 
the  generalized  "t"  model  when  the  parameter  9  is  assigned  an 
appropriate  value  in  the  range  2  <  9  <  4,  and  the  parameter  y 
is  chosen  according  to  an  appropriate  moment  of  the  envelope  given 
9  as  discussed  above.  It  is  interesting  to  note  that  the  choice  of 
9  is  dependent  on  the  location,  number,  and  intensity  of  noise 
sources,  since  in  general,  the  smaller  values  of  9  are  required  to 
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fit  the  Canal  Zone  data  (Figs.  8  and  9)  characterized  by  strong  local 
thunderstorm  activity,  while  the  larger  values  of  £  are  appro¬ 
priate  for  the  Alaskan  data  (Figs,  4  and  5)  characterized  by  much 
less  local  thunderstorm  activity. 

In  addition  to  the  fit  of  the  generalized  "t"  model  at  vlf,  it  is  oi 
interest  to  consider  the  applicability  of  the  model  as  the  band  center 
frequency  f  is  increased.  This  subject  is  considered  in  Figs.  10  to  14 
and  leads  to  the  following  conclusions: 

1.  Although  the  improvement  shown  by  the  generalized  "t"  model  as 

ompared  with  a  gaussian  model  is  most  dramatic  at  vlf  (Fig.  10) 
and  at  If  (Fig.  ll)  it  is  seen  that  the  generalized  "t"  moael 
may  have  application  at  frequencies  up  to  and  including  the  hf  band 
(Figs,  12  and  13).  The  fact  that  the  importance  of  individual  noise 
source  components  decreases  with  increasing  band  center  frequency 
is  evident  in  the  plotted  data  however,  and  Fig.  14  is  presented  to 
emphasize  that  the  choice  between  the  gaussian  model  and  the  gener¬ 
alized  "t"  model  is  not  as  clear  at  hf  as  it  is  at  vlf,  i.e., 

Fig.  14  shows  that  the  statistics  of  hf  atmospheric  noise  are 
strongly  a  function  of  local  thunderstorm  act*  ity. 

2.  It  is  seen,  consistent  with  the  fact  that  the  dynamic  range  of  the 
received  noise  decreases  with  increasing  band  center  frequency,  that 
the  value  of  9  required  to  fit  the  measured  data  increases  with 
band  center  frequency  with  0  >  4  being  appropriate  at  hf  (Figs.  12 
and  13) 

In  summary  of  the  above  discussion,  it  is  concluded  that  the  first-order 
statistics  of  the  generalized  "t"  model  are  in  good  agreement  with  a 
usefully  wide  range  of  measured  data  on  received  atmospheric  noise.  The 
fit  of  the  generalized  "t"  model  is  seen  to  be  especially  good  at  vlf 
and  If  in  those  situations  characterized  by  low-to-modc-rate  local  thunder¬ 
storm  activity. 

With  this  result  in  mind,  it  is  of  interest  to  consider  modifi¬ 
cations  to  the  generalized  "t"  model  that  will  improve  its  fit  when  either 
the  measured  data  is  characterized  by  strong  local  thunderstorm  activity 
or  the  band  center  frequency  of  interest  is  increased  to  hf.  The  plotted 
results  discussed  above  are  useful  in  suggesting  the  appropriate  modifica¬ 
tions  when  the  following  observations  are  made: 

1.  It  is  noted  above  that  the  generalized  "t"  model  produces  a  prob¬ 
ability  distribution  for  the  envelope  that  diverges  from  some  experi¬ 
mental  results  for  cumulative  probabilities  above  about  1-10-5  (see, 
for  example,  Figs.  10  to  13).  It  is  emphasized  that  this  diverge  .  :e 
is  somewhat  questionable  at  vlf  because  of  the  tremendous  dv  ’.a.nic 
range  required  for  the  receiver  to  observe  these  probabilit -  <■  s  without 
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distortion,  so  that  it  may  be  receiver  nonlinearity  that  causes 
this  apparent  divergence.  In  fact,  it  is  seen  that  the  relatively 
wideband  data  in  Figs.  4  to  9  adapted  from  Watt  and  Maxwell  do  not 
exhibit  this  divergence,  whereas  if  is  evident  in  the  relatively 
narrowband  data  m  Fig.  10  adapted  from  Clarke,  e_t  jU  as  well  as  in 
relatively  narrowband  data  presented  by  Watt  ana  Maxwell  [Ref.  9], 
Thus,  despite  the  questionable  existence  of  the  divergence  between 
measured  and  generalized  "t"  model  curves  for  bandwidths  greater 
than  about  1000  cps  at  vlf,  the  fact  that  such  a  divergence  is  seen 
in  narrower  bandwidths  for  cumulative  probabilities  greater  than 
about  1-10”°,  plus  the  fact  that  the  divergence  may  become  more  ev- 
ident  as  the  band  center  frequency  of  interest  is  increased  (see 
Fig.  I'*'  indicate  the  desirability  of  modifying  the  generalized  "t" 
model  tu  accommodate  this  behavior. 

2.  It  is  pointed  out  above  that  the  fit  of  the  generalized  "t"  model 
to  measured  vlf  data  is  optimized  if  0  is  chosen  appropriately  in 
the  range  2  <  0  <  4.  Now,  as  mentioned  in  Section  111,  0  in  the 

range  2  <  0  <  3  means  that  E[y2]  =  a>,  so  that  the  generalized 
"t"  model  with  9  in  this  range  does  not  correspond  to  a  physical 
noise.  This  is  not  disturbing,  since  it  is  expected  [Refs.  18  and 
2l]  that  models  for  "impulsive"  phenomena  will  be  characterized  by 
barely  convergent  or  divergent  second  moments  because  of  the  large 
dynamic  range  which  these  models  must  possess.  Despite  this  argu¬ 
ment,  it  is  desirable  to  consider  modifications  to  the  generalized 
"t"  model  which  satisfy  the  condition  that  every  physical  noise  has 
finite  average  power. 

c.  Development  of  the  Modified  Generalized  "t"  Model 

With  the  above  observations  in  nund,  the  modification  proposed 
here  is  merelj’  that  the  ensemble  of  values  that  the  modulating  process 
a(t)  can  assume  be  truncated  as  follows: 


Pa(a) 


-P  <  a  <  3  , 


(2.32) 


where  k  is  chosen  to  satisfy 


The  physical  justification  for  this  modification  is  simply  that  since  ca.h 
lightning  stroke  must  emit  finite  energy  the  measured  data  must  diverge  at 


\ 


i 


some  level  from  the  generalized  "t"  model  with  0  <  3.  The  above  trun¬ 
cation  of  the  range  of  a(t)  is  perhaps  the  simplest  way  to  introduce 
into  the  generalized  "t"  model  a  change  in  trend  that  will  follow  the 
change  in  trend  of  the  measured  data  for  cumulative  probabilities  above 
about  1-1CT®  (see  Figs.  10  to  14).  Furthermore,  it  is  clear  in  agreement 
with  the  measured  data  that  this  trend  change  will  occur  at  lower  noise 
levels  if  either  the  receiver  bandwidth  is  reduced  or  the  band  center  fre¬ 
quency  is  increased,  since  either  procedure  results  in  reducing  the  impor¬ 
tance  of  any  single  lightning  stroke  to  the  received  noise.  Now,  applying 
the  same  procedure  used  in  conjunction  with  the  generalized  "t"  model 
in  order  to  calculate  the  distribution  of  the  envelope  for  the  modified 
generalized  "t"  model,  it  is  found  that 


t 


j 


(2.33) 


This  integral  is  easily  evaluated  in  closed  form  for  m  even,  and  we  will 
proceed  by  considering  m  =  2,  which  corresponds  to  the  important  case 
0=3.  (See  Appendix  A  for  evaluation  of  Eq.  (2,33)  with  m  odd.)  Thus, 
taking  m  =  2  and  substituting  the  proper  value  for  the  normalizing  con¬ 
stant  k,  it  is  found  using  Dwight,  item  567.9,  that 


P,,(v)  [  =  2V  exp 

0=3 
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0  <  V  <  00  , 
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Now,  calculating  the  complement  of  the  probability  distribu- 
tion  function  for  comparison  with  the  measured  data, 

P0<V>  !  “  2  2  GXp  (  2  2  2  /  '  (2.35) 

9=3  V  r  y  \  &  <r  7  / 

This  result  is  plotted  in  Fig.  10,  where  it  is  seen  that  the  modification 

has  the  desired  effect  of  improving  the  fit  of  the  model  to  measured  data 

-6 

at  cumulative  probabilities  above  about  1-10  without  significantly 
altering  the  fit  at  lower  cumulative  probabilities.  Thus  it  is  concluded 
that  the  modified  generalized  "t"  model  has  flexibility  not  possessed 
by  the  generalized  "t"  model.  This  flexibility  may  be  of  practical 
significance  when  the  dynamic  range  of  the  received  noise  becomes  suffi¬ 
ciently  limited  because  of  either  increasing  band  center  frequency  or 
decreasing  receiver  bandwidth. 

Finally,  it  is  conjectured  that  the  model  can  be  further 
improved  with  regard  to  its  ability  to  fit  measured  data  characterized 
by  strong  local  thunderstorm  activity  by  taking  the  first-order  statis¬ 
tics  of  a(t)  to  be  given  by 


1 

1 


m^  +  i 


m 


P  (a) 

a 


hi  <  h 

(2.3s) 

h  <  M  <  e2 . 


This  second  modification  will  not  be  further  considered  here,  because  it 
does  not  fundamentally  change  the  form  of  the  model  although  it  substan¬ 
tially  complicates  the  calculations  of  interest. 

2 .  Average  Intensity  of  Envelope  Level  Crossings 

As  discussed  earlier,  it  is  of  interest  to  calculate  the  average 
number  of  crossings  of  a  fixed  level  by  the  envelope  of  the  noise  per  unit 
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time,  both  as  a  further  check  on  the  first-order  statistics  of  the  gen¬ 
eralized  "t"  modi)  and  as  a  demonstration  of  the  analytical  usefulness 
of  this  model.  Letting  N y (VQ )  denote  the  average  number  of  crossings 
of  the  level  V  by  the  envelope  per  second,  it  is  well  known  [Ref.  24 j 
that  for  an  ergodic  random  process  V(t) 


p  -(v,v)  |  dV  , 

'  v=v„ 

0 


(2.37) 


where  v(t)  =  d/dt  v(t).  Now,  it  was  shown  in  Section  Cl  that  with  the 
assumption  of  negligible  spectral  overlap  of  a(t)  and  n(t), 


V(t)  ~  J  a  ( t )  |  E(  t )  ,  (2.4) 

where  E(t)  is  the  envelope  of  the  narrowband  gaussian  process  n(t). 
Thus 

V(t)  *  j  a(t  )j  E(t)  +  Ja(t)|  E(t)  .  (2.38) 

However,  noting  that  the  available  experimental  measurements  of  N^Tv^T, 
which  are  presented  in  Fig.  15,  were  obtained  in  a  bandwidth  of  1100  cps, 
it  is  interesting  to  consider  >.he  assumption  that  the  modulating  function 
a(t)  is  sufficiently  slowly  varying  so  that  with  high  probability 

j  a(t)|E(t)  »  |  a(t)|  E(t)  .  (2.39) 


This  assumption  is  of  interest  since  it  means  that  with  high  probability 

V(t)  «  J  a ( t ) |  E(t)  ,  (2.40) 


so  that  one  can  write 
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(2.41) 
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where 


0  <  x  <  oo 


follows  from  Eq.  (2.8)  and  when  p  "(v/x,  V/x)  is  calculated  as  fol- 

£  i  £ 

lows:  Since  E(t)  is  the  envelope  of  a  stationary  gaussian  process  n(t) 
with  covariance  function  R^(t),  it  can  be  shown  [Ref.  26],  making  the 
reasonable  assumption  that  n(t)  has  finite  second  moment,  that 


•  exp 


=  PE(E)  Pj(E)  , 


(2.42) 


2  ~ 

where  R  (o)  =  <t  in  che  previous  notation,  and  R  (t)  is  the  Hilbert 
n  1  r  n 

transform  of  R  (t).  Now,  it  is  interesting  to  note  that  if 
n 


Rn(T)  =  RC(T)  cos  ^0T 


(2.43) 
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AVERAGE  NUMBER  OF  POSITIVE  CROSSINGS/SECOND 


ENVELOPE  PROBABILITY  DENSITY 


FIG.  15.  AVERAGE  RATE  OF  ENVELOPE  LEVEL  CROSSINGS:  COMPARISON  OF  MODEL 
RESULTS  WITH  MEASURED  VLF  DATA. 

LEGEND: 

Measured  data:  A-  rage  rate  of  envelope  level 
crossings.  (Watt  and  Maxwell  [Ref.  9]) 

Point  Barrow,  Alaska,  1450  A.S.T.,  Sept.  27,  1956 
Band  center  frequency,  V  =  2' 

Receiver  6-db  bandwidth  =  1100  cps 
6'  antenna,  N.E.-S.W.  plane 

Measured  data:  Envelope  probability  density  function. 
(Watt  and  Maxwell  [Ref.  9]) 

(Same  conditions  as  above) 

Generalized  ’’t"  model:  Average  rate  of  envelope 
level  crossings. 

9  =  2 

7  =  1.5  X  10-6 
rms  bandwidth,  Bc  =  500  cps 

Generalized  "t"  model:  Envelope  probability  density 
function. 

(Same  conditions  as  above) 


whore  S  (f ),  *ho  Fourier  transform  of  R  (t,,  is  zero  for  ! f I  >  f 
c  c  '  '  —  0 

(which  is  in  fact  the  case  of  interest  in  this  development),  then 


-s„(0>-rfo  • 


Also,  it  is  easily  shown  that  in  this  situation 


Rc(0)  =  (2n)  Rc(0)  =  (2rt)  Rn(o)  BT  , 


is  defined  to  be  the  rms  bandwidth  of  n(t).  Thus,  substituting  these 
results  into  Eq.  (2.42): 


(2.44) 


(2.45) 
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(2 nfT*  3  CT. 
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exp< 
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(2«)2  b2c 


(2.46) 


Now,  substituting  this  result  into  Eq.  (2.41 ),  making  a  change  of  variable, 
and  using  Dwight,  item  860.17,  it  is  found  that 


_.(e+2)/2 


9  >  1  . 


(2  «r  b 


(2.47) 


Finally,  substituting  Eq.  (2.47)  into  Eq.  (2.37),  there  is  obtained  the 
desired  result 


0  <  vo  <  co  ,  e>i  . 

(2.48) 

Now,  the  measured  data  presented  by  Watt  and  Maxwell  [Ref.  9]  and  repro¬ 
duced  in  Fig.  15  are  for  crossings  of  the  level  VQ  in  the  positive 
direction  only,  and  are  seen  to  correspond  to  the  case  0=2,  for  which 
formula  (2.48)  becomes 
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vo +  y 


(2.49) 


a.  Discussion  of  Plotted  Results 

Although  the  basic  agreement  of  the  results  derived  using  the 
generalized  "t"  model  with  the  measured  level-crossing  results  is  clear 
from  inspection  of  Fig.  15,  the  following  comments  and  conclusions  are  in 
order:  Recalling  the  result  due  to  Rice  [Ref.  24]  that  if  the  noise  were 
gaussi an  the  average  number  of  level  crossings  by  the  envelope  per  second 
would  be  equal  to  an  appropriate  bandwidth  factor  multiplied  by  the  enve¬ 
lope  probability  density  function,  it  is  once  again  clear  that  the  gen¬ 
eralized  "t"  model  is  far  more  appropriate  than  a  gaussian  model  for  vlf 
atmospheric  noise.  In  fact,  it  is  seen  that  the  generalized  "t"  model 
gives  the  asymptotic  result 
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whereas 


1  im 

V  OO 
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= s" 


1/2 


(2.51) 


In  addition  to  being  in  good  agreement  with  the  experimental  data  from 
Watt  and  Maxwell  [Ref.  9],  the  asymptotic  result  given  by  Eq.  (2.5l)  is 
also  in  good  agreement  with  the  vlf  level-crossing  rate  measurements 
made  by  Clarke,  e_t  a_l  [Ref.  10],  who  find  for  thresholds  in  excess  of  the 
average  value  of  the  noise  envelope  that  )  is  given  by 

NV(VQ)  «  (Vq/a)'C  ,  (2.52) 

where  A  and  C  are  constants  such  that 

A  «  V  +  30  db  ,  1 . 3  <  C  <  2  . 

avg  -  - 

Taking  now  a  closer  look  at  the  results  plotted  in  Fig.  15,  it 
is  noted  that,  although  the  measured  data  correspond  well  with  the  general¬ 
ized  "t"  model  with  0=2,  the  experimental  determination  of  Pv(Vq) 
contains  no  data  points  for  envelone  values  larger  than  that  corresponding 
to  2.4  crossings  per  second.  Now,  consistent  with  the  arguments  presented 
earlier  in  the  discussion  of  the  distribution  of  the  envelope,  it  is  seen 
that  the  slope  of  Pv(Vq)  must  change  at  high  values  of  VQ,  since  the 
noise  has  finite  average  power.  Thus  it  is  conjectured,  although  it  will 
not  be  considered  further  here,  that  the  fit  of  the  generalized  "t" 
model  to  the  measured  level-crossing  data  can  be  further  improved  by  em¬ 
ploying  the  modified  generalized  "t"  model  proposed  in  subsection  1c. 

In  support  of  this  conjecture,  it  is  noted  that  the  measured  data  presented 

~c 

by  Clarke,  £t^  a_l^  also  shows  a  departure  from  the  power  law  (^q/A)  at 
crossing  rates  below  1  per  second,  with  this  departure  becoming  less  ap¬ 
parent  as  the  value  of  C  increases.  This  behavior  is  consistent  with 
the  result  found  from  envelope  distribution  considerations — namely  that 
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the  generalized  "t"  model  provides  a  better  fit  to  the  received  noise 
at  vlf  as  the  required  value  of  0  increases  corresponding  to  the  situa¬ 
tion  of  low-to-moderate  local  thunderstorm  activity. 

Finally,  in  summary,  it  is  concluded  that  the  data  plotted 
from  Watt  and  Maxwell  in  Fig,  15,  as  well  as  comparisons  with  the  experi¬ 
mental  results  of  Clarke,  et  aJL  indicate  good  agreement  between  the  aver¬ 
age  number  of  level  crossings  of  the  envelope  given  by  the  generalized 
"t"  model  and  by  the  measured  vlf  data.  The  significance  of  this  agree¬ 
ment  is,  of  course,  that  it  is  further  verification  of  the  applicability 
of  the  generalized  "t"  model  as  far  as  the  first-order  statistics  of  the 
noise  are  concerned.  From  a  physical  point  of  view,  this  agreement  in¬ 
creases  confidence  in  the  generalized  "t"  model,  since  it  indicates  that 
on  the  average  the  envelope  of  the  model  varies  with  time  at  the  same  rate 
(at  least  at  vlf)  as  the  envelope  of  the  received  noise.  To  complete 
verification  of  the  applicability  of  the  generalized  "t"  model  as  a  model 
for  vlf  atmospheric  noise,  it  remains  to  investigate  the  manner  in  which 
these  variations  with  time  occur.  This  is,  of  course,  a  function  of  the 
second-  and  higher-order  statistics  of  the  noise,  which  will  be  the  next 
topic  of  discussion.  Before  leaving  the  present  discussion  of  the  average 
level-crossing  rate,  however,  it  is  noted  that  the  fit  of  the  model  to  the 

measured  data  in  Fig.  15  was  achieved  by  taking  the  rms  bandwidth  B  to 

c 

be  500  cps.  On  the  other  hand,  the  measured  data  were  taken  in  a  6-db 
bandwidth  of  1100  cps,  which  corresponds  to  an  rms  bandwidth  of  the  order 
of  400  cps.  This  discrepancy  is  not  large,  but  may  be  a  result  of  the  as¬ 
sumption  that  the  term  ja(t)jE(t)  is  negligible  in  Eq.  (2.38).  Thus  it 
is  conjectured  that  the  fit  of  the  calculated  level-crossing  rate  to  the 
experimental  results  may  be  improved  by  relaxing  the  above  assumption, 
which  is  equivalent  to  the  assumption  that  the  modulating  process  a(t) 
has  a  much  smaller  bandwidth  than  the  envelope  E(t)  of  the  narrowband 
gaussian  process  n(t).  Finally,  it  is  pointed  out  that  the  bandwidth 
factor  is  absorbed  into  the  constant  A  in  the  work  of  Clarke,  ejt  £l  where 
all  of  the  measurements  of  level-crossing  rate  were  made  in  a  fixed  power 
bandwidth  of  425  cps. 
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3.  The  Distribution  of  Envelope  Level  Crossings 


Having  demonstrated  that  the  first-order  statistics  of  the  gener¬ 
alized  "t"  model  are  in  good  agreement  with  measured  data  on  received 
atmospheric  noise  over  an  interesting  range  of  receiver  bandwidth  and  band 
center  frequency,  completion  of  the  verification  of  the  applicability  of 
the  model  to  received  atmospheric  noise  requires  investigation  into  the 
higher-order  statistics  of  the  noise.  Physically,  this  means  that  although 
the  consideration  of  the  average  rate  of  level  crossings  shows  that  the  en¬ 
velope  of  the  generalized  "t"  model  fluctuates  at  the  same  average  rate 
as  the  envelope  of  observed  atmospheric  noise,  it  remains  to  investigate 
the  fashion  in  which  these  fluctuations  occur.  That  is,  il  remains  to 
verify  that  the  higher-order  statistics  of  the  generalized  "t"  model  can 
be  chosen  so  that  the  relationship  between  the  process  at  various  distinct 
instants  in  time,  a3  predicted  by  the  model,  is  consistent  with  measured 
results  on  this  relationship.  As  mentioned  earlier,  the  available  experi¬ 
mental  data  [Ref.  9]  dependent  on  the  higher-order  statistics  of  the  noise 
consist  of  measurements  of  the  probability  distribution  function  of  the 
interval  between  successive  crossings  of  a  fixed  level  by  the  envelope  of 
the  noise.  Inspection  of  these  data  indicates  tl  t  at  vlf  the  noise  pulses 
do  not  in  general  occur  in  a  Poisson  fashion,  but  rather  that  received 
noise  pulses  are  usually  statistically  dependent  on  preceding  ones.  Now, 
the  calculation  of  the  probability  distribution  of  the  inter-level-crossing 
interval  is  difficult,  requiring  numerical  techniques  to  obtain  even  an 
approximate  solution  [Refs.  24,  28,  29].  This  is  true  even  for  the  special 
cast  of  gaussian  processes,  which  is  the  only  case  that  has  been  treated 
in  any  detail  in  the  literature  [Refs.  28,  29j.  Thus  we  will  not  attempt 
to  find  an  exact  solution  for  the  nongaussian  situation  of  interest  here, 
but  will  resort  to  simplifying  assumptions,  based  on  an  understanding  of 
the  generalized  "t"  model,  in  order  to  obtain  ar.  approximate  solution. 
Noting,  as  discussed  in  the  introduction,  that  the  received  atmospheric 
noise  of  interest  can  be  considered  to  be  a  narrowband  random  process,  it 
follows  that  its  envelope  can  at  least  be  assumed  to  be  bandlimited.  Thus 
it  is  proposed  that  a  useful  representation  of  the  envelope  in  an  interval 
of  observation  -At  <  t  T  is  given  by 
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(2.53) 


N 

v(t)  ^  Y  V(tt)  (Ati)1/2  qi(t)  , 

i=0 

where  q  ^  ( t; ) ,  i=0,...,  N,  Is  defined  in  Fig.  16  and  where  the  value  of 
N  required  to  make  this  representation  reasonable  depends  upon  the  band¬ 
width  through  which  the  signal  of  duration  T  +  At  is  observed. 


FIG.  16.  DEFINITION  OF  q± (t) ,  i  =  0,...,  N,  WHERE  t.  =  (i  -  At. 

Now,  recalling  that  with  the  reasonable  assumption  of  negli¬ 
gible  spectral  overlap  of  a(t)  and  n(t)  that 

V(t)  *  { a ( t ) j  E(t)  ,  (2.4) 

where  E(t.)  is  the  envelope  of  the  narrowband  gaussian  process  n(t), 
the  representation  in  Eq.  (2.53)  becomes 

N 

y(t)  'S  j  a(t .  ) j  E(t.)  (At  )1//2  q  (t)  .  (2.54) 

i-0 

Let  T  now  be  defined  as  the  interval  between  a  down-crossing  of  the  level 

V  by  the  envelope  of  the  received  noise,  and  the  next  up-crossing  of  the 
0 
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same  level.  (This  is  the  quantity  whose  statistics  have  been  measured  by 
Watt  and  Maxwell  (Ref.  9].)  Applying  the  representation  of  V(t)  pro¬ 
posed  in  £q .  (2.53),  it  is  seen  that  the  probability  distribution  actually 
measured  by  Watt  and  Maxwell  is  given  by 

P0(T0)  -  Probability  that  T  exceeds 
=  Pr  |T  >  T0| 

=  Prjno  up-crossings  of  VQ  in  [t,  t+T^j  j  down-crossing  of  VQ  at  t| 

Pr  | no  up-crossings  of  in  £o , TQ J ,  down-crossing  of  VQ  at  t  =  0} 

Pr jdown-crossing  of  at  t  =  0| 

Pr(v(t0)>v0,  V(.1)<V0 . V(V<V°}  (2.55) 

=  Prlv'to^  >  vo'  v(t,,<voi 

where  Nq  =  smallest  integer  >  TQ/At . 

With  this  simplified  formulation  of  the  solution  for  Pq(Tq)>  one 
can  now  proceed  to  investigate  the  possibility  of  specifying  the  higher- 
order  statistics  of  a(t)  and  n(t)  in  such  a  way  as  to  reproduce  the 
distribution  of  the  inter-level-crossing  interval  observed  experimentally. 
Two  special  cases  will  be  considered,  as  follows. 

a.  Case  of  Independent  Samples 

In  this  special  case,  which  is  the  simplest  case  of  interest, 
it  is  assumed  that  the  values  taken  on  by  the  envelope  at  the  (N^+l) 
sampling  instants  are  statistically  independent.  It  is  obvious  that  thir 
assumption  of  independence  greatly  simplifies  the  evaluation  of  Pq(Tq) 
given  by  Eq.  (2.55);  but  it  remains  to  verify  that  the  representation  of 
v(t)  given  by  Eq.  (2.53)  is  reasonable  in  the  light  of  the  independence 
assumption.  Assuming  that,  the  observed  atmospheric  noise  is  bandlimited 
by  the  receiver  to  an  rf  bandwidth  2B,  it  follows  that  the  envelope  of  the 
received  noise  is  bandlimited  to  a  frequency  band  of  width  23  centered 
about  zero  frequency  [Ref.  30 ].  Thus  it  can  be  shown  [Ref.  31 ]  that  this 


SEL-60-052 


-  64 


bandlimitcd  envelope  v(t)  is  approximately  described  in  the  time  inter¬ 
val  [O.Tq]  by  its  samples  at  equidistant  instants,  where  r.  2T^B 

is  required  if  the  samples  are  to  be  reasonably  assumed  both  independent 
and  sufficient  in  number  to  yield  a  good  approximation.  Such  a  represen¬ 
tation  is  given  by: 


V(t.)  ?.(t)  , 


(2.56) 


where  the  are  defined  as  in  Fig.  16  and  where  the  best  approximation 

to  V(t)  comes  from  taking  the  (t)  to  be  the  prolate  spheroidal  wave 
functions  [Ref.  31 ].  However,  because  the  spectral  shape  of  the  received 
atmospheric  noise  is  determined  by  the  filter  through  which  it  is  observed 
in  the  situations  of  practical  interest,  it  is  proposed  that  the  represen¬ 
tation  in  Eq.  (2.53)  in  terms  of  the  pulse  "bases''  *  lions  q  (t), 

i=0,  ...,  2TqB,  is  in  fact  a  reasonable  approximation  to  the  envelope  of 
this  filtered  "white"  noise. 

Now,  making  use  of  the  assumption  that  the  (2TQB+1 )  samples 
in  Eq.  (2.55)  are  statistically  independent,  it  is  found  directly  that 


r  -i  (2T  B-l ) 

po(V  *  [Pr  iv(*)<  v0|J  “ 


(2.57) 


where,  from  Eq.  (2.27), 


Pr|v(t)  <  V'0|  =  1  - 


/,,2  2 

\V0  +  ?  . 


(2.58) 


Thus,  for  the  case  of  independent  samples, 
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O'  0' 


r 

/  \9_i  ^ 
[JL\ 

Vo) 

P 

i  i 

2-(e-iJ/2> 

/JL\1 

l  + 

(v  J  j 

(2TQB-l) 


(2.59) 


Special  cases  of  this  result  corresponding  to  0=2  and  0-3  are 
plotted  in  Fig.  17  along  with  data  taken  from  Watt  and  Maxwell  [Ref.  9], 
and  the  special  case  corresponding  to  @  =  2  is  also  plotted  in  Fig.  18 
along  with  data  taken  from  Watt  and  Maxwell  and  along  with  data  from  the 
constant  a(t)  case  considered  below.  It  is  seen  that  the  assumption  of 
independent  samples,  like  the  assumption  that  noise  pulses  occur  in  a 
Poisson  fashion  (exponentially  distributed  intervals,  Fig.  17),  neglects 
certain  dependencies  between  adjacent  level  crossings.  In  fact,  it  is 
seen  that  for  the  high  levels  of  most  interest,  i.e.,  for  VQ  »  y, 


lim  P0(Tq)  =  exp j“  2B^  TQ  , 

7 


(2.60) 


so  that  the  assumption  of  independent  samples  produces  exponentially  dis¬ 
tributed  inter-level-crossing  intervals  at  high  levels.  It  is  noted  that 
the  value  of  Pq^To^  Biven  by  Eq .  (2,59)  and  plotted  in  Figs.  17  and  18 
is  actually  consistent  with  Eq .  (2.55)  only  for  those  values  of  TQ  such 
that  2TqB  is  an  integer.  A  smooth  curve  is  fitted  through  these  points 
in  preference  to  the  step  function  specified  by  Eq.  (2.55)  because  the 
probability  distribution  function  of  the  interval  is  of  course  continuous 
in  nature.  In  any  case,  the  discrepancy  between  the  smooth  curve  and  the 
associated  step  function  is  negligible  at  large  values  of  2TqB,  the  sit¬ 
uation  of  most  interest.  Now,  the  goal  of  the  calculation  of  the  probabil¬ 
ity  distribution  of  the  inter-level-crossing  interval  here  is  to  investigate 

^It  is  noted  that  the  discrepancy  between  the  calculated  and  measured  re¬ 
sults  at  low  values  of  T0  is  due  in  part  to  the  breakdown  of  the  assump¬ 
tion  that  Nq  =  2TQB  is  sufficient  at  low  values  of  the  TqB  product. 
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our  ability  to  specify  the  higher-order  statistics  of  the  generalized 
"t"  model  in  order  to  reproduce  the  distribution  of  inter-level-crossing 
intervals  observed  experimentally.  The  assumption  of  independent  samples 
is  seen  to  give  a  distribution  which  agrees  with  some  experimental  data 
but  which,  in  general,  neglects  some  observed  dependence  between  adjacent 
level  crossings.  A  second  special  case  will  now  be  considered  in  an 
attempt  to  introduce  this  observed  dependence  into  the  calculated  results. 

b.  Constant  aft)  Case 

In  this  special  case  the  slowly  varying  modulating  function 
a(t)  is  in  fact  assumed  to  be  constant  for  the  duration  of  observation, 
so  that  v(t)  as  represented  by  Eq.  (2.54)  takes  the  form 


N 

V(t)  *=  ^  |a|  E(t.)  (At.)1/2  qi(t)  .  (2.6l) 

i=0 

In  order  to  complete  the  specification  of  V(t)  in  the  interval  L-At,TQJ, 
which  is  the  interval  of  interest  in  Eq,  (2.55),  it  is  assumed  that  the 
(Nq+i)  samples  of  E(t)  are  statistically  independent  [recalling  that 
E(t)  is  the  envelope  of  the  narrowband  gaussian  process  n(t)].  Thus  it 
follows  from  the  development  of  the  independent-samples  case  that 
Nq  =»  2TqB  is  required,  so  that  the  representation  of  v(t)  in  the  inter¬ 
val  (-At,  Tq]  takes  in  this  case  the  form 

2ToB 

V(t )  as  j  a |  ^  E(ti)  (Ati)1/2  q.(t)  .  (2.62) 

i=0 


Now,  turning  to  the  calculation  of  Pq(Tq)  as  given  by  Eq.  (2.55),  it  is 
seen  that  the  solution  requires  evaluation  of  the  joint  probability  dis¬ 
tribution 


v(tx)  <  V 


0' 


v(t„)  <  V0}  $  PS(V0)  .  (2.63) 
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INTERVAL 


FIG.  17.  PROBABILITY  DISTRIBUTION  OF  THE  TIME  INTERVAL  BETWEEN  ENVELOPE  LEVEL  CROSSINGS*.  COMPARISON 
OF  MODEL  RESULTS  WITH  MEASURED  VLF  DATA. 
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TIME  INTERVAL,  Tn  (ms) 


FIG.  18.  PROBABILITY  DISTRIBUTION  OF  THE  TIME  INTERVAL  BETWEEN  ENVELOPE  LEVEL  CROSSINGS:  COMPARISON 
OF  MODEL  RESULTS  WITH  MEASURED  VLF  DATA. 
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Now,  in  this  constant  a(t) 


case 


in  which 


Thus 


where 


v(ti) 

.  -  , 

Vi| 

Eil 

Py(V)  =  |  dx  p(a(  (x)  P£(f)  *  (2 . 65 ) 


follows  from  Eq,  (2.10),  since  E(t)  is  the  envelope  of  n(t),  and  where 
p,  1 (x)  is  given  by  Eq.  (2.4l).  Evaluating  Eq.  (2.65)  by  making  a  change 

I  al 

of  variable  so  that  Dwight,  item  860.17,  is  appropriate: 


S EL-66-052  -  72  - 


PV(X> 


--^  2 


N+l 


r/2N  +  m  -t  2j 

r(f) 


m 


N 

n  (v.) 

i=0 


y  7 


N 


-,(2N+m+2)/2 


(2.66) 


2ti) 


Li-0 


Substituting  this  result  into  Eq.  (2.64), 


W  - 2 


/2N  +  m  +  2 1 


m 

7 


(2.67) 

Now,  Eq.  (2.67)  can  be  evaluated  by  recursively  applying  Dwight,  item 
201,9,  which  gives  the  result 


/2N  +  m  +  2 ' 


W 


=  2 


N+l 


'(f) 


7  (2N  +  ra)(2N  +  m  -  2$  ...  (mj 


'I 


<-»k  C) 

1 

[(k+1  )Vq  +  72] 

m/2 

(2.68) 


where 


is  the  binomial  coefficient. 
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Finally,  substituting  Eq.  (2.68)  with  the  appropriate  value 
for  N  into  the  numerator  and  denominator  respectively  of  Eq.  (2.55) 
gives  the  desired  result 


po(To>  - 


4 


« 


•  4 

*  *' 

;  i 

i 

! 


(2.69) 


where  NQ  =  2TqB.  In  the  interesting  special  case  0  =  3  (m  =  2),  this 
gives 

2TB  ,  / 2T  B\ 

y  °) 

k“°  *  &) 

po(To)  1 3  ■  -  —  - 0  •  <2-™> 

9  y  (-i) 

k=0  k.l  +  fa) 

This  result  is  plotted  in  Fig.  18  for  various  values  of  -y/v  ,  along  with 

0 

results  from  the  independent-samples  case  and  experimental  data  taken  from 
Watt  and  Maxwell.  These  plotted  results  indicate  that  the  constant  a(t) 
case  is  characterized  by  the  fact  that  very  short^  and  very  long  inter¬ 
level-crossing  intervals  occur  with  a  higher  probability  at  high  levels 


This  conclusion  follows  when  we  neglect  Tq  corresponding  to  2TqB  <  1, 
where  the  Nq  -  2TqB  assumption  breaks  down. 
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than  is  actually  observed  experimentally.  This  means  that  the  constant 
a(t)  case  predicts  more  "clustering"  of  noise  pulses  than  is  consistent 
with  measured  data,  as  contrasted  to  the  independent-samples  case,  which 
predicts  less  "clustering"  than  is  actually  observed.  The  usefulness  of 
this  result  is  that  it  indicates  the  versatility  of  the  generalized  "t" 
model,  and  it  is  conjectured  that  the  observed  distribution  of  inter-level- 
crossing  intervals  can  be  obtaineo  with  the  generalized  "t"  model  with 
the  proper  specification  of  the  higher-order  statistics  of  the  modulating 
process  a(t).  Although  it  will  not  be  further  considered  here,  it  is 
suggested  that  a  reasonable  approach  to  this  problem  follows  from  con¬ 
sidering  che  special  case  9=2  for  which  the  generalized  "t"  model 
takes  the  form 


Kt)  = 


bin 


n(t )  , 


(2.71) 


where  n(t)  is  a  zero-mean  narrowband  gaussian  process  with  covariance 
function  R  (t),  and  b(t)  is  a  zero-mean,  slowly  varying  gaussian  pro¬ 
cess,  independent  of  n(t),  with  covariance  function  R^x).  Thus  in 
this  case  the  higher-order  statistics  of  the  modulating  process  are  com¬ 
pletely  specified  by  the  covariance  function  R^(x),  so  that  investiga¬ 
tion  of  the  proper  choice  of  these  statistics  appears  to  be  analytically 
feasible. 


D.  SUMMARY  AND  CONCLUSIONS 

The  work  in  this  chapter  has  been  concerned  with  the  development  of 
an  analytical  model  for  "impulsive"  phenomena  and  with  verification  of  the 
applicability  of  this  model  as  a  description  of  received  atmospheric  noise. 
The  generalized  "t"  model  proposed  here  takes  the  received  atmospheric 
noise  y(t)  to  be  given  by 


y(t)  =  a(t)  n(t )  , 


(2.3) 


where  n(t)  is  a  zero-mean,  narrowband  gaussian  process  with  covariance 

function  R  (x),  and  a(t)  is  a  stationary  slowly  varying  random  process, 
n 
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independent  of  n(t),  which  modulates  n(l).  This  modulating  process 


is  further  described  as 


a  ( t)  = 


(2.72) 


where  the  first-order  statistics  of  b(t)  arc  specified  by  the  "two-sided" 
chi  distribution  with  parameters  m  and  cr  given  by  Eq.  (2.8). 

The  applicability  of  the  generalized  "t"  model  as  a  model  for  re¬ 
ceived  atmospheric  noise  has  been  investigated  in  some  detail  in  this 
chapter,  and  the  pertinent  results  can  be  summarized  as  follows:  The 


first-order  statistics  of  the  generalized 


model  are  in  good  agree¬ 


ment  with  experimental  results,  with  this  agreement  be^ng  particularly 
good  at  vlf  and  If  in  those  situations  characterized  by  low-to-moderate 
local  thunderstorm  activity.  This  conclusion  follows  from  results  ob¬ 
tained  in  Section  Cl  and  C2,  where  it  is  shown  that: 

1.  The  probability  distribution  of  the  envelope  of  the  generalized  "l" 
model  is  in  good  agreement  with  a  large  amount  of  measured  data  on 
received  atmospheric  noise.  This  is  demonstrated  in  Figs,  4  to  13. 

It  is  noted  that  the  agreement  with  measured  data  achieved  by  the 
generalized  "t"  model  compares  favorably  with  that  achievable  by 
either  the  filtered  impulse  or  the  empirical  models  discussed  in 
Section  A2. 

2.  The  phase  of  the  generalized  "t"  model  is  distributed  uniformly 
in  the  interval  (o,  2n]  in  agreement  with  intuition. 

3.  The  average  rate  of  crossings  of  a  fixed  level  by  the  envelope  of 
the  generalized  "t"  model  is  in  good  agreement  with  the  available 
measured  data.  This  is  demonstrated  in  Fig.  15,  and  means  physi¬ 
cally  that  the  envelope  of  the  generalized  "t"  model  fluctuates 

at  the  same  average  rate  as  does  the  envelope  of  received  atmospheric 
noise.  It  is  noted  that  no  calculation  of  this  average  rate  of  level 
crossings  has  beer,  reported  for  any  of  the  empirical  models,  although 
Nakai  [Ref.  25 J  has  obtained  numerical  results  in  agreement  with  the 
measured  data,  using  a  f iltered-impulse  model. 

The  higher-order  statistics  of  the  generalized  "t"  model  can  be  specified 
to  give  good  agreement  with  experimental  results.  This  conclusion  follows 
from  the  calculation  in  Section  C3  of  the  probability  distribution  of  the 
interval  between  successive  crossings  of  a  fixed  level  by  the  envelope  of 
the  model.  This  is  demonstrated  in  Fig.  18,  where  the  plotted  data  indi¬ 
cate  that  the  higher-order  statistics  of  the  modulating  process  a(t) 
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[it  is  to  be  noted  that  these  statistics  are  not  yet  specified  in  Eq .  (2.72) 
above]  can  be  specified  to  produce  the  observed  dependence  between  adjacent 
level  crossings.  Physically  this  means  that  the  generalized  "t"  model 
can  bo  specified  to  imitate  the  tendency  of  received  atmospheric  noise  to 
consist  of  clusters  of  noise  pulses.  This  is  an  advantage  of  the  gener¬ 
alized  "t"  model  over  any  of  the  empirical  models  or  any  of  the  filtered- 
impulse  models  which  assume  the  noise  pulses  to  occur  <n  a  Poisson  fashion. 

Finally,  it  is  concluded  that  the  generalized  "t"  model  is  an  appro¬ 
priate  model  for  received  atmospheric  noise  that  may  be  useful  in  the  study 
of  signal  detection  and  estimation  problems  in  the  presence  of  additive 
atmospheric  noise.  This  will  be  demonstrated  in  Chapter  III,  where  the 
detection  problem  is  examined  in  detail. 


e 
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III.  APPLICATION  OF  THE  GENERALIZED  "t"  MODEL 


A.  INTRODUCTION 

In  Chapter  II  a  new  model  for  impulsive  phenomena  was  (leveloDed  a.'d 
was  demonstrated  to  be  applicable  to  the  representation  of  received  atmo¬ 
spheric  noise.  The  work  in  this  chapter  is  concerned  with  the  application 
of  this  model  to  the  signal  detection  problem  in  the  presence  of  additive- 
atmospheric  noise.  Thus*,  the  communication  channel  to  be  considered  here 
is  the  additive-noise  channel  shown  in  Fig,  19,  where  m(t)  is  the  trans¬ 
mitted  signal,  y(t)  is  the  additive  atmospheric  noise,  and  x(t)  is  the 
received  sum  of  signal  plus  noise.  Now,  the  goal  of  the  analysis  is  che 
determination  of  the  recei\er  form  which  detects  the  transmitted  signal 
in  the  presence  of  additive  atmospheric  noise  in  an  optimal  manner  with 
respect  to  a  performance  criterion  to  be  specified.  Thus,  while  the  addi¬ 
tive  noise  y(t)  is  given  by  the  generalized  "t"  model,  completion  of 
the  statement  of  the  statistical  detection  problem  requires  specification 
of  both  the  character  of  the  transmitted  signal  and  the  performance  cri¬ 
terion  to  be  used. 

1 .  Summary  of  Known  Results 

Before  proceeding  wJ  .h  this  specification,  however,  a  brief  summary 
of  known  results  on  the  detection  problem  in  the  presence  of  additive  atmo¬ 
spheric  noise  is  presented.  To  the  author's  knowledge,  no  analytical  re¬ 
sult  specifying  an  optimal  receiver  form  for  use  in  the  presence  of  additive 


FIG.  19.  DESCRIPTION  OF  ADDITIVE-ATMOSPHERIC-NOISE 
COMMUNICATION  CHANNEL. 


atmospheric  noise  is  available  in  the  literature.  Rather,  the  receiving 
techniques  that  have  been  proposed  stem  from  observation  of  the  distinc¬ 
tive  impulsive  nature  of  the  noise  and  can  be  conveniently  classified  in 
two  categories: 

1.  Linear  receivers.  At  least  two  linear  receiving  techniques  have 
been  proposed [Refs.  17,  32]  for  use  in  the  presence  of  atmospheric 
noise.  One  of  these  is  the  "smear-desmear"  technique  proposed  by 
several  workers  and  discussed  by  Lerner  [Ref.  17).  This  technique 
consists  of  linearly  smearing  the  received  noise  pulses  in  time,  so 
that  the  smeared  noise  can  be  modeled  as  a  gaussian  process.  Having 
done  this,  the  optimal  receiver  in  the  presence  of  the  smeared  noise 
is  just  the  well-known  matched-filter  receiver  which  must,  of  course, 
be  matched  to  the  smeared  signal.  It  is  clear  that  this  technique 

lc  suboptimal  in  the  presence  of  atmospheric  noise.  This  can  be  seen 
quantitatively  by  comparing  the  error  curve  resulting  from  its  use 
with  the  optimal  error  curve,  if  known.  In  fact,  for  the  known  sig¬ 
nal  situation  considered  in  this  chapter,  it  will  be  seen  that  use 
of  the  smear-desmear  technique  corresponds  to  a  loss  of  12-15  db  in 
input  SNR. 

2.  Nonlinear  receivers.  Most  proposed  techniques  [Refs.  17,  33,  34], 
and  in  fact  the  techniques  generally  used  in  practice,  employ  non¬ 
linear  processing  of  the  received  signal.  The  purpose  of  this  non¬ 
linear  processing  is  to  reduce  the  effects  of  noise  pulses  on  the 
receiver  decision  by  amplifying  the  received  .'ignal  nonlinearly  prior 
to  the  performance  of  any  linear  filtering.  (This  of  course  refers 
to  matched  filtering,  and  not  to  the  bandpass  filtering  discussed  in 
Section  IIB,  which  is  required  in  practice  in  the  front  end  of  the 
receiver  to  suppress  adjacent  channel  interference.)  It  will  be 
shown  that  this  procedure  is  supported  by  the  "logarithmic-correlator" 
receiver  shown  to  be  optimum  in  the  known  signal  situation  analyzed 

in  this  chapter.  It  is  noted  that  the  log-correlator  receiver  is 
compared  in  detail  with  specific  nonlinear  receivers  used  in  practice 
in  Chapter  IV. 

Although  no  analytical  determination  of  optimal  receiver  forms  is  available 
in  the  literature,  there  are  several  analytical  results  of  interest  on  the 
performance  of  specific  receiving  techniques.  Several  workers  [Refs.  35, 

36,  37,  38]  have  computed  the  probability  of  error  resulting  from  the  use 
of  a  matched  filter  receiver  in  the  presence  of  additive  atmospheric  noise. 
Bello  [Ref.  37]  and  Conda  [Ref.  38]  have,  in  addition,  included  the  effects 
of  fading  on  these  error  curves.  It  will  be  seen  in  Section  D4  below  that 
the  results  obtained  by  these  investigators  in  the  absence  of  fading,  using 
various  models  for  the  atmospheric  noise,  are  in  good  agreement  with  the 
results  obtained  using  the  generalized  "t"  model.  These  results  lead  to 
the  conclusion  that  the  matched  filter  is  a  poor  receiver  choice  in  the 
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presence  of  additive  atmospheric  noise.  Finally,  we  note  the  work  of 
Bowen  (Ref.  7 J  who  uses  an  SNR  criterion  to  investigate  the  use  of  hard 
limiting  in  the  presence  of  atmospheric  noise.  He  shows  that  the  use  of 
a  hard  limiter  followed  by  a  zonal  filter  produces  a  much  greater  improve¬ 
ment  in  the  ratio  of  output  SNR  to  input  SNR  when  the  additive  noise  is 
impulsive,  than  is  obtained  when  the  additive  noise  is  a  gaussian  process. 
This  result  supports  the  use  of  nonlinear  processing  of  the  received  sig¬ 
nal  as  mentioned  above. 

With  this  brief  summary  of  available  results  in  mind,  the  specifi¬ 
cation  of  the  detection  problem  to  be  considered  here  will  now  be  completed. 
There  are,  of  course,  several  possible  specifications  of  the  transmitted 
signal  of  practical  importance.  Consideration  will  be  focused  here  on  the 
simple  hypothesis-testing  situation,  wherein  the  transmitted  signal  is  one 
or  the  other  of  two  a  priori  equally  probable  known  signals,  This  sit¬ 
uation  is  described  in  Fig.  20,  where  m^^(t)  and  m^2^(t)  are  known 
signals  described  bv 

m^(t)  =  m(t)  ,  0  <  t  <  T  , 

m^(t)  =  0  ,  (3.1) 


with 


dt  =  E 


and  where  y(t)  is  the  additive  atmospheric  noise  described  by  the  gener¬ 
alized  "t"  model.  The  performance  criterion  to  be  used  here  is  the 
probability-of -error  criterion,  so  that  the  statistical  detection  problem 
at  the  receiver  becomes  that  of  choosing  between  the  £  priori  equally 
probable  hypotheses 


h^1^:  x(t)  =  m(t)  +  y(t) 

> 

h^:  x(t)  =  y(t) 


0  <  t  <  T 


(3.2) 
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FIG.  20.  DESCRIPTION  OF  "KNOWN ''-SIGNALS  DETECTION  PROBLEM. 

and  nr  '(t)  are  equally  probable,  a  priori. 

with  the  smallest  probability  of  error.  Now,  it  is  well  known  [Refs.  5, 
39,  40,  41 ]  that  the  decision  rule  which  minimizes  the  probability  of 
error  in  this  situation  is  the  Bayes  rule  given  by 

Choose  h^  iff 


L[x(t)] 


£  >  i 

~  w[x(t)|h(2)]  - 


(3.3) 


where  w[x(t)j  h^  ^  ] ,  i=l,  2,  is  the  probability  density  of  x(t), 

0  <  t  <  T,  under  hypothesis  (i),  and  where  L[x(t)]  is  the  likelihood 
ratio.  Thus  the  solution  of  che  detection  problem  requires  the  calculation 
of  w[x(t)|h^^J,  i=l,  2,  which  will  be  done  here  by  expressing  x(t), 

0  <  t  <  T,  in  vector  notation  with  respect  to  a  suitable  set  of  ortho¬ 
normal  basis  functions.  In  order  to  express  x(t)  ir.  this  manner,  how¬ 
ever,  it  is  necessary  to  specify  the  temporal  behavior  of  the  additive 
atmospheric  noise;  i.e,,  it  is  necessary  to  specify  the  higher-order  sta¬ 
tistics  of  the  generalized  "t"  model. 

It  is  recalled  that  the  generalized  "t"  model  takes  the  received 
atmospheric  noise  to  be  given  by 
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y ( t )  =  a(t)  n(t)  , 


(2.3) 


where  n(t)  is  a  zero-mean  narrowband  gnussian  process  with  covariance 
Rn(T),  and  a(t)  is  a  stationary  slowly  varying  random  process,  indepen¬ 
dent  of  n(t),  whose  first-order  statistics  are  given  by  Eq,  (2.72), 

Thus,  as  in  Section  IIC3,  it  is  convenient  to  distinguish  two  cases,  de¬ 
pending  on  the  temporal  behavior  of  the  modulating  process  a(t). 


B.  CASE  I:  SHORT -DURATION  SIGNALS,  a(t)  =  a,  0  <  t  <  T 

This  case  is  similar  to  the  constant  a(t)  case  considered  in  Section 
IIC3,  in  that  it  assumes  the  slowly  varying  modulating  process  a(t)  to 
remain  constant  for  the  signal  duration  T.  Physically,  this  corresponds 
to  the  situation  in  whieh  the  signal  duration  is  sufficiently  short  that, 
with  high  probability,  the  modulating  process  a(t)  can  be  assumed  con¬ 
stant  in  the  signaling  interval.  (The  practical  applicability  of  this 
assumption  will  be  discussed  further  in  Section  B5.)  In  this  case  the 
decision  problem  reduces  to  choosing  between  the  two  hypotheses 


h^:  x(t)  =  m(t)  +  a  n(t) 
hv2‘\  x(t)  -  a  n(t) 


0  <  t  <  T  . 


(3.4) 


An  often-used  procedure  for  developing  the  desired  vector  formulation, 

which  is  applicable  to  this  problem,  is  to  expand  the  gaussian  process 

n(t)  in  terms  of  a  Karhunen-Lolve  expansion.  Assuming  that  R  (t)  is 

n 

continuous  in  T,  it  is  possible  to  write  [Refs.  30,  41  ] 


n(t)  =  l.i.m.  'S'  n  <p  (t)  ,  0  <  t  <  ^  » 

v  ^  ..  *  * 


N  oo 


i=l 


where  l.i.m.  denotes  "limit  in  the  mean," 

T 


n.  =  J  n(t)  cpj,  ( t )  dt  , 


(3.5) 


(3.6) 
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and  the  »t(t)  arc  the  eigenfunctions  of  «n(«,s),  1...,  they  are  the 

solutions  to 


JL 

j  Rn(t,s)  ^(s)  ds  =  9A(t)  ,  0  <  t  <  T 


(3.7) 


This  expansion  has  the  desirable  properties: 

1.  The  cp^t)  are  orthogonal,  i.e.,  they  can  be  normalized  to  satisfy 


X 

J  <^(0  (Pj(t)  dt  =  8ij  * 


(3.8) 


where  6^  is  the  Kronecker  delta. 
2.  The  nA  are  uncorrelated,  i.e., 


Et”i  "j1  -  5ij  • 


(3.9) 


If  R^(t,s)  is  strictly  positive  definite,  the  cp^t),  1=1,  ...,  », 
will  span  L2(t)  t^ie  (Hiltjert)  space  of  square  integrable  functions  or 
[o,T].  Then  the  signal  m(t)  can  be  expressed  in  terms  of  the  p^(t) 


as 


m(t)  =  l.i.m.  }  m  q>.  (t)  ,  0  <  t  <  T  ,  (3.10) 

N  -  co  ^  1 


i=l 


If  Rn(t,s)  is  not  strictly  positive  definite,  it  suffices  to  assume 
that  m(t)  lies  in  the  space  spanned  by  the  eigenfunctions  of  Rn(t,s). 
If  this  is  not  true,  error-free  detection  can  obviously  be  obtained. 
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where 


m. 


(3.11) 


Having  expressed  n(t)  and  m(t)  in  terms  of  the  orthonormsl  functions 

9  (t),  i=l . N,  in  Eqs.  (3.5)  and  (3.10),  the  problem  can  now  be 

formulated  in  vector  notation  as  follows: 


h 


x  =  m  +  a  n  , 


h 


x 


a  n  , 


(3.12) 


where,  for  example,  n  =  n  ,  ...,  n  is  an  element  of  Euclidean  N-space 

1  i_i _ Nj 

E  .  It  is  noted  that  this  formulation  is  written  in  terms  of  N  dimen- 
sions,  so  that  N  -*  oo  will  be  taken  at  the  conclusion  of  the  problem. 

For  a  rigorous  justification  of  this  procedure,  see  Grenander  [Ref.  42 j. 

1  •  Calculation  of  the  Likelihood  Ratio 

One  can  write 


00 

w(x|h^)  =  J  w(x|h^,  b)  Pb(b)  db  , 

-00 


(3.13) 


where  p.  (b)  is  given  by  Eq,  (2.7),  and  where  Eq.  (3,12)  plus  the  fact 

k  A 

that  a(t)  =  l/b(t)  gives 


w(2ilh^>  b)  =  lb|N  Pn(b(i  "  » 


(3.14) 


denotes  the  transpose  of  the  column  vector 


n . 
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. . .  vi  w»  n  u 


Similarly, 


w(x|h<2)) 


rjrjz)  y) 

r(|)  ,"/2  |A| 


Hill2.! *  % 

La  or  J 


[m+N  )/2  ' 


m  >  0  . 


(3.19) 


Now,  substituting  Eqs.  (3.18)  and  (3.19)  into  the  vector  formulation  of 
Eq.  (3.3),  there  results 


r  9  m  1 {m 

ll-H  _!  +  ■% 

La  ct  J 

lli'iall2-!  + 

A  cr  J 


(m+N )/2 


.m+N)/2 


(3.20) 


Thus,  the  optimal  decision  rule  is  given  by: 
Choose  h^  iff 


-i|2-l  +  ^2  -  fe‘^l|2-l  +  ^2  » 

Act  Act 


(3.21) 


which  reduces  to  the  familiar  rule  for  the  detection  of  known  signals  in 
additive  gaussian  noise: 


Choose  h^1^  iff 


<  2.  »  >  _!  >  £  IIhII2-! 

A  1  1  A 


(3.22) 


where 


^  ^  A  -1 

<  x,  m  >  _  =  x  A  m 
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Now,  letting  N  -*  oo,  it  is  seen  that  in  terms  of  the  time  functions 
available  to  the  receiver,  rule  Eq.  (3.22)  goes  over  into 
( l ' 

Choose  h‘  1  iff 


A  1 

|  x(t)  c(t)  at  >\j  m(t)  c(t)  dt  , 


(3.23) 


where  c(t)  is  the  solution  to  the  Fredholm  integral  equation 


/  R„(t's 


s)  c(s)  ds  =  m ( t )  , 


0  <  t  <  T 


(3.24) 


2 .  Discussion  of  the  Optimal  Receiver  Rule 

It  is  seen  from  Eq.  (3.23)  that  the  receiver  that  minimizes  the 
probability  of  error  for  the  additive  generalized  "t"  noise  channel  in 
the  case  where  a(t)  =  a  is  constant  for  the  duration  of  the  signal  is 
the  well-known  correlator  or  matched-filter  receiver.  Before  proceeding 
to  the  calculation  of  the  probability  of  error  resulting  from  the  use  of 
this  optimal  receiver,  there  are  several  points  of  interest  to  note.  First, 
the  decision  rule  is  independent  of  the  parameters  m  and  <j  of  the  gen¬ 
eralized  "t"  model;  i.e.,  the  same  receiver  is  applicable  in  the  presence 
of  additive  generalized  "t"  noise  regardless  of  the  value  of  the  param¬ 
eter  0  >  1.  (it  will  be  shown  in  Section  B3  that  similar  statements  are 
also  true  for  additive  modified  generalized  "t"  noise.)  Also,  it  is 
worthwhile  to  note  that  the  receiver  rule  Eq .  (3,23)  could  have  been  de¬ 
duced  directly  by  inspection  of  p  (j^),  the  probability  density  of  the 
additive  noise  vector.  This  deduction  proceeds  as  follows:  Noting  that 
w(x|n  )  given  by  Eq.  (3.19)  is  just  p^(x),  there  results 


pi(Jd 


?)  _ , 

r(f) 


m+N )/2  ’ 


m  >  0 


(3.25) 
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where  y  =  y, ,  . . . ,  y„  .  However,  assuming  R  (t,s)  to  be  strictly 
t  i  1  N  i  n 

positive  definite,  there  exists  n  nonsingulur  "whitening"  transformation 
W  such  that 


m  >  0  , 


(3.26) 


where  y  =  Wjr.  In  fact,  in  this  case  it  is  easily  seen  that 


W  = 


(3.27) 


Now,  in  terms  of  the  whitened  noise  vector  the  detection  problem 

becomes  that  of  choosing  between  the  £  priori  equally  probable  hypotheses 


(1) 


x  =  m  +  y 
— w  — w 


(2) 

I  *  • 


*w  =  *w  • 


(3.28) 


It  can  easily  be  shown  that  the  correlator  receiver  given  by  Eq.  (3.23) 
is  indeed  the  optimal  receiver  in  this  case,  and  in  fact  this  result  can 
be  generalized  as  specified  in  the  following  well-known  proposition: 


Proposition  3.1.  Given  an  additive  noise  channel  where  the  two  hypotheses 
can  be  written  in  the  form  of  Eq.  (3.28),  and  given  that  the  probabil¬ 
ity  density  function  of  the  additive  noise  can  be  written  as  a  monotoni- 
cally  decreasing  function  of  the  norm  of  the  noise  vector,  then  the 
receiver  that  minimizes  the  probability  of  error  is  a  correlator  re¬ 
ceiver  . 
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Proof  of  Proposition  3.1. 

The  proof  of  this  proposition  is  so  short  that  it  will  be  given  here. 
The  likelihood  ratio  is  given  by 


w(x  |h^) 


w(x  ! hv  ) 
— w' 


P  (x 
Jw_7 
~p  (x  ) 


-  m  ) 
— w 


f(ll\  -  aj 

~ fTiijFT 


(3.29) 


where  iU|x'i)  is  a  mono  ton  leal  ly  decreasing  function  of  jjxj|.  Therefore, 
from  Eq.  (3.3),  the  optimal  decision  rule  is  given  by: 

Choose  h^1^  iff 


f  ( 


|x  -  m  jj)  >  f(||x  i|) 
— w  — w  —  — w 


<  =  >  Choose 


iff 


m  I*  < 
— ■ w'  ~ 


!x  ; 

— w 


t 


(•  .30) 


(3.31) 


since  f(  |x||)  is  r  monotonically  decreasing  function  of  j|xjj. 
<  =  >  Choose  h^1 ^  iff 


< 


X 

— w 


m 

— w 


> 


This  completes  the  Proof  of  Proposition  3.1. 


(3.32) 


3.  Case  Im:  Modified  Generalized  "t"  Noise 

As  an  application  of  Proposition  3.1,  the  results  of  Case  I  will 
now  be  extended  to  the  case  of  additive  modified  generalized  "t"  noise. 
In  this  situation  the  additive  noise  is  given  in  the  vector  notation  of 
Eq.  (3.12)  by 


£  =  an 


(3.33) 
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where 


(K-)  ■  K-* + p)  ■ 


(3.36) 


For  m+N-2  >  0  and  odd,  Eq.  (3.34)  is  evaluated  in  Appendix  A,  where  it 


is  shown  that 


py(l)  =  ^7^!  ’A  |K'  2p2 


4  l»-+  »  -  2) 

pin+N-2  f  ju *  pm+K-4  *  - 


,(m+N-3)/2  (m  + 


)3PK-)j 


m+N-3)/2 


,(m+N-l)/2  (m+N-2)  (m  +  N  -  4) 


KK-4 


2 

(m+N-l)/2 


•••  (i)<o 


K 


<i/z 

_ A  A  i 


(3.37) 


Nov:,  assuming  that  R  (t.s)  is  strictly  positive  definite,  the  whitening 

n 

transformation  W  given  by  Eq.  (3.27)  can  be  applied  to  the  modified  gen¬ 
eralized  "t"  noise  vector  £  to  produce  the  whitened  noise  vector 
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v  =  Wj£.  Furthermore,  inspection  of  p  (jf)  given  by  Eq.  (3.35)  or  Eq. 
(3.37)  shows  that  the  probability  density  p^  (^w)  of  the  whitened  noise 
vector  will  be  a  monotonically  decreasing  function  of  Thus,  ap¬ 

plying  Proposition  3.1,  the  receiver  that  minimizes  the  probability  of 

error  in  the  presence  of  Case  I  modified  generalized  "t"  noise  is  a 

m 

correlator  receiver.  In  fact,  inspection  of  the  proof  of  Proposition  3.1 
shows  that  this  optimal  receiver  is  the  same  correlator  receiver,  speci¬ 
fied  by  Eq.  (3.23),  that  is  optimal  in  the  presence  of  Case  I  generalized 
"t"  noise. 


4.  Calculation  of  the  Probability  of  Error  for  Case  1  and  Case  I 


m 


As  given  by  Eq.  (3.22),  the  receiver  that  minimizes  the  probability 

of  error  in  Case  I  (and  also  in  Case  I  )  is  the  one  that  implements  the 

m 

rule 


Choose 


h^  iff 


<  is  >  >  ^•!®l!2_1 

A  A 


(3.22) 


Since  the  hypotheses  are  £  priori  equally  probable,  the  probability  of 


error  P  is  given  by 
e 


=  Pr|<  x,  m  >  _x  >  ^itmt|2_1 

l  A  A 

“  pr  j<  m  >  _!  >  •|i!iE!|2_1| 


.(2) 


is  true 


(3.38) 


However,  this  expression  is  identical  tc 


pe  =  Pr|<  1  5L  > 


— w 


(3.39) 


where  ^  =  W£  and  m^  =  Wm,  since  as  mentioned  above,  the  nonsingular 

transformation  W  is  known  to  exist  when  R  (t,s)  is  strictly  positive 

n 
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definite,  and  in  fact  is  given  by  Eq.  (3.27).  Now,  it  is  seen  from  Eqs. 
(3,36),  (3.35),  and  (3.37)  that  p^,  )  f^r  both  generalized  "t"  noise 

and  modified  generalized  "t"  noise  is  a  function  only  of  the  length  of 


the  noise  vector;  i.e., 


(4)  ■  • 


(2.40) 


Thus,  because  of  the  spherical  symmetry  of  this  probability  distribution, 
it  follows,  as  shown  geometrically  in  Fig.  21,  that  the  calculation  of 
reduces  to  the  one-dimensional  calculation  independent  of  direction: 


P  =  Pr 
e 


(y  >  ill  • 

|w  —  2  — w  [ 


(3.41) 


That  is,  the  spherical  symmetry  of  p^  (jr^)  means  that  the  probability 
that  the  component  of  the  noise  vector  in  the  direction  of  m^  is  greater 
than  or  equal  to  ( l/2 ) j|m  ||  is  independent  of  the  direction  of  m  . 


CONTOURS  OF 
P,  (y J  s  CONST 


t  ll  mwll  A 


y*.  \ 


<2,=  o 


FIG.  21.  GEOMETRICAL  PICTURE  OF  SPHERICAL  SYMMETRY  OF  THE 
PROBABILITY  DISTRIBUTION  OF  THE  CASE  I  ADDITIVE-NOISE 
VECTOR . 
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5 .  Genera llzed  "t"  Noise 


Applying  the  above  results  to  the  case  of  generalized  "t"  noise, 
it  is  found  by  setting  N  =  1  in  p  (j^,)  given  by  Eq.  (3.26)  that 


- - -J75  .  -»<*<“.  (3.42) 

:  *  M 


where  Q  =  m  +  1  >  1  and  =  m  l/c.  Thus,  using  this  probability 
density  to  evaluate  Eq.  (3.41 ),  there  results: 


P  = 
e 


(3.43) 


Now,  it  is  necessary  to  specify  0  in  order  to  evaluate  this  integral  in 
closed  form.  Initially,  consider  the  important  case  0=3  (see  Section 
IICl)  which  gives 


(3.44) 


In  order  to  evaluate  the  significance  of  this  result,  it  is  convenient  to 
consider  two  subcases  of  interest,  corresponding  to  high  and  low  average 
SNR: 

a.  high  SNR;  i.e.,  jjmjl  »  7y/ 

In  this  case  there  is  found  the  asymptotic  result 


— * . — -  00 


(3.45) 


* 


SEL-66-052 


94 


b.  Low  SNR;  i.e.,  ||mj|  «  y w 

la  this  case  there  is  found  the  asymptotic  result 


Discussion  of  these  results  for  the  generalized  "t"  model  will  be  post¬ 
poned  until  the  associated  results  for  the  modified  generalized  "t" 
model  have  been  determined.  Plots  of  P^  vs  li^ll/^  are  presented  in 
Fig.  22  for  both  the  generalized  "t"  and  modified  generalized  "t" 
models.  These  error  curves  are  presented  to  show  the  effects  of  various 
types  of  noise  on  the  optimal  performance  achievable  using  "short  duration" 
signals.  For  purposes  of  comparison,  the  gaussian  noise  curve  is  also 
presented  for  a  noise  whose  average  power  is  equal  to  that  typically  mea¬ 
sured  in  atmospheric  noise  situations  where  the  generalized  "t"  model 

$ 

_ith  parameters  0  =  3  and  y  is  appropriate. 

'  w 


erfc(x)  =  I  - exp 

denotes  the  area  under  the  tail  of  the  normal  probability  density  with 
zero  mean  and  unity  variance. 

This  normalization  is  accomplished  by  using  the  fact  that  for  the  gen¬ 
eralized  "t"  model 

I  »  §  7W  .  (2.28) 

0=3  *  W 

together  with  the  empirical  result  at  temperate  latitudes  that,  typically 

at  vlf  and  If,  V  /V  «  3  [Refs.  9,  10], 
rms'  avg  ' 
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PROBABILITY  OF  ERROR, 


WWftrt'B*.! 


FIG,  22.  “SHORT-DURATION"  SIGNAL  (CASE  i)  PERFORMANCE:  COMPARISON  OF 
ERROR  CURVES  ACHIEVABLE  BY  OPTIMAL  MATCHED  FILTER  RECEIVER  IN  PRESENCE 
OF  VARIOUS  ADDITIVE  NOISES. 

(Note:  Vertical  scale  is  log| log  Pj  . ] 


LEGEND: 

— — — —  Modified  generalized  "t  "  noise  model 

0  =-  2 

e  =  io3 

Generalized  “t"  noise  model 
0=3 


Gaussian  noise  model 

Same  average  power  as  typically 
observed  when  the  generalized 
"t"  model  with  parameters  0=3 
and  ^  is  appropriate. 
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6.  Modified  Generalized  ,!t"  Noise 


Turning  attention  to  modified  generalized  !,t"  noise  and  taking 
0  =  2,  since  this  is  the  simplest  case  of  interest,  it  is  found  by  set¬ 
ting  N  =  1  in  P^(^)  given  by  Eq.  (3.35)  that 


p  (x)  |  =  — 

yw  0=2  (2 

m 


K _ 2  f  1/2  2.' 

\l/2  2  2  esp  “  2  +  yj  * 

it)  '  x  +  y  L  20  J 


where  y  =  l/cr  and 


’  J? 


1  7w 

H  da 

a 


-00  <  X  <  00  9 

(3.47) 


(3.48) 


Thus,  evaluating  Eq.  (3.4l)  in  this  case, 

00  , 

f  /Ov1/2 


1  (!)  exp[-  ~2 (x2 + *] dx  •  (3-49) 
x  +  y  L  20  J 


x!in>  ‘i 

2  — w 


Now,  this  integral  cannot  be  evaluated  in  closed  form  as  it  stands;  but 
it  can  be  approximated  in  two  subcases  of  interest,  corresponding  to  high 
and  low  SNR: 

a.  High  SNR;  i.e., 

In  this  case  is  found  the  approximate  result 
f  1/2  K  1  I  7w  2 

I  /!>  \  '  —  — —  pyn  I—  .  y  dy 


J  (!) 


7  2 

'  w  y 


1  I'm  '! 

2  — w 


r  „  .-2" 

i 

l\ 

2 

!  ''V 

-  2 ^  orfpf 

,  exp 

!m  » 

— w 

!  802 

T  erfc\ 

t  20  / 

j 

(3.50) 
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In  order  to  complete  this  evaluation,  K  must  be  determined  from  Eq. 
(3.48),  which  gives 


2(2«)l/2erfc^j 


(3.51) 


However,  this  expression  can  be  simplified  by  using  the  fact  that  (3  »  y 

is  satisfied  in  the  application  of  the  modified  generalized  "t"  model 

t 

to  atmospheric  noise.  Now 


lim  K  =  - 

ZiL  _  0  (2n)] 

3 


-&) 


Therefore,  substituting  this  result  into  Eq.  (3.50), 


(3.52) 


7  o  m 

„  i  'w  2  -v, 

P  as  -  - fr  exp - — 

ee=2  n  '^w'l  8^2 


■(¥))  ■ 


(3.53) 


This  follows  from  the  results  of  Section  IIC1.  From  their  definitions 


(3  fr1  3  C'(3  So  ’ 

since  m  =  1  here.  However,  in  the  application  of  the  modified  general¬ 
ized  "t"  model  in  Fig,  10  it  is  seen  that  (3c r=  5  X  10^  is  required  to 
fit  the  measured  data.  Thus,  for  this  case,  which  is  typical  of  results 
at  vlf  and  If, 

7W  -3 

T  =  2  X  10  «  1  . 
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In  order  more  easily  to  compare  the  form  of  this  high  SNR  result  with 
that  given  by  Eq .  (3.45)  for  the  generalized  "t"  model,  it  is  illumi¬ 
nating  to  consider  two  further  subcases  as  follows: 

(0  i. 

In  this  case  there  results 


2  7W 

n  II™  j 


(3.54) 


»  1: 


In  this  case  there  results 


r2 

8  P  7w 

P  |  = - r  exp  - 

e9=2  n  ||m  f 

— w  L 


b.  Low  SNR;  i.e.,  ilmjj  «  y^ 


1/2  40  y 


1  erfcf— ). 

2  1C\  2p  } 


(3.55) 


Rewriting  P  |  given  by  Eq.  (3.49)  in  the  form 
eQ=2 


1  'V 


f  (2  \ly^2  K  1  7w  /  2  . 

J  U  T~2 — ; exp  -  “2  (y  + 1}  dy  - 

q  7w  y  +  1  L  23 


(3. .5) 


it  is  seen,  using  the  fact  7W/P  «  1  in  conjunction  with  the  low  SNR 

condition  '|m  i '/y  «  1,  that 

— w  w 


.  m 
1  — w 


'  w 

/ 


(3.5/) 
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Thus,  for  the  low  SNR  case: 


l 

n 


t 


-1 

an 


lim  erfc  (( 2 n ) 


1/2 


•'v 


(3.58) 


7.  Discussion  of  Probabillty-of-Error  Results 

Before  proceeding  to  the  discussion  of  the  Pe  results  calculated 
above,  it  is  recalled  that  these  Case  I  results  apply  to  the  situation 
where  the  additive  atmospheric  noise  is  given  by  y(t)  =  a  n(t).  0  <  t  <  T; 

i.e.,  they  apply  when  the  transmitted  signal  duration  T  is  short  enough 
that  a(t)  =  a  can  be  assumed  constant  for  the  duration  of  the  signal. 

With  this  in  mind,  the  results  for  Case  I  lead  to  the  following  con¬ 

clusions  : 


a.  Generalized  "t"  Noise 

The  description  of  the  lightning  discharge  presented  in  Section 
IIA  indicates  that  the  Pg  results  for  Case  I  apply  at  vlf  when  the  dura¬ 
tion  of  the  transmitted  signal  is  less  than  about  0.1  msec.  In  addition, 
inspection  of  the  plotted  results  in  Chapter  II  indicates  that  these  results 
apply  oractically  to  any  situation  at  vlf  or  If  where  the  duration  cf  the 
transmitted  signal  is  less  than  or  equal  to  the  inverse  of  the  receiver 
bandwidth  for  bandwidths  greater  than  a  few  hundred  cycles  per  second. 


This  follows  as  an  engineering  approximation  from  the  fact  that  the  rise 
time  of  a  signal  passed  through  a  bandpass  filter  of  bandwidth  2B  is 
about  l/2B  seconds;  i.e.,  2TB  <1  in  these  situations.  This  is  consis¬ 
tent  with  the  short-duration  signal  idea,  since  receiver  bandwidth  is  re¬ 
stricted  in  practice  at  both  vlf  and  If  because  of  the  lack  of  spectrum 
availability  (see  Section  IIA). 
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To  the  author's  knowledge,  there  are  no  vlf  communication  systems  in  oper¬ 
ation  with  either  of  these  characteristics.  There  is  good  reason  for  this 
(in  addition  to  the  difficulty  ol'  obtaining  transmitting  antenna  bandwidths  v 

of  this  size  at  vlf),  as  indicated  by  the  P  results.  In  particular, 

e 

the  high  SNP.  result  given  by  Eq .  (3.45)  indicates  that,  even  when  using 

the  optimal  receiver,  the  ..siprcvement  in  P  with  increased  signal  energy 

e 

at  low  error  rates  is  v  >ry  small  compared  to  the  improvement  achieved  in 
the  case  of  additive  gauss^an  noise  (see  Fig,  22).  This  indicates  that 
short-duration  signals  are  not  an  optimal  choice  for  use  in  the  presence 
of  vlf  or  If  atmospheric  noise  if  it  is  desired  to  achieve  economically  a 
low  probability  of  error;  this  is  physically  reasonable  because  of  the  dis¬ 
tinctly  impulsive  nature  of  atmc  pheric  noise. 

b.  Modified  Generalized  "t"  Noise 

Inspection  of  the  plotted  results  in  Chapter  II  indicates  that 

the  P^  results  for  Case  I-  apply  practically  in  two  situations:  those 

situations  at  mf  or  hf  in  which  the  duration  of  the  transmitted  signal  is 

less  than  or  equal  to  the  inverse  of  the  receiver  bandwidth,  and  those 

situations  at  vlf  or  If  in  which  the  duration  of  the  transmitted  signal  is 

less  than  or  equal  to  the  inverse  of  the  receiver  bandwidth  for  bandwidths 

less  than  a  few  hundred  cycles  per  second.  Once  again,  the  most  interesting  * 

results  are  those  pertaining  to  the  high  SNR  case,  given  by  Eqs.  (3.54)  and 

(3.55).  These  results  are  demonstrated  in  Fig.  22,  where  it  is  seen  that 

the  improvement  in  P  with  increasing  ||m  ||  is  greatly  enhanced  when  the 

e  —w 

inequality  »  1  is  satisfied. 

In  regard  to  th"  practical  applicability  of  these  results  for 
the  modified  generalized  "t"  model,  the  plotted  data  in  Chapter  II  Lni  J  - 
cate  that,  for  a  fixed  receiver  bandwidth,  the  value  of  P/Vw  decreas.es 
as  the  opc  ‘ating  frequency  is  increased.  Thus  the  use  of  short-duration 
signals  is  more  attractive  at  hf  than  at  lower  operating  frequ  ncies. 

► 

C.  CASE  II;  LONG-DURATION  SIGNALS;  THE  GENERAL  CASE 

The  probability-of-error  results  obtained  in  Case  I  indicate  that  short- 
duration  signals  are  not  an  optimal  choice  in  the  presence  of  additive  at-  * 

mospheric  noise.  This  follows  from  the  fact  that,  particularly  at  vlf  and 
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If  where  the  generalized  "t"  model  is  appropriate,  the  Pg  resulting 
from  the  use  of  short-duration  signals  behaves  in  a  much  less  desirable 
fashion  than  would  the  Pe  if  the  additive  noise  were  gaussian.  Thus 
it  Is  desirable  to  investigate  the  use  of  long-duration  signals,  where 
"long  duration"  means  that  the  Case  I  assumption  that  a(t)  =  a  is  con¬ 
stant  for  the  duration  of  the  signal  mus'  now  be  relaxed.  This  consider¬ 
ation  of  long-duration  signals  is  well  motivated  physically,  since  we 
would  like  to  be  able  to  disregard  the  received  signal  when  it  i  largely 
the  result  of  a  pulse  of  noise,  and  base  our  decision  on  the  relatively 
noise  free  signal  received  between  t^ese  noise  pulses.  This,  in  f.at,  is 
the  procedure  typically  employed  in  practice,  where  long-duration  signals 
ar-5  used  in  conjunction  with  nonlinear  receiving  techniaues: 

In  this  case,  the  decision  problem  at  the  receiver  becomes  that  of 
choosing  between  the  two  £  priori  equally  probable  hypotheses 


h^:  x(t)  =  m(t)  +  a(t)  n(t) 
h(2):  x(t)  =  a(t )  n(t) 


0  <  t  <  T  . 


(3.59) 


The  desired  vector  formulation  in  this  case  is  developed  as  follows:  De¬ 
fining  q^t),  i=l,  ...,  N,  as  ^hown  in  Fig.  16,  and  assuming  the  additive 
noise  y(t)  =  a(t)  n(t)  to  have  a  continuous  covariance  function,  it  is 
possible  to  write  [Ref.  41  ] 


y ( t )  =  l.i.m. 

N  -»  co 


N 

£  y(‘i)  (c‘i)l/2  *,(*>■ 

i=l 


0  <  t  <  T  .  (3. 60) 


Since  it  is  also  possible  to  write 


N 


i(t)  =  l.i.m.  ^  m(ti)  (At^)1/2  qi'.t), 


i=l 


0  <  t  <  T  (3. 61) 
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for  all  sqi'-'re-integrable  signals  m(t),  it  is  seen  that  the  vector  for¬ 
mulation  of  the  decision  problem  can  be  written  as  follows  with  respect 
to  the  orthonormal  basis  functions  q^(t),  i  =  1,  ...»  N:  (As  in  Case  I, 
N  will  be  taken  to  infinity  at  the  conclusion  of  the  problem.) 


.  (l) 

h  *  =  H  +  1  i 


where,  e.g.,  =  ,yi >  •  •  * »  yN,  with  yi  =  y(ti)  (At^1 

1,  Calculation  of  the  Likelihood  Ratio 


(3.62) 


w(*lh^)  =  "  ®)  , 


(3.63) 


where  the  probability  density  function  p^(^)  the  additive-noise  vector 

is  specified  as  follows:  Recalling  that  a(t)  is  now  allowed  to  vary 
"slowly”  on  tUe  interval  [0,t]  one  can  write 

yi  =  y^i^  =  a(ti^  n(ti)  (^ti)1/2 


=  a4  n4. 


i  =  1,  ....  N  , 


(3.64) 


where  n.^  =  n(t^)  (At^^)1^2  and  =  a{t^).  (it  is  noted  that  this  is 
consistent  with  the  interpretation  in  Chapter  II  of  a(t)  as  a  dimension¬ 
less  weighting  factor.)  Since,  as  in  Case  T,  n(t)  is  a  zero-mean  gaus- 

sian  process  with  covariance  function  R  (t,s),  the  distribution  p  (n) 

n  n  — 

of  the  vector  n  =  n.,  ....  n  is  given  by 
-  i_i _ Hi 

pn(n)  = - JJ75 - TK  exp[-  ^  n  r”1  n]  ,  (3.65) 

2  (2„)N/2  |R  i  1/2  L  2-t  n  -j 
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where 


» 


R 

n 


t  )  (At  At  ) 

J  A  J 


1/2 


(3.66) 


In  order  to  complete  the  specification  of  the  problem,  the  distribution 
of  the  rector  a  of  coefficients  must  be  specified.  However,  recalling 
that  a(t)  =  l/b(t),  the  distribution  p^(b)  of  the  vector  b  = 
b  ,  ....  b,  can  be  equivalently  specified.  Now,  the  first-order  statis- 
tics  of  b(t)  are  given  by  the  two-sided  chi  distribution  X2(m,cj).  How¬ 
ever,  in  the  special  case  9=2  (m  =  l)  this  reduces  to  the  normal 
distribution  N(o,a  ),  so  that  b(t)  is  a  gaussian  process  in  this  spe¬ 
cial  case.  Thus,  initially  considering  the  case  0=2, 


(3.67) 


where 


R 


b 


(3.68) 


/  \  2 

in  which  it  is  noted  that  R^t^,  )  =  cr  in  the  notation  previously 

used.  Now 


Si)  =  j  P^zlb)  Pb(b)  db  . 


where  setting  =  ni/bi*  i  =  1,  ...,  N,  gives 


p ^i\k)  = 


N 

n  b. 

i=i  1 


>n(ni  =  bi  V  1  =  1 . N) 


(3.69) 


(  2  jt ) 


N/2 


1/2 


b, 


(3.70) 
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in  which  b,  y.  =  b  y  ,  ....  b  y  .  Thus,  substituting  into  Eq.  (3.69), 

11  11  ii  ri , 

» - 1  l _ I 

0°  00 

P!(i)  =  (2„>»  |H  .  ^  I  ** .  /dxA|bl1) 


exp[-  5  (fiji  C  bi  yi)  +  fit  "b1  fi)] 


(3.71) 


In  order  to  evaluate  this  N-fold  integral,  it  seems  reasonable  to  diag¬ 
onalize  simultaneously  the  two  quadratic  forms  'n  the  exponent.  This  can 
be  done  (Ref.  43],  but  requires  the  use  of  linea.*  transformations  depen¬ 
dent  in  an  unknown  way  on  the  vector  As  a  result,  it  does  not  appear 

feasible  to  pursue  this  general  formulation  without  specification  of  one 
or  the  other  (or  both)  of  Rn(t,s)  and  Rjj(t,s).  Actually,  it  is  easily 
shown  that  both  Rn(t,s)  and  Rb(t,s)  must  be  specified  if  p^(j)  is 
to  be  evaluated  as  a  known  function  of  hence  we  will  proceed  to  con¬ 

sider  special  cases  of  interest  as  follows: 

The  special  case  that  first  comes  to  mind  is  a  "white"  noise  case: 


Rn(t,s)  =  NQ/2  6(t-s)  , 


(3.72a) 


Rb(t,s)  -  Bq/2  S(t-s)  , 


(3.72b) 


but  this  case  will  not  be  further  considered  here  since  it  is  not  appli¬ 
cable  to  the  communication  channels  of  interest.  This  follows  from  the 
results  of  Chapter  II,  where  it  is  shown  that  the  generalized  "t"  model 
is  in  practice  an  appropriate  model  for  received  atmospheric  noise  when 
n(t)  is  a  narrowband  random  process  and  a(t)  =  l/b(t)  is  a  slowly- 
varying  random  process.  While  these  results  discourage  the  "whiteness" 
assumption,  they  do  suggest  alternate  ways  of  specifying  the  higher-order 
statistics  of  a(t)  and  n(t).  One  such  specification  which  is  of 


Note  that  b  y  is  the  transpose  of  the  column  vector  b 

t  *  _  1 
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particular  interest,  since  it  provides  a  formulation  retaining  much  of 
the  simplicity  of  the  white  noise  case,  is  considered  below  as  Case  Ila. 

D.  CASE  11a:  COMPLEX  ENVELOPE  REPRESENTATION 

Recalling  from  the  discussion  in  Section  1IB  that  in  practice  atten¬ 
tion  can  be  restricted  to  the  situation  in  which  the  receiver  bandwidth 
is  substantially  less  than  the  band  center  frequency,  the  Case  II  hypothe¬ 
ses  can  be  further  described  by  noting  that  x(t),  m(t),  and  y ( t )  =  n(t) 
n(t)  can  be  considered  to  be  narrowband  functions  of  time.  Thus  it  fol¬ 
lows  that  the  Case  II  decision  problem  can  be  written  directly  m  terms 
of  the  slowly  varying  complex  envelopes  [Ref.  40 ]  of  x(t),  m(t),  and 
y(t).  In  terms  of  these  complex  envelopes,  the  decision  problem  becomes 
that  of  choosing  between  the  hypotheses 


hU;:  x(t)  =  p(t)  +  n(t) 
h(2):  x(t)  =  T,(t) 


0  <  t  <  T  , 


(3.73) 


where,  e.g. , 


x(t)  =  *c(t)  cos  2«fpt  -  xg(t)  sin  2nf^t 


f  i2«f0t\ 

!|*(0  O  0  ) 


(3.74) 


and  x(t)  =  x£(t)  +  ixg(t)  is  the  slowly  varying  complex  envelope  of 
x(t).  Proceeding  as  in  the  Case  II  development,  one  can  write 


■(  t )  =  l.i.m.  ^  x  ( t  )  (At  J1/2  q  (t), 

N  -»  oo  ^  1  1  1 

i=l 


0  <  t  <  T  .  (3.75) 
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Furthermore,  it  is  proposed  that  a  practically  useful  model  results  from 
assuming  that  the  slowly  varying  complex  envelope  is  given  by 

N 

x(t)  «  ^  x(ti)  (Ati)i//2  qt(t)  ,  0  <  t  <  T  ,  (3.76) 

i=l 


where  the  value  of  N  required  will  be  discussed  further  below  and  is 
dependent  on  the  bandwidth  in  which  the  signal  of  duration  T  is  observed. 
Proceeding  again  as  in  Case  II,  the  vector  formulation  of  this  problem 
with  respect  to  the  orthonormal  basis  functions  q  (t),  i  =  1,  ....  N, 
can  be  written  as 


h^ 1  ^ :  *■  =  t±.  +  H 


h<2>: 


(3.77) 


1/2 


where,  e.g.  ,  ^  ....  nN|  with  =  q^)  (At^ 

1 .  Calculation  of  the  Likelihood  Ratio  for  Case  Ila 


Now 


w(xjh^)  =  p^(x  -  £)  , 


(3.78) 


where  the  probability  density  function  P^(d)  of  the  complex  additive 
noise  vector  is  specified  as  follows:  Recalling  that  a(t)  is  itself  a 
slowly  varying  random  process,  one  can  write 


=  n<t1)  (Att)l/2  =  „(t4)  v(tt)  (At4)l/2  =  v4  (3.79) 

where  a4  =  aft  )  and  v4  =  vft^)  (At^)1  ,  in  which  v ( t )  is  the  com- 

plex  envelope  of  n(t).  Since  v(t)  is  the  complex  envelope  of  a  gaus- 

sian  process,  it  follows  that  it  is  a  complex  gaussian  process.  Thus  the 

probability  density  function  p  (.v)  of  the  vector  v  =  ,  ...,  v  can 

—  ~~  *  *  1  £Li 

be  written  (assuming  the  power  spectral  density  of  the  narrowband  process 
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n(t)  to  be  symmetric  about  the  band  center  frequency)  [Refs.  40,  4l]  in 


the  form: 


“yW  -  7^7  -l  -]  1 


(3. 80) 


where 


0  =  E[v  v  1  , 
v  - 1 


(3. 81) 


E[v  v*]  =  E[v(  t  )  (At  At  )l/ 

A  J  •  J  **■ 


(3.82) 


Motivated  by  the  difficulties  encountered  in  the  Case  II  calculations,  it 
is  proposed  that  an  interesting  case  for  initial  consideration  is  given 


E[v(ti)  V*(t.)J  =  Nq  61J( 


i,j  =  1,  ....  N . 


(3.83) 


Substituting  this  specification  into  Eq.  (3.8O), 


Pv(*}  =  ~ 


it  (N0At)' 


1  * 

’['  NoAt  - 


(3.84) 


where  t’  e  fact  that  At.  =  At  for  all  i  =  1,  ...,  N  has  been  used. 

It  now  remains  to  specify  the  probability  density  function  of  the 

/ector  a  =  a  ,  . . . ,  a„  of  coefficients.  Following  Case  II  above,  the 

1  t_l - b 

equivalent  specification  of  p^b)  will  again  be  made  where  9  =  2  is 
assumed  so  that  b(t)  is  a  gaussian  process.  Thus  pb(b)  is  given  by 
Eq.  (3.67).  However,  again  motivated  by  the  Case  II  results,  it  is  pro¬ 
posed  that  a  case  of  interest  follows  from  the  assumption 


E[b(t, )  b(t  .)]  =  —  6  ,  i, j  =  } ,  . . . ,  N  . 


(3.85) 
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Although  the  representation  for  the  received  signal  given  by  Eq.  (3.76) 

is  particularly  convenient  in  the  detection  problem  being  considered  here, 

it  remnins  to  verify  the  validity  of  this  representation  in  the  light  of 

the  assumption  on  v  and  R  maue  in  Eqs .  (3..  )  and  (3.85),  respec- 

v  0 

tively.  In  order  to  do  this,  it  will  be  recalled  that  the  received  signal 
x(t)  has  one  of  the  forms 

h*'1^:  x(t)  =  m(t)  +  y(t )  j 

l  0  <  t  <  T  ,  (3.59) 

h^2^:  x ( t )  =  y ( t )  j 

where  y(t)  =  (l/b(t)j  n(t).  Then  the  following  observations  can  be  made: 

1.  As  f ar  as  the  observed  noise  y(t)  is  concerned,  the  representation 

proposed  by  Eq.  (3.76)  with  <t>v  and  given  by  Eqs.  (3.83)  and 

(3.85)  is  identical  in  form  to  the  representation  considered  in  the 
case  of  independent  samples  in  Section  IIC3.  It  follows,  from  the 
results  shown  there,  that  this  representation  is  a  reasonable  first 
approximation  to  the  complex  envelope  ’l(t)  when  N  ar  2TB,  where 
2B  is  the  rf  bandwidth  in  which  the  received  signal  of  duration  T 
is  observed. 

2.  As  far  as  the  known  signal  m(t)  is  concerned,  it  will  be  found 
that  in  practice  this  signal  of  duration  T  has  a  bandwidth  signifi¬ 
cantly  le3s  than  the  receiver  bandwidth  2B.  Thus  the  representation 
proposed  by  Eq.  (3.76)  with  N  =»  2TB  gives  a  good  approximation  to 
the  complex  envelope  p(t). 

Therefore,  combining  observations  (l)  and  (2),  it  is  concluded  that  the 
representation  given  by  Eq.  (3.76)  with  N  =  2TB  is  a  reasonable  repre¬ 
sentation  for  the  received  signals  of  interest  in  Case  Ila  when  4^  and 

R.  are  given  by  Eqs.  (3.83)  and  (3.85).  Now,  substituting  Eq.  (3.85) 
b 

into  Eq,  (3.67), 

Pb(^}  =  ■n/21;n72  exp 


=  bt  rti,  i  =  1,  ....  N)  ,  (3.87) 


Now 


— •  b  b 

“o”1  " 


(3.86) 
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where  use  has  been  made  of  the  fact  that  in  terms  of  real  variables 


lV(vi)  *  V  n  (”ic-  ;\s>  ■  V  Kc’  pn  (nl.>  •  <3'8‘' 

i  ic,  is  ic  is 


so  that  the  Jacobian  of  the  transformation  is 


j(v  -*  Ti  )  =  b 
v  i  i '  i 


(3.89) 


Therefore, 


Pa<l|£> 


1 


N 

Jt 


(«0At)' 


N 


_ 1__ 

"oAt 


1 


I  b. 


(3.90) 


so  that 


3N  / 2 


(»oat)N 


00 


—00 


(2  nf 


N 

II 

i=l 


1 


NoAt 


3/2  • 


(3.91) 


t 
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(  2 ) 

Now,  using  Eq .  (3.78)  and  the  similar  expression  for  h'  ,  the  likeli¬ 
hood  ratio  can  be  written  directly  as  follows: 


w(x)h^)  N 

l(x)  = - nry-  =  n 

w(x| h^  )  i=l 


L  ,>  tNoAtl 

3/2 

J  1  V 

Xi  -  + 


„  N  At 
2  0 


T572 


Thus  the  Bayes  rule  in  this  case  is  given  by: 


Choose  h^1^  iff 


N 

n 

i=l 


.  2  N0At 


>  [l*,  - 


2  NoAt 

+ 


B 


0  J 


(3.92) 


(3.93) 


However,  this  is  equivalent  to  the  rule 
Choose  h^1)  iff 


I* 

i=l 


2  N0 


l^i>!  ♦sj 


N 

*  2  *■ 


i=l 


(3.94) 


Therefore,  noting  that  it  has  been  assumed 

x(t)  =  x(tt)  V  t  e  AtA  , 
m(0  =  n(tt)  V  t  e  , 


(3.95a) 

(3.95b) 


it  follows  that,  in  terms  of  the  time  functions  available  to  the  receiver, 
the  Bayes  rule  is  given  by 

Choose  h^^  iff 


T 

f 

N  ' 

T 

r 

O  " 

in 

!x(t)|2 

dt  >  1  in 

MO  -„<«)|2  +  £ 

J 

0 

L  oj 

J 

0 

0. 

(3.96) 


f 
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Noting  that  jx(t)j  is  the  envelope  of  x(t),  and  |x(t)  -  ,-(t)|  is 
the  envelope  of  x(t)  -  m ( l ) ,  it  is  seen  that  one  receiver  structure 
which  implements  this  rule  is  that  shown  in  Fig.  23. 


m(t)-»4  2 
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FIG.  23.  BLOCK  DIAGRAM  OF  OPTIMAL  LOGARITHMIC-CORRELATOR  RECEIVER. 


2.  Discussion  of  Optimal  Receiver  Rule  (Case  Ha, 

In  this  case,  which  is  proposed  as  perhaps  the  implest  formulation 
of  the  long-duration  signal  problem  of  practical  interest,  several  comments 
are  warranted.  First,  it  is  interesting  to  compare  the  receiver  rule  given 
by  Eq.  (3.96),  and  its  implementati  -r  shown  in  Fig.  23,  with  intuitive  no¬ 
tions  concerning  receiver  operation  in  the  presence  of  additive  atmospheric 
noise.  It  is  seen  from  Eq.  (3.96)  that  this  nonlinear  receiver,  which  will 
be  called  here  the  "logarithmic-correlator"  receiver,  does  indeed  agree 
with  intuition,  since  it  suppresses  consideration  of  the  received  signal 
at  ihose  times  when  the  quantities  |x(t)|  and  |x(t)  -  p(t)|  are  both 
large,  i.e.,  those  times  when  the  received  signal  is  predominantly  due  to 
a  pulse  of  noise,  ard  bases  its  decision  strongly  on  the  received  signal 
at  those  times  wher  one  of  the  quantities  jx(t)|  or  |x(t)  -  p(t)j  is 
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small.  Thus,  heuristically,  the  log-correlator  receiver  ignores  the  re¬ 
ceived  signal  when  it  is  largely  the  result  of  a  pulse  of  noise  and  buses 
its  decision  on  the  relatively  quiet  periods  between  the  occurrence  of 
these  noise  pulses.  Furthermore,  it  is  noted  that  the  receiver  becomes 


increasingly  nonlinear  as  the  average  SNll  increases,  which  (again  in  agree¬ 
ment  with  intuition)  means  that  the  receiver  has  increasing  confidence  in 

2  2 

the  received  signal  at  those  times  when  either  jx(t)j  or  |x(t)  “  p(0| 
is  small  as  the  average  SNR  increases. 

Although  the  calculation  of  the  receiver  rule  Eq.  (3.96)  was 
carried  out  for  the  special  case  9=2,  it  is  seen  that  the  Case  Ila 
calculations  could  have  equally  easily  been  carried  out  without  specifi¬ 
cation  of  0,  by  making  the  assumption  that  the  received  noise  is  statis¬ 
tically  independent  at  the  instants  t^  and  t ^  for  all  i,j  =  1,  ...,  N. 
When  the  calculation  of  the  likelihood  ratio  is  carried  out  with  this  as¬ 


sumption,  it  is  seen  that  the  decision  rule  becomes: 
Choose  h^1 ^  iff 


where  it  is  recalled  that  9  =  m  +  1  >  1.  Thus  it  is  seen  that  the  form 
of  the  optimal  decision  ruie  is  not  changed  as  m  is  varied,  the  only 
effect  on  the  realization  shown  in  Fig.  23  being  a  change  of  the  bias  on 
the  logarithmic  amplifier.  It  is  noted  that  such  a  modification  is  rea¬ 
sonable,  since  the  average  noise  power  is  proportional  to  "'(Nq/Hq),  with 
m  =  1  leading  to  Eq,  (3. 96). 

It  is  clear  that  the  realization  of  the  optimal  receiver  rule  given 
in  Fig.  23  is  net  the  only  realization  possible.  It  is  desirable  to  find 
realizations  that  are  simpler  in  the  sense  that  the  number  of  zero-memor> 
nonlinear  operations  required  is  reduced.  Specialization  to  the  high  SNR 
case  does  not  appear  to  lead  to  any  such  simplifications,  but  specializa¬ 
tion  to  the  low  SNR  case  does  as  follows: 
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Low  SNR  case,  i.e., 


fell2  «  I2  . 


whore 


2  A  N0 
5  =  m  —  At  . 

0 


In  this  case  it  is  easily  shown  that  the  optimal  decision  rule  is 


given  approximately  by 
Choose  h^  iff 


i  r _ *(o  /it) 

1°  ix(t)i2  + 


dt\  > 


l\M\' 


J  r  . 

o  jx(t 


dt  .  (3.98) 


Thus  it  is  seen  that  the  low  SNR  assumption  gives  an  interesting 
rule  which  is  in  fact  the  correlator  receiver  modified  by  the  weighting 
of  each  instantaneous  contribution  by  a  factor  dependent  on  the  value  of 
the  envelope  of  the  received  signal  at  that  instant.  However,  closer 
examination  of  the  derivation  of  Eq.  (3.98)  shows  that  the  low  SNR  assump¬ 
tion  can  actually  be  relaxed,  since  the  condition  required  to  produce 

2  2 

Eq.  (3.98)  is  precisely  that  |  pj  «  £  V  i  =  1,  . ..,  N.  Now,  this  con¬ 
dition  does  not  necessarily  imply  the  low  SNR  case,  and  in  fact  may  be 

2  2 

compatible  with  the  high  SNR  case,  i.e.  vith  j|yj|  »  |  ;  so  that  the  re¬ 
ceiver  rule  given  by  Eq.  (3.98)  may  be  of  practical  interest.  This  possi¬ 
bility  will  be  examined  further  in  Section  E  and  is  of  importance  because 
Eq.  (3.98)  has  a  realization  similar  to  the  matched-filter  realization  of 
the  correlator-receiver  rule. 

Finally,  it  is  interesting  to  apply  the  Case  Ila  assumptions  to 
the  situation  where  the  additive  noise  is  given  by  the  modified  generalized 
"t"  model,  rather  than  by  the  generalized  "tM  model  considered  above. 

The  calculation  of  the  Bayes  rule  appropriate  to  this  situation  can  be 
carried  out  in  the  same  manner  as  that  demonstrated  above  for  additive 
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generalized  "t"  noise.  The  decision  rule  for  the  special  case  9  =  3 
is  found  to  be: 

Choose  h^^  iff 


r  „ 


* 

l*(0|2  +2^ 
L  oJ 

2 

2 

|x(t)|2  +  N0^2  H 

■d 

*  2 
Nop 

Re[x(t)  p*(t) ]  \  dt 


>  j  hn 


I  jjx(t 


i2  1 


)  “  n(0| ' 


+  2 


I 

- 


x(t)  -  p(t)|  + 


•■p-% 


V 


|n(t)|: 


dt 


(3.99) 


This  result  is  interesting,  since  it  allows  a  quantitative  check  as  to 
whether  or  not  the  log-correlator  receiver  is  appropriate  in  a  given  situa¬ 
tion.  It  clear  that  the  rule  given  by  Eq.  (3.99)  reduces  to  the  log- 

2 

correlator  receiver  rule  as  p  -»  and  inspection  of  the  experimental 
dnta  presented  in  Section  IIC1  shows  that,  in  fact,  the  log-correlator 
receiver  is  the  receiver  of  practical  interest  at  vlf  and  If  frequencies 
of  operation.  On  the  other  hand,  no  such  statement  can  be  made  when  the 
operating  frequency  is  increased  to  hf.  In  this  case,  Eq.  (3.99)  must  be 
examined  for  each  particular  situation.  This  follows  from  the  fact  noted 
in  Section  IIC1  that  the  nature  of  the  received  noise  at  hf  is  strongly  a 
function  of  local  thunderstorm  activity. 


3.  Calculation  of  the  Probability  of  Error  (Case  Ila) 


The  receiver  that  minimizes  the  probability  of  error  in  Case  Ila 
is  the  one  that  implements  the  rule  given  by  Eq.  (3.93).  Thus  the  prob¬ 
ability  of  error  resulting  from  the  use  of  the  optimal  receiver  is  given 
by 


% 
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f 

N 

n 

2  No  1 

hi!  +  m  At 

* 

>  1 

N 

i=i 

\\  ~  Mil  +  m  b“  At 
0 

4 

(3.100) 


However,  (3.100)  is  equivalent  to  the  more  manageable  formulation 


P  =  Pr( 
e 


N 


2* 


i=l 


,2  N0 

\\  +  m  At 


2  N0 

\\  -  Hj  +  m  —  At 


>  0)  =  Pr< 


N 

l n  z^  >  0 
i=l 


(3.101) 


where 


Zi  = 


2  N0 

Mi  -  ^  |  +  m  — At 

0 


(3.102) 


Since  the  z are  statistically  independent,,  it  follows  that  the  central 

limit  theorem  can  be  invoked  to  obtain  a  va.id  estimate  of  the  probability 

of  error  at  high  values  of  P  ,  where  the  range  of  validity  of  this  esti- 

e 

mate  is  dependent  upon  the  size  of  N  (Ref.  *14],  Therefore,  the  calcu¬ 
lations  of  initial  interest  are  those  of  E[in  z4  J  and  var(iJn  z^,). 
However,  the  difficulties  involved  in  making  even  these  calculations  are 
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clear,  so  that  it  is  reasonable  to  consider  bounds  on  the  probability  of 

error  obtainable  in  closed  form.  One  such  bound  of  interest,  since  it 

tightens  with  increasing  P  (i.e.,  tightens  as  the  SNR  decreases),  is 

e 

given  as  follows:  Noting  that 


£n  z  > 


1 


with  equality  iff  =  1,  it  follows  that 


P 

e 


>  0>  . 


(3.103) 


(3.104) 


Invoking  the  central  limit  theorem  as  discussed  above  in  order  to  evaluate 

this  lower  bound  on  P  .  it  is  found  that 

e 


OO  00 

J  dr'lc  J 


fm 

«u¥r 


2n  ll  +  2n  uJ 
‘ic  ^ic  'is  ^is 


r  2  2  ,2 

Pic  +  ^is  +  s 


,(m+4)/2 


(3.105) 


where  use  has  been  made  of 


(V  = 


mg 

2jt 


m 


[W 


~2  2  i(m+2  )/2 

+  1  J 


=  Pi,.  ,n  <\c’  \s)  '  (3*106) 

ic  is 


in  which  r,.  =  +  iqi£j,  ^  =  pic  +  ipis,  and  g2  *  m  N0/bq  At.  [it 

is  noted  that  Eq.  (3.106)  is  obtained  by  setting  N  =  1  in  P^(jl),  which 
is  given  by  Eq.  (3.9l)  for  the  special  case  m  =  1.]  Taking  now  the  spe¬ 
cial  case  0=3,  since  it  is  the  cose  of  most  interest, 
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Similarly, 


so  that 


Therefore,  since  the  z^  are  statistically  independent, 


i  Rj  R  [#)2  *  (¥)4.  • 


Now,  knowing  that 


(3.109) 


(x  - 1)  “  s  (¥)  +n(¥)  •  (3-no) 


(3.111) 


R) 


becomes  normally  distributed  as  N  »,  an  estimate  of  the  lower  bound 
given  by  Eq.  (3,104)  can  be  obtained  as  follows: 


CO 

I  s  f— 4 
5=3  o  K) 


(x  - 

exp - 5 -  dx  , 


(3.112) 
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where  p  is  given  by  £q.  (3.108)  and  cr  is  given  by  Eq.  (3,lll),  Thus 

2 *  Z 


>  erfc 


'9=3 


i=l 


V 


12 


i  [■(¥)’•  (¥)1 


i=l 
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(3.113) 


Since  this  bound  was  obtained  by  applying  the  central  limit  theorem  to  the 
sum  of  N  independent  random  variables,  it  is  well  known  that  the  result¬ 
ing  error  is  small  if  attention  is  restricted  to  those  situations  where  P 

l/o  e 

is  sufficiently  large— i.e.,  greater  than  about  C/N  '  ,  where  C  is  a 
proportionality  constant  [Kef.  44).  Thus  for  t.he  finite  N  of  interest, 
there  is  negligible  error  in  the  result  given  by  Eq.  (3.113)  when  thi_  re¬ 
sult  is  specialized  to  the  low  SNR  ca3e;  i.e.,  there  is  negligible  error 
when 


* 


>5 


and  there  is  obtained  the  asymptotic  result 


where  use  has  been  made  of  the  fact  that 


iklf 


(3.114) 
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However,  in  addition  to  the  fact  that  the  accuracy  of  the  result  obtained 
using  the  central  limit  theorem  improves  as  jffcJj/S  -*  0  for  fixed  N,  it 
is  also  true  that  the  lower  bound  on  in  z used  to  obtain  Eq,  (3.113) 
tightens  as  |pj  /§  -*  0,  so  that,  asymptotically,  we  have  the  low  SNR 


result 


lim  P  I  = 


(3.115) 


It  is  interesting  to  note  that,  as  expected,  this  low  SNR  result  is  depen¬ 
dent  only  upon  a  ratio  of  signal  energy  to  noise  energy  and  is  not  depen¬ 
dent  on  signal  shape. 

Now,  if  the  low  SNR  condition  |jyj|  «  £  is  not  satisfied,  Eq. 
(3.115)  does  not  hold  and  attention  must  be  returned  to  the  bound  given 
by  Eq.  (3.113).  However,  as  mentioned  above,  Eq.  (3.113)  is  not  a  useful 
bound  on  the  probability  of  error  at  low  values  of  P£  with  the  finite 
N  of  interest  here.  Therefore,  it  is  necessary  to  find  a  better  estimate 
of  P  for  the  high  SNR,  low  situation.  One  method  for  accomplishing 

this  is  to  apply  the  technique  known  as  distribution  "tilting"  [Ref.  45], 


Noting  that  the  results  found  for  low  P  using  the  central  limit  theorem 

e 

are  poor  because  of  the  fact  that  it  is  the  tail  of  the  probability  dis¬ 
tributed  that  must  be  evaluated,  the  tilting  procedure  will  be  applied 
here  with  the  purpose  of  moving  the  mean  of  the  tilted  distribution  to  the 
point  at  which  the  cumulative  probability  is  to  be  calculated.  Having 
accomplished  this,  the  desired  cumulative  probability  can  be  determined 
accurately  using  the  central  limit  theorem;  and,  in  fact,  an  exponentially 
correct  bound  on  this  probability  is  given  by  the  appropriate  Chernoff 
bound  [Ref.  41 ].  Now,  Eq.  (3.101)  states  that 


P  =  Pr 
e 


in  z >  0 


(3.101) 
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where 


Zi"  \\ -mJ2  +  I2  ' 


(3.102) 


Also,  there  is  the  Chernoff  bound  given  by 


Pr 


f  N 

in  >  0 


i=l 


l  <  exp[nN(s)  -  s  |^(s)] 


.(3.116) 


s=s. 


♦ 


where 


HN(s)  a  in  gN(s) 


(3.117) 


and 


N 

gN(s)  =  II  g.U)  , 

i=l 


(3.118) 


A 


in  which  the  g^s),  i  =  1,  . ..,  N,  are  the  moment-generating  functions 
of  the  Independent  random  variables  in  z^,  i  =  1,  ...»  N,  respectively. 
Furthermore,  the  point  at  which  the  right  hand  side  of  Eq.  (3.116) 

is  to  be  evaluated  specifies  the  amount  of  tilt  employed.  It  is  given 
here  by  the  value  of  s  >  0  which  satisfies  p^(s)  =  0»  where 

^(s)  ’  ^  ^  ■  (3-n9) 


Thus  for  the  case  of  interest  here,  there  is  obtained  the  bound 


(3.120) 
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Now,  although  the  Chernoff  bound  given  by  Eq,  (3.120)  is  an  upper  bound, 
it  can  be  shown  [Ref.  41 ]  through  application  of  the  Berry-Esseen  theorem 
[Ref,  44 J  that  this  bound  is  "exponentially  correct",  i.e.,  that 


N  N 

Pr.y  in  JjJO.aftf1/2)  n  «1(s)  |  , 

iti  1=1  *=ao 


where  f(N  is  a  polynomial  in  N  Therefore 


r1/2 


(3.121 


pe  »  f  (N  1//2)  n  gi(s )  j  , 
i=l  s=s0 


(3.122; 


where  f(N  1^2)  is  a  polynomial  in  N  and  sQ  >  0  is  that  value 

of  s  satisfying  jj  *  (s )  =  0.  Thus  the  calculation  of  P  reduces  to 

Ii  c 

the  determination  of  g  (s)  |  ,  i  =  1,  ....  N.  This  determination  is 

1  s=so 

accomplished  as  follows: 


gi(s)  =  E[exp(s  £n  z^]  =  E[z“j 

rH1)  £  r  f  [1 
■  T75T  *  J  J  7 


1 2  2 

+  ^ 


(2s-m-2)/2 


(3.123) 


Now,  in  order  to  simplify  the  calculations  the  following  assumptions  will 
be  made: 


Assumption  1.  Let  pis  =  0  V  i  -  1,  ...,  N,  so  that  the  known  signal 
m(t)  takes  the  form  m(t)  =  pc(t)  cos  o>gt.  This  assumption  is  not 
essentia]  but  greatly  simplifies  the  resultant  integrals. 

Assumption  2.  Take  0=3  (m  -  2),  since  this  case  is  easily  evaluated 
and  in  fact  is  the  case  of  most  interest. 
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With  these  assumptions, 


2 

‘<s,i8  =  ^  /  ax  I 


r  2  2  ,2, 

lx  +  y  +  £  ) 


0  -oo 


oc  r  /  \2  2  2  •» 

((*  -  nlc>  +  y  +  !  ) 


(3.124) 


where  sQ  >  0  must  satisfy  the  equation 


y  Ms^  i  =  o 

Z,  ds  i  SLS 


(3.125) 


That,  is,  sQ  >  0  must  be  the  solution  of 


7  7  ,2  2  ,2.S-2  r  2  2  ,2 

f  to  f  dy  ^  -*1  - - 

N  ^  Jr ,  ,2  2  2,  (x  -  U  )  .  y  +  r 

^  -oo  -oc  L(x  -  Pic)  +  y  +  5  J  [/  ic' 


—00  — oc 


00  00 

/  -  / 

-00  —00 


=  0  * 


t*2  4-  y2  *2J 
[(*  -h1c)2  +  y'  +  l2) 


(3.126) 


Conveniently,  it  can  be  shown  (see  Appendix  B)  that  Eq.  (3.126)  is  in 
fnct  uniquely  satisfied  by  s  =  1.  Thus,  plugging  this  fortunate  result 
into  Eq.  (3.124): 


2  00  CO 

,(*)  I  =  I  dx  I  dy  — - - 2 - 2 - " - 2 - 5 - 

1  sii  *  2  2  (x2  +  y  +  5  i[(x  -  J  +  y2  + 


I2]  ' 


(3.127) 

Now,  the  integration  with  respect  to  x  can  easily  be  done  as  a  contour 
integral  (see  Appendix  A),  and  it  is  found  that 
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This  integral  can  be  evaluated  using  Erdelyi  [Ref.  46],  item  6.2(35), 
which  gives 


(3.129) 


where  (a,  b;  c;  x)  is  Gauss*  hypergeometric  function.  Finally, 

substituting  Eq.  (3.129)  into  Eq.  (3.122), 


(3.130) 


where  it  is  recalled  that  N  =  2TB. 

At  this  point  it  should  be  mentioned  that  the  estimate  of 

given  by  Eq.  (3.130)  can  perhaps  be  further  refined  by  approximating  the 
—  l/2 

f(N  '  )  factor.  This  possibility  will  not  be  pursued  at  this  time,  how¬ 
ever,  since  Eq.  (3.130)  gives  the  essential  behavior,  particularly  at  low 
values  of  Pg.  Rather,  we  will  consider  optimization  of  the  exponentially 
correct  estimate  given  by  Eq.  (3.130)  with  respect  to  parameters  under  the 
control  of  the  design  engineer,  where  it  is  noted  that  this  estimate  is 
dependent  on  both  signal  shape  and  receiver  bandwidth  in  addition  to  sig¬ 
nal  energy  and  noise  energy.  In  order  to  proceed  with  this  optimization, 
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it  is  necessary  to  set  down  the  pertinent  constraints  on  the  transmitted 
signal  and  on  the  receiver  bandwidth.  These  can  be  specified  as  follows: 

1.  The  only  constraint  to  be  placed  on  the  transmitted  signal  in  the 
interval  [0,T]  is  an  average-power  constraint,  although  it  is 
noted  that  this  will  lead  to  the  same  result  as  if  a  peak-power 
constraint  were  also  imposed. 

2.  In  order  to  specify  the  constraint  on  the  receiver  bandwidth,  it 
is  noted  that  the  dependence  of  Eq.  (3.130/  on  the  receiver  band¬ 
width  enters  through  N  =  2TB.  This  dependence  follows,  of  course, 
from  the  signal  representation  used  in  this  case  (Case  Ha),  which 
assumes  that  the  atmospheric  noise  is  observed  in  an  rf  bandwidth 
2B.  This  assumption  stems  from  the  discussion  of  limited  spectrum 
availability  in  Section  IIB,  and  means  that  the  log-correlator  re¬ 
ceiver  shown  in  Fig,  23  is  in  practice  preceded  by  a  bandpass  filtei 
of  bandwidth  2B.  Therefore,  it  follows  that  the  pertinent  con¬ 
straint  on  the  receiver  bandwidth  is  that  it  be  small  enough  to 
suppress  strong  adjacent  channel  signals,  and  yet  large  enough  to 
pass  the  transmitted  signal. 

Proceeding  now  with  the  optimization  problem,  it  is  first  of  all  noted 
that  the  optimal  signal  shape  must  be  achieved  when  the  transmitted  signal 
energy  is  distributed  uniformly  in  a  portion,  say  T,  <  T,  of  the  avail- 

tv 

able  signaling  interval  of  length  T.  This  follows  directly  from  the 
Case  xla  representation  of  the  complex  envelope  of  the  received  signal  in 
terms  of  the  pulse  basis  functions  q^(t),  i  =  1,  . ..,  N.  With  this  re¬ 
sult  in  mind,  it  is  convenient  to  rewrite  Eq.  (3.130)  as 


’  y  b  /  y  b  \  j 

M1’ * f‘ 1  ‘ •  (3-131) 


where 


A  2  ^  ^ 

Qi  =  "ic  B  =  “2— 


(3.132) 


is  the  signal  power  received  under  h^  in  the  interval  At^,  and 


(3.133) 
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is  proportional  to  the  average  noise  energy  received  in  each  interval  of 
length  £t;  so  that 


v  = 


(3.134) 


is  proportional  to  the  average  noise  power  received  in  the  bandwidth  2B. 
These  identifications  are  seen  to  be  consistent  with  the  assumption  tha* 
in  practice  the  atmospheric  noise  can  be  considered  to  be  "white"  priox 
to  filtering  by  the  receiver. 

It  is  next  observed  that,  in  the  optimal  situation, 

V 


must  be  independent  of  the  value  of  T,  <  T,  since  both  Q,  and  B  . 

k  —  i  min 

are  directly  proportional  to  l/Tk  when  the  constraints  imposed  above 
are  satisfied.  Thus  it  follows  directly  that  the  optimal  choice  of  T^ 
is  given  by  T^  =  T.  Plugging  this  result  into  Eq.  (3.13l)  and  making 
use  of  the  fact  discussed  above  that  Q  =  Q  V  i  =  1,  ...»  2TB  is  opti- 

x  K 

mal,  there  results 


(3.135) 


Thus,  as  far  as  the  signal  shape  is  concerned,  it  is  concluded  that  the 
minimum  value  of  follows  from  distributing  the  signal  energy  uni¬ 

formly  in  the  signaling  interval  [o,T],  Furthermore;  as  expected,  P 

Q 

is  minimized  by  making  T  as  large  as  possible. 

It  now  remains  to  choose  the  optimal  receiver  bandwidth  2B.  It 
is  shown  in  Appendix  C  that,  in  fact,  B  ^  =  <».  This  result  is  disturbing, 
of  course,  since  the  noise  model  from  which  it  is  derived  has  been  demon¬ 
strated  to  be  valid  only  for  bandwidths  up  to  a  few  thousand  cycles  per 
second.  However,  examination  of  the  value  of  Pg  resulting  from  large  B 
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shows  that  the  estimate  given  by  Eq.  (3,135)  is  of  practical  interest, 
since  there  is  found  the  result  (see  Appendix  C): 

If  B  »  q/Yq,  then 


(3.136) 


Furthermore,  the  practical  usefulness  of  this  result  is  demonstrated  by  the 
fact  that  the  asymptote  is  approached  quite  rapidly  with  increasing  re¬ 
ceiver  bandwidth,  since,  e.g.,  B  =  3Q/Yg  gives 


(3.137) 


Examination  of  the  error  curve  given  by  Eq.  (3.137)  (see  Fig.  24)  shows 

~5 

that  error  rates  of  the  order  of  10  can  be  achieved  in  practice  using 
receiver  bandwidths  of  the  order  of  a  thousand  cycles  per  second.  Thus 
it  is  concluded  that  Eq.  (3.136)  gives  a  valid  estimate  of  the  error  per¬ 
formance  achievable  in  the  presence  of  additive  atmospheric  noise.  This 
estimate  is  conveniently  written  as 


(3.138) 


(l' 

where  E  =  QT  is  the  signal  energy  under  h'  '  distributed  uniformly 
in  the  signaling  interval.  This  Pg  result  is  plotted  vs  SNR  in 
Fig.  24,  along  with  the  error  curve  resulting  from  the  use  oi  a  matched- 
filter  receiver  in  the  presence  of  additive  generalized  "t”  noise.  In 
addition,  the  error  curve  achieved  by  a  noncoherent  FSK  system  operating 
at  vlf  is  also  presented  using  measured  data  obtained  by  Watt,  _et  £l 
[Ref.  47] . 
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In  Fig.  24,  the  parameter  A  =  V  /v  is  introduced  and  the 

rms'  avg 

error  curves  plotted  for  2  <  A  <  4  because  of  the  fact,  discussed  in 
detail  in  Chapter  II,  that  ■  V^s  Is  determined  primarily  by  the  tail  of 
the  probability  distribution,  whose  observation  is  made  difficult  by  the 
large  dynamic  range  exhibited  by  atmospheric  noise.  Thus,  noting  that 
the  noise  envelope  distributed  according  to  the  generalized  "t"  model 
actually  has  a  divergent  second  moment  for  Q  <  3,  the  parameter 
2  <  A  <  4  is  introduced  to  show  the  variation  in  the  error  curve  due  to 
the  range  ol  values  of  A  actually  observed  in  practice.^ 

4 .  Discussion  of  Probability-of -Error  Results  (Case  Ila) 

Before  proceeding  to  the  discussion  of  the  results  calculated 
above  it  is  recalled  that  these  Case  II  results  apply  to  long-duration 
signals,  where  long  duration  means  that  the  slowly  varying  modulating 
process  a(t)  is  allowed  to  vary  on  the  signaling  interval.  Furthermore, 
the  Case  Ila  results  apply  specifically  to  the  Case  II  situation  where 
the  complex  envelope  of  the  received  signal  can  be  described  by  Eq.  (3.76) 
in  terms  of  N  =  2TB  statistically  independent  samples.  It  is  proposed 
that  this  is  the  simplest  representation  of  Case  II  signals  of  practical 
interest.  It  is  seen  that  it  essentially  includes  the  Case  I  signals  as 
a  special  case.  With  this  in  mind,  inspection  of  the  error  curves  plotted 
in  Fig.  24  shows  that  the  performance  predicted  for  the  log-correlator 
receiver  in  the  presence  of  additive  atmospheric  noise  is  significantly 
better  than  either  that  predicted  for  the  matched-filter  receiver  or  that 
actually  achieved  by  the  FSK  system.  Although  the  error  curves  are  fairly 
self  explanatory,  several  comments  concerning  their  derivation  are  war¬ 
ranted: 


Note  that  the  square  of  the  nns  carrier-to-noise  ratio  plotted  on  the 
abscissa  in  Fig.  24  is  related  by  a  proportionality  constant  to  the 
signal-to-noise  energy  ratio  E/Yq.  This  constant  is  proportional  to 
the  ratio  of  signal  bandwidth  to  receiver  bandwidth.  For  the  case  in 
Fig.  24  (using  A  =  3)  it  is  about  5  X  10“3,  although  this  value  may 
be  too  small  by  a  factor  l/2.  This  error  would  result  in  the  performance 
of  both  the  matched-filter  receiver  and  the  log-correlator  receiver 
being  3  db  worse  than  that  shown  in  Fig.  24.  This  cannot  be  resolved 
due  to  the  fact  that  E^g  as  given  by  Watt,  et  £l  in  their  Fig.  18  does 
not  appear  to  be  consistent  with  the  scale  on  the  abscissa. 


129 


SEL-66-052 


t  .‘ijiffltfrirrtinwnmim,  „  ,  ,  . ! 


20  40  60 

)-N0ISE  RATIO  IN  A  I  KC 
BAND  -  db 


FIG.  24.  RECEIVER  PERFORMANCE  IN  PRESENCE  OF  ADDITIVE  ATMOSPHERIC  NOISE: 
COMPARISON  OF  PREDICTED  ERROR  CURVES  WITH  MEASURED  PERFORMANCE  OF  AN 
FSK  SYSTEM  OPERATING  AT  VLF. 


Specification  of  transmitted  signal: 

Frequency  shift  keyed,  ±  25  cps  shift 
60  WPM  start-stop  teletype 


[Note:  Vertical  scale  is  log| log  P  |  or  log] log  Pq(v)| 
where  Pq(v)  *s  t^ie  probability  that  the  envelope 
V  exceeds  the  value  specified  by  abscissa.] 


LEGEND: 


Measured  data:  Complement  of  the  probability  dis¬ 
tribution  function  of  the  envelope  of  the  received 
vlf  noise.  (Watt,  et  _al,  [Ref.  47]) 

Receiver  6-db  bandwidth  =  120  cps 


Measured  data:  Noncoherent  receiver  binary  error 
curve  (derived  by  Watt,  £t  al  from  measured  character 
error  rates). 

Receiver  specifications: 

120  cps  IF  bandwidth 
70  cps  base-bandwidth 


Predicted  matched-filter  receiver  error  curves 
(matched  to  mark  frequency  only). 

Generalized  "t"  noise  model,  9=3 


■  /v 

=  2 

rms' 

avg 

■  /V 

=  4 

rms 

avg 

(see  page  129) 


Predicted  log-correlator  receiver  error  curves  (coher¬ 
ent  detection  of  mark  frequency,  space  frequency  re¬ 
moved  ) 

Generalized  "t"  noise  model,  9=3 


V  /V  =2 

rms'  avg 

V  /V  =4 

rms  avg 


(see  page  129) 
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Log-correlator  receiver:  The  error  curve  presented  for  the  log- 
correlator  receiver  was  obtained  by  using  the  <ponentially  correct 
estimate  given  by  Eq.  (3.138)  in  conjunction  with  the  low  SNR  result 
given  by  Eq.  (3.115).  The  low  SNR  result  gives  a  good  estimate  of 
the  probability  of  error  at  high  values  of  Pe.  This  portion  of  the 
error  curve  was  obtained  by  writing  Eq.  (3.115)  in  terms  of  the  no¬ 
tation  introduced  in  Eq.  (3.13l).  This  gives  the  convenient  form: 


(3.139) 


Note,  as  mentioned  earlier,  that  this  low  SNR,  high  Pe  result  de¬ 
pends  only  upon  a  ratio  of  signal  energy  to  noise  energy  and  is 
independent  of  signal  shape.  In  contrast  to  this,  the  estimate  of 
Pe  given  by  Eq.  (3.138),  which  is  exponentially  correct  at  all 
values  of  Pe,  was  shown  above  to  depend  crucially  on  signal  shape. 
Now,  it  is  interesting  to  investigate  the  claim  that  the  performance 
predicted  by  Eq.  (3.138)  is,  in  fact,  guaranteed  only  if  the  signal 
energy  is  distributed  uniformly  in  the  signaling  interval,  Obser¬ 
vation  of  the  result  that  B  =  <»  is  optimum  is  actually  seen  to 
discount  this  claim  when  the  only  constraint  on  th-j  transmitted  sig¬ 
nal  is  an  average  power  constraint,  since  a  given  Pe  performance 
can  be  obtained  (when  the  signal  energy  is  not  uniformly  distributed 
in  the  signaling  interval)  merely  by  sufficiently  increasing  the 
receiver  bandwidth  2B.  However,  it  must  be  remembered  that  chese 
results  are  obtained  using  the  Case  Ila  assumption  that  the  samples 
of  the  complex  envelope  of  the  received  noise  are  statistically  in¬ 
dependent — an  assumption  that  certainly  breaks  down  for  receiver 
bandwidths  larger  than  about  104  cps.  Thus  it  is  concluded  that, 
practically  speaking,  the  optimal  situation  is  indeed  that  where  the 
signal  energy  is  uniformly  distributed  in  the  signaling  interval, 
since  this  achieves  a  given  level  of  performance  with  the  smallest 
receiver  bandwidth. 

2.  Matched-filter  receiver:  The  error  curve  presented  for  the  matched- 
filter  receiver  was  obtained  by  noting  that  if  attention  is  re¬ 
stricted  to  a  class  of  signals  whose  matched  filter  can  be  assumed 
to  be  a  narrowband  filter  of  bandwidth  l/T,  then  the  calculation 
of  the  probability  of  error  reduces  to  a  one-dimensionai  calculation 
[Ref.  37]  similar  to  that  required  in  Case  I.  It  may  be  noted  that 
this  class  of  signals  contains  many  signals  of  interest,  including 
the  FSK  signal  used  in  the  system  studied  by  Watt,  et_  al. 

Now,  as  in  Case  I,  the  high  SNR  case  is  the  most  interesting  and 
gives  the  result 
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Thus  the  improvement  possible  using  the  log-correlator  receiver 
instead  of  the  matched-filter  receiver  in  the  presence  of  additive 
generalized  "t"  noise  is  clear.  This  is  demonstrated  in  Fig.  24 
where  the  entire  erxor  curve  is  plotted. 

3.  Noncoherent  FSK  system:  The  error  curve  presented  for  the  FSK  sys¬ 
tem  was  derived  by  Watt,  et  a_l  from  measured  character  error  rates. 
Once  again,  the  advantage  to  be  gained  by  using  a  log-correlator 
receiver  instead  of  the  receiver  actually  employed  is  obvious.  Al¬ 
though  all  of  the  details  of  operation  of  this  receiver  are  not 
presented  by  Watt,  it  is  one  of  a  class  of  receivers  often  used  in 
practice  and  will  be  considered  in  more  detail  in  Chapter  IV. 

The  fact  that  the  predicted  performance  of  the  log-correlator 
receiver  is  much  better  than  that  predicted  for  the  matched-filter  receiver 
indicates  the  necessity  of  using  nonlinear  receiving  techniques  in  the 
presence  of  additive  atmospheric  noise.  In  fact,  this  is  the  procedure 
used  in  practice  in  conjunction  with  long-duration  signals.  However,  the 
relatively  poor  performance  of  the  FSK  system  studied  by  Watt,  et  al^ 
indicates  that  some  of  the  receiving  techniques  used  in  practice  are  far 
from  optimum.  For  example,  for  the  60  wpm  start-stop  teletype  system  con¬ 
sidered  in  Fic.  24,  the  plotted  error  curves  indicate  that  an  error  rate 
4 

of  1  error  per  10  bits  can  be  achieved  using  an  appropriate  log-correlator 
receiver  with  an  SNR  that  is  20  to  30  db  less  than  the  SNR  required  by  the 
FSK  system  analyzed  by  Watt,  et  _al.  This  means  that  there  is  an  order-of- 
magnitude  difference  in  the  transmitter  power  required  by  the  two  systems 
in  order  to  reach  this  particular  level  of  performance;  and  this,  of  course, 
has  strong  economic  overtones.  As  a  result,  it  is  of  Interest  to  apply 
the  generalized  "t"  model  in  an  attempt  to  predict  the  probability  of 
error  resulting  from  the  use  of  the  various  receiving  techniques  that  have 
been  proposed  for  use  in  additive  atmospheric  noise.  For  example,  the 
"smear-desmear"  technique  mentioned  earlier  can  be  investigated  using  the 
generalized  "t"  model.  It  appears  from  cursory  examination  chat  this 
technique  leads  to  considerably  poorer  error  performance  than  that  calcu¬ 
lated  for  the  log-correlator  receiver.  This  follows  from  the  fact  that 
the  performance  predicted  by  Eq,  (3.138)  for  the  log-correlator  receiver 
is  significantly  superior  to  the  performance  of  an  optimal  matched-filter 
receiver  in  the  presence  of  additive  gaussian  noise  with  the  same  average 
noise  power.  Quantitatively,  this  superiority  is  equivalent  to  12  to 
15  db  in  input  SNR;  and  the  fact  of  this  superiority  is  consistent  with 
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the  notion  that  the  additive-atmospheric-noise  channel  has  a  higher  ca¬ 
pacity  than  the  additive-gaussian-noise  channel  having  the  same  average 
noise  power. 

This  investigation  of  proposed  receiving  techniques  will  be  pur¬ 
sued  further  in  Chapter  IV,  where  proposed  nonlinear  receiving  techniques 
will  be  considered. 

E.  SUMMARY  AND  CONCLUSIONS 

The  work  in  this  chapter  has  been  concerned  with  the  detection  of 
known  signals  in  additive  atmospheric  noise.  In  particular,  detailed 
consideration  has  been  given  to  the  problem  of  deciding  between  the  two 
&  priori  equiprobable  hypotheses 


h 


x(t)  =  m(t)  +  y(t) 


x(t)  =  y(t) 


0  <  t  <  T 


(3.2) 


with  the  smallest  probability  of  error.  The  additive  atmospheric  noise 
y(t)  is  represented  in  this  work  by  the  generalized  "t"  model  developed 
in  Chapter  II.  Two  cases  are  identified: 

1.  Case  I:  Short-duration  signals.  This  case  is  characterized  by  the 
fact  that  the  slowly  varying  modulating  process  a(t)  in  the  gen¬ 
eralized  "t"  model  con  be  assumed  constant  for  the  duration  of 
the  transmitted  signal.  The  correlator  receiver  described  by  Eq. 
(3.23)  is  found  to  be  optimum  in  this  case,  and  the  probability  of 
error  achieved  by  this  receiver  is  calculated  in  Eq.  (3,44).  The 
probability  of  error  for  high  SNR  is  the  result  of  most  interest, 
and  is  given  by 


(3.45) 


noise 

energy  and  E,^  is  the  energy  of  the  known  signal  tnU)  at  the  out¬ 
put  of  the  whitening  filter.  The  slow  decrease  with  increasing  sig¬ 
nal  energy  of  this  Pe  result,  relative  to  the  exponential  decrease 
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where  yw  is  proportional  to  the  average  received  "whitened" 
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achieved  in  the  presence  of  additive  goussinn  noise  (see  Fig.  22), 
indicates  that  short-duration  signals  are  not  an  optimal  choice  in 
the  presence  of  atmospheric  noise.  This  is,  of  course,  physically 
v  reasonable  because  of  the  distinctive  impulsive  nature  of  atmo¬ 

spheric  noise.  This  suggests  the  use  of  signals  of  long  duration 
relative  to  that  of  the  typical  noise  pulse. 

In  addition  to  these  short-duration  signal  results,  the  Case  I  as¬ 
sumption  was  also  applied  to  the  situation  (Case  Im)  where  the 
additive  noise  is  described  by  the  modified  generalized  "t"  model. 
The  correlator  receiver  given  by  Eq.  (3.23)  continues  to  be  optimum 
in  this  case,  but  the  probability  of  error  is  now  given  approximately 
in  the  high  SNR  case  by 


(3.53) 


This  result  reduces  in  form  to  that  given  by  Eq.  (3.45)  above  when 
(E^/p2)  «  1  is  satisfied,  which  is  the  case  of  practical  interest 
at  vlf  and  If  frequencies  of  operation.  On  the  other  hand. 

(3.53)  predicts  an  exponential  decrease  in  Pe  with  increas-..0  sig¬ 
nal  energy  when  (E-^/p2)  »  1  is  satisfied;  this  condition  may  be 
approached  in  practice  at  hf.  Thus  it  is  concluded  that  short- 
duration  signals  are  more  attractive  for  use  at  hf  than  at  lower 
frequencies  of  operation.  This  is  consistent  with  the  fact  that  the 
received  noise  loses  its  impulsive  appearance  as  the  frequency  of 
observation  increases. 

2.  Case  II:  Long-duration  signals,  the  general  case.  This  case  is 
characterized  by  the  fact  that  the  slowly  varying  process  a(t)  is 
now  allowed  to  vary  on  the  signaling  interval.  In  particular, 

Case  Ila  generalized  "t"  noise  is  considered  in  detail  and  is 
specified  to  be  that  Case  II  noise  whose  complex  envelope  can  be 
described  by  Eq.  (3.76)  in  terms  of  2TB  statistically  independent 
samples.  This  is  proposed  as  the  simplest  representation  of  Case  II 
signals  of  practical  interest,  and  the  optimal  receiver  in  this  case 
is  found  to  be  a  "logarithmic-correlator"  receiver  (see  Fig.  23) 
which  implements  the  rule 

(l' 

Choose  h'  '  iff 


T  T 

J  ^n [ | x ( t ) j  2  +  72]dt  >  j  £n [ | x ( t )  -  p(t)J2  +  ?2Jdt  ,  (3.97) 
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where  |x(t)|  is  the  envelope  of  x(t),  |x(t)  -  p(t)|  is  the 
envelope  of  x(t)  -  m(t)  and  y ^  is  proportional  to  the  average 
received  noise  power.  An  "exponentially  correct"  estimate  of  the 
probability  of  error  achieved  by  this  receiver  when  the  transmitted 
signal  energy  is  distributed  uniformly  in  the  signaling  interval 
is  given  by 


lim  P  |  »  f  (l/y/ira)  exp 1-y5  E1  ,  (3.138) 

YqB  e0“3  L  0  . 


where  f  ( l/^TB)  is  a  polynomial  in  (l/y/2TB) ,  E  =  QT  is  the 
energy  of  the  known  signal  m(t),  and  Yq  =  2(^2/b)  is  propor¬ 
tional  to  the  average  received  noise  energy.  The  exponential  be¬ 
havior  of  this  Pe  result  is  strikingly  superior  to  that  given  by 
Eq.  (3.45)  for  short-duration  signals.  Furthermore,  it  is  shown 
(see  Fig.  24)  that  this  performance  of  the  log-correlator  receiver 
is  significantly  better  than  either  that  achievable  by  a  matched- 
filter  receiver  or  that  actually  achieved  oy  a  typical  FSK  vlf 
communication  system.  Also,  to  demonstrate  that  the  asymptotic  per¬ 
formance  predicted  by  Eq.  (3.138)  can  be  closely  approached  in  prac¬ 
tice,  the  special  case  YqB/Q  =3  is  noted.  It  is  shown  in  Section  • 

D3  that  this  special  case  is  reasonably  achieved  in  practice,  and 
that  it  gives  the  result 


(3.141) 
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With  the  above  probability-of -error  results  in  mind,  it  is  inter¬ 
esting  to  note  that  there  is  an  alternate  realization  of  the  log- 
correlator  receiver  given  by  Eq.  (3.98)  that  is  approximately 
equivalent  to  Eq.  (3,97)  when  the  condition  (YqB)/(4Q)  »  1  is  sat¬ 
isfied.  This  realization  takes  the  form: 

Choose  h^^  iff 


I  f _ *(il  h*(.U  . 
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dt  V  > 
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This  receiver  rule  is  of  interest,  since  it  may  be  simpler  to  imple¬ 
ment  than  Eq,  (3.97),  and  since  it  has  a  realization  closely  related 
to  the  matched-filter  realization  of  the  correlator  receiver  rule. 
However,  it  is  important  to  note  that  the  condition  (YqB)/(4Q)  »  1 
is  a  low  signal-to-noise  power  ratio  condition  whose  satisfaction  in 
practice  may  introduce  more  system  complexity  than  it  removes,  given 
a  desired  level  of  system  performance. 

Finally,  it  is  concluded  that  the  error  performance  predicted  by  Eq.  (3.138) 
for  the  log-correlator  receiver  in  the  presence  of  additive  atmospheric 
noise  is  of  particular  interest,  since  it  is  the  only  theoretical  estimate 
known  to  the  author  of  the  performance  achievable  in  the  presence  of  addi¬ 
tive  atmospheric  noise.  Furthermore,  the  simplicity  of  the  log-correlator 
receiver  indicates  that  this  level  of  performance  can  be  approached  by 
practically  realizable  receivers. 


IV.  RELATED  TOPICS 


A.  COMPARISON  OF  THE  LOG-CORRELATOR  WITH  ANOTHER  RECEIVER  FORM 

In  the  Case  II  "known  signal"  considerations  in  Chapter  III,  it  was 
shown  that  the  receiver  which  is  optimal  in  tne  presence  of  additive  at¬ 
mospheric  noise,  in  the  sense  that  its  use  minimizes  the  probability  of 
error,  is  the  log-correlator  receiver.  During  the  discussion  of  this 
Case  II  result,  the  desirability  of  comparing  the  performance  predicted 
for  the  log-correlator  receiver  with  that  obtainable  using  other  receiving 
techniques  was  mentioned.  In  particular,  it  is  desired  to  compare  its 
performance  with  that  of  receivers  commonly  used  in  practice.  Inspection 
of  the  log-correlator  receiver  (Fig.  23)  leads  to  several  qualitative  con¬ 
clusions  on  this  comparison.  It  is  seen  that  the  log-correlator  receiver 
performs  essentially  four  operations  as  it  processes  the  received  signal. 
In  the  order  of  their  performance,  these  operations  can  be  described  as 
follows: 

1.  The  receiver  makes  use  of  the  "known"  nature  of  the  transmitted 
signal,  i.e.,  the  received  signal  is  compared  with  each  of  the 
possible  transmitted  signals,  each  of  which  is  known, 

2.  The  results  of  the  comparisons  in  (l)  are  nonlinearly  processed  in 
a  fashion  which  suppresses  the  effect  on  the  resultant  decision  of 
the  comparisons  at  those  times  when  the  received  signal  is  largely 
the  result  of  a  pulse  of  noise. 

.3.  The  results  of  the  zero-memory  nonlinear  processing  of  step  (2)  are 
summed  over  the  signaling  interval. 

4.  The  resultant  sums  in  step  (3)  are  compared  against  one  another  .nd 
a  decision  made  on  the  basis  of  this  comparison. 

It  is  important  to  note  that  the  Case  II  considerations  of  Chapter  III 
make  the  assumption  that  the  received  signal  has  been  appropriately 
bandpass-filtered  prior  to  step  (l)  above.  Several  factors  enter  into 
the  choice  of  this  predetection  bandwidth,  which  is  discussed  in  detail 
in  Chapter  III;  but  in  practice  the  selection  procedure  can  be  summarized 
as  follows:  At  operating  frequencies  where  atmospheric  noise  is  impor¬ 
tant,  the  predetection  bandwidth  should  be  chosen  as  large  as  possible 
consistent  with  the  exclusion  of  strong  interfering  narrowband  signals. 
This  means  that  the  limited  spectrum  availability  in  the  frequency  range 
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of  interest  is  typically  the  factor  that  determines  the  predetection  re¬ 
ceiver  bandwidth. 


The  receiver  form  whose  performance  will  be  compared  here  with  the 
performance  of  the  log-correlator  receiver  is  that  shown  in  Fig.  25.  This 
receiver  is  of  interest  because  of  its  simplicity  of  construction,  which 
makes  it  perhaps  the  most  commonly  used  receiver  form  in  the  presence  of 
additive  atmospheric  noise.  Typical  zero-memory  nonlinear  devices  used 
in  the  configuration  shown  in  Fig.  25,  depending  upon  the  type  of  modula¬ 
tion  employed  at  the  transmitter,  are: 

1.  The  wideband  "clipper"  which  operates  to  limit  the  amplitude  of  the 
received  signal  when  this  signal  amplitude  exceeds  a  specified  level 
indicating  that  it  is  predominantly  the  result  of  a  pulse  of  noise. 

2.  The  wideband  "limiter",  often  used  with  angle  modulation  schemes, 
which  provides  an  output  signal  whose  amplitude  is  fixed  independent 
of  the  amplitude  of  the  input  signal,  and  whose  sign  is  the  same  as 
that  of  the  input  signal. 


DECISiOW 


FIG.  25.  BLOCK  DIAGRAM  OF  RECEIVER  COMMONLY  USED  IN  PRESENCE  OF 
ADDITIVE  ATMOSPHERIC  NOISE. 
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Comparison  of  the  performance  of  the  receiver  form  described  in 
Fig.  25  with  the  performance  previously  outlined  for  the  log-correlator 
receiver  proceeds  qualitatively  as  follows:  Noting  that  the  predetec¬ 
tion  filtering  required  in  conjunction  with  the  log-correlator  receiver 
is  also  required  in  practice  by  the  receiver  described  by  Fig.  25,  it  is 
seen  that  the  two  receivers  perform  similar  operations  on  the  received 
signal,  with  the  exception  that  the  order  in  which  operations  (l)  and  (2) 
are  performed  is  interchanged.  In  addition  to  the  fact  that  the  two  re¬ 
ceivers  employ  different  zero-memory  nonlinear  devices,  it  is  clear  that 
the  receiver  described  by  Fig.  25  must  be  inferior  to  the  log-correlator 
receiver  because  of  this  interchange  of  operations.  This  conclusion  fol¬ 
lows  from  the  fact  that  the  receiver  described  by  Fig.  25  destroys  part 
of  the  information  known  about  the  possible  transmitted  signals  before 
it  makes  use  of  this  information,  and  thus  must  perform  less  well  than 
the  optimal  log-correlator  receiver  in  the  "known  signal"  situation  being 
considered  here.  Experimental  support  for  this  conclusion  was  given  in 
Fig.  24,  where  the  performance  predicted  for  the  log-correlator  receiver 
is  seen  to  be  far  superior  to  that  actually  achieved  by  an  FSK  communica¬ 
tion  system  operating  at  vlf .  Although  Watt,  et  £l  do  not  give  all  of 
the  details  of  the  particular  FSK  receiver  employed,  they  indicate  that 
it  basically  has  the  form  described  in  Fig.  25,  where  a  limiter  performs 
the  zero-memory  nonlinear  operation. 

B.  DETECTION  AND  ANALYSIS  OF  WHISTLER-MODE  SIGNALS 

In  addition  to  the  "known  signal"  detection  problem  considered  so  far, 
it  is  of  interest  to  consider  the  application  of  the  generalized  "t" 
model  to  other  signal  analysis  problems  in  the  presence  of  additive  atmo¬ 
spheric  noise.  One  such  problem  is  the  estimation  problem  in  which  it  is 
desired  to  find  the  "best"  (in  some  sense)  estimate  to  a  transmitted  sig¬ 
nal  that  is  not  completely  known  £  priori .  given  observations  at  the 
receiving  site— i.e.,  given  observations  only  after  the  atmospheric  noise 
has  been  added.  As  on  example  of  practical  importance,  we  will  consider 
the  problem  of  detecting  the  presence  of  signals  propagating  in  the 
whistler  mode  (Ref.  48]  and  in  addition  of  measuring  the  defining  charac¬ 
teristics  of  the  detected  whistler-mode  signals.  Although  no  analytical 
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results  are  as  yet  available  on  this  estimation  problem  in  the  presence 
of  additive  atmospheric  noise,  it  appears  worthwhile  to  note  the  following: 
The  whistler-mode  signals  which  it  is  desired  to  analyze  are  described  in 
detail  by  Helliwell  [Ref.  48],  and  can  be  described  for  the  purposes  of 
this  discussion  as  vlf  signals  possessing  a  very  large  time-bandwidth 
product,  i.e.,  bandwidths  of  the  order  of  10  kc  and  durations  of  the  order 
of  one  second.  Furthermore,  the  name  "whistler"  stems  from  the  fact  that 
these  signals  typically  consist  of  a  gliding  tone  which  sweeps  across  the 
vlf  band  from  high  to  low  frequencies  producing  an  audible  "whistle," 


Analysis  of  Typical  Whistler  Receivers 


The  purpose  of  this  investigation  of  the  whistler  estimation  prob¬ 
lem  is  twofold:  to  consider  the  usefulness  of  the  generalized  "t"  model 
in  formulating  a  solution  to  this  problem,  and  to  use  the  insights  gained 
in  the  known-signal  problem  for  evaluating  the  performance  of  receiver 
forms  presently  used  in  whistler  analysis.  Thus  in  order  to  proceed,  the 
receiver  form  typically  used  to  study  whistlers  in  the  presence  of  additive 
atmospheric  noise  is  presented  in  Fig.  26.  The  zero-memory  nonlinear  de¬ 
vice  normally  used  in  this  configuration  is  the  wideband  clipper  described 
in  connection  with  Fig.  25,  so  that  this  receiver  will  henceforth  be  called 


FIG.  26.  BLOCK  DIAGRAM  OF  CLIPPER  RECEIVER  COMMONLY  USED  TO 
ANALYZE  WHISTLER  SIGNALS. 


a  "clipper  receiver."  In  order  to  proceed  with  the  analysis  of  this  re¬ 
ceiver,  it  is  first  of  all  noted  that  the  Case  II  known-signal  results 
cannot  be  applied  directly,  since  the  information  known  about  the  whistler 
and  of  course  that  utilized  by  the  clipper  receiver  is  not  sufficient  to 
make  the  whistler  a  known  signal.  Rather,  we  can  conclude  that  if  the 
whistler  to  be  analyzed  can  be  further  specified  so  that  the  problem  be¬ 
comes  a  known-signal  problem,  then  the  log-correlator  receiver  is  the 
optimal  receiver.  Otherwise  (as  is  usually  the  case  since  it  is  the  de¬ 
tailed  measurement  of  whistler  parameters  that  is  of  scientific  importance) 
the  optimal  receiver  is  not  obvious,  and  the  clipper  receiver  is  a  candi¬ 
date  whose  use  can  be  supported  on  physical  grounds.  In  order  to  show 
this,  it  is  noted  that  whereas  the  whistler  is  characterized  by  both  a 
large  bandwidth  and  long  duration,  the  received  whist’.er  signal  is  con¬ 
taminated  by  additive  noise  which  primarily  comes  from  two  distinct  sources 

1.  Additive  atmospheric  noise:  This  noise,  which  is  discussed  in  detail 
in  Chapter  II,  has  a  bandwidth  of  the  order  of  10  kc  (no  narrowband 
assumption  here)  and  can  be  characterized  as  a  sequence  of  noise 
pulses,  each  of  which  has  a  duration  of  much  less  than  one  second. 

2.  Narrowband  interfering  signals:  In  contrast  to  the  communication 
situation  where  a  narrowband  assumption  is  reasonable,  the  wide  band¬ 
width  of  the  whistler  means  that  strong  narrowband  interfering  com¬ 
munication  signals  may  exist  within  the  frequency  band  of  interest. 

With  this  description  of  the  additive  noise  in  n  .nd,  the  perfor¬ 
mance  of  the  clipper  receiver  can  be  investigated  by  examining  the  opera¬ 
tions  performed  by  the  individual  receiver  components.  This  examination 
shows  that  in  fact  the  clipper  receiver  is  a  good  candidate  for  use  in  the 
presence  of  the  additive  noise  described  above,  since  it  suppresses  the 
effect  of  the  noise  on  the  spectrum  analyzer  output  in  the  following  ways: 

1.  The  wideband  clipper  serves  to  reduce  the  effect  of  atmospheric 
noise  pulses  before  these  pulses  can  be  "smearea"  in  time  by  the 
narrowband  filters  comprising  the  typical  spectrum  analyzer.  It  is 
clear  that  the  usefulness  of  the  clipper  depends  upon  the  ability 
of  the  operator  to  set  the  clip  level  close  to  the  maximum  level  of 
the  received  whistler.  This  ability  is  a  function  of  knowledge  both 
of  the  level  of  the  whistler  signal  and  of  the  level  of  narrowband 
interfering  signals.  (The  problem  of  optimally  setting  the  clip 
level  is  discussed  further  below,  ) 

2.  The  effect  of  narrowband  interfering  signals  is  suppressed  by  the 
fact  that  they  contaminate  only  a  small  portion  of  the  frequency 
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band  examined  by  the  spectrum  analyzer.  Thus  the  observer  at  the 
output  of  the  spectrum  analyzer  is  able  to  pick  out  the  whistler 
signal  by  simply  ignoring  the  spectral  lines  produced  by  the  narrow- 
band  interference. 

It  is  important  to  note  that  although  the  above  discussion  indi¬ 
cates  that  narrowband  interfering  signals  are  not  directly  a  problem,  they 
may  be  an  important  factor  in  preventing  an  optimum  setting  of  the  clip 
level.  Thus  it  is  seen  that  the  actual  interfering  signal  environment  is 
important  in  the  analysis  of  the  performance  of  the  clipper  receiver,  and, 
of  course,  is  also  important  in  any  analytical  formulation  of  the  whistler 
estimation  problem.  In  fact,  the  above  discussion  dictates  that  the  clip¬ 
per  receiver  performs  best  when  strong  narrowband  interfering  signals  are 
removed  by  selective  linear  filtering  prior  to  the  performance  of  the  clip¬ 
ping  operation.  Having  removed  these  strong  interfering  signals  ("strong" 
refers  here  to  the  situation  where  the  received  interfering  signal  power 
is  significantly  greater  than  that  of  the  received  whistler),  the  perfor¬ 
mance  of  the  clipper  receiver  becomes  a  function  of  how  well  the  level  of 
the  whistler  signal  to  be  analyzed  is  known.  In  fact,  if  this  level  is 
known  quite  precisely,  and  if  the  setting  of  the  clip  level  Just  above  this 
known  signal  level  results  in  clipping  for  only  a  small  percentage  of  the 
time,  then  the  clipper  appears  to  be  a  good  choice  of  nonlinearity  that  is 
effective  in  suppressing  the  additive  atmospheric  noise.  On  the  other  hand, 
if  the  whistler  signal  level  is  known  only  approximately  (perhaps  only  sta¬ 
tistically),  or  if  the  setting  of  the  clip  level  just  above  the  known  sig¬ 
nal  level  results  in  clipping  a  large  percentage  of  the  time,  then  it 
appears  that  there  must  certainly  be  another  nonlinearity  which  performs 
better  than  the  wideband  clipper. 

2.  Conclusions 

The  foregoing  discussion  has  been  presented  to  specify  an  estimation 
problem  of  practical  interest  in  the  presence  of  additive  atmospheric  noise. 
The  clipper  receiver  typically  used  in  the  analysis  of  whistlers  has  been 
examined  in  some  detail,  and  is  seen  to  operate  on  the  received  signal  in 
a  fashion  that  agrees  both  with  intuition  and  with  the  results  of  the  known- 
signal  detection  problem.  However,  the  discussion  also  indicates  that  the 
clipper  receiver  is  certainly  suboptimal  if  either  the  power  level  of  the 


A 

whistler  signal  is  unknown  or  the  whistler  signal  is  a  very  weak  signal — 
i.e.,  if  there  is  a  very  low  SNR  characterized  by  clipping  for  a  large 
percentage  of  the  time.  Since  usually  the  first  and  often  both  of  these  ¥ 

conditions  are  characteristic  of  the  whistler  estimation  problem  in  prac¬ 
tice,  it  Is  proposed  that  improved  receiving  techniques  may  result  from 
applying  the  generalized  "t"  model  to  the  whistler  estimation  problem. 

This  possibility  will  be  considered  further  in  Chapter  V. 
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V.  SUMMARY  AND  CONCLUSIONS 


— 1 v  - 


A.  THE  GENERALIZED  "t"  MODEL  FOR  RECEIVED  ATMOSPHERIC  NOISE 

This  work  has  been  concerned  with  the  development  and  application  of 
an  analytical  model  for  impulsive  phenomena.  In  particular,  the  work  has 
focused  on  a  specific  impulsive  phenomenon,  atmospheric  raaio  noise. 

In  Chapter  II  a  new  model,  called  the  generalized  "t"  model,  was 
developed  and  was  verified  to  be  appropriate  as  a  description  of  received 
atmospheric  noise,  i.e.,  atmospheric  noise  as  observed  through  a  given 
receiver  passbana.  This  generalized  "t"  model  describes  the  received 
atmospheric  noise  y(t)  as 


y(t)  =  a(t)  n(t)  ,  (2.3) 

where  n(t)  is  a  zero-mean  narrowband  gaussian  process  with  covariance 
function  Rn(T)  and  a(t)  is  a  stationary,  slowly  varying  random  process, 
independent  of  n(t),  that  modulates  n(t).  This  modulating  process 
a(t)  is  further  described  as 


a(t)  =  l/b(t)  ,  (2.72) 

where  the  first-order  statistics  of  b(t)  are  specified  by  the  two-sided 
chi  distribution  given  by  Eq.  (2.8). 

The  applicability  of  the  generalized  wt"  model  as  a  description  of 
received  atmospheric  noise  was  considered  in  detail  in  Chapter  II,  and  the 
pertinent  results  can  be  summarized  as  follows:  The  first-order  statistics 
of  the  generalized  "tM  model  are  in  good  agreement  with  measured  data  on 
received  atmospheric  noise.  This  agreement  is  demonstrated  in  Figs.  4  to 
13  and  in  Fig,  15,  and  it  is  noted  that  this  agreement  is  particularly  good 
at  vlf  and  If  in  those  situations  characterized  by  low-to-moderate  local 
thunderstorm  activity.  In  addition  to  this  good  agreement  of  first-order 
statistics,  the  higher-order  statistics  of  the  generalized  ”t"  model  can 
be  specified  to  give  good  agreement  with  measured  data  on  these  statistics. 
This  result  is  demonstrated  in  Fig.  18,  where  the  plotted  data  indicate 
that  the  higher-order  statistics  of  the  modulating  process  a(t)  can  be 
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specified  (it  is  noted  that  these  statistics  are  not  yet  specified  in 
Eq.  (2.72)  above]  to  reproduce  the  experimentally  observed  probability  dis¬ 
tribution  of  the  interval  between  successive  crossings  of  a  fixed  level  by 
the  envelope  of  the  noise. 

The  investigation  into  the  applicability  of  the  generalized  "t"  model 
leads  to  the  conclusion  that  this  new  model  is  an  analytically  attractive 
model  that  is  appropriate  for  received  atmospheric  noise.  Furthermore, 
it  v.  s  pointed  out  in  Chapter  II  that  the  generalized  "t"  model  will  be 
useful  in  the  study  of  signal  detection  and  estimation  problems  in  the 
presence  of  additive  atmospheric  noise.  This  follows  from  the  fact  that 
limited  spectrum  availability  generally  dictates  that  the  received  signal 
be  observed  through  a  bandwidth  that,  in  fact,  is  substantially  smallt? 
than  the  band  center  frequency. 


B.  DETECTION  OF  KNOWN  SIGNALS  IN  ADDITIVE  ATMOSPHERIC  NOISE 

Having  developed  and  checked  the  applicability  of  the  generalized  "t" 
model  in  Chapter  II,  this  model  was  applied  in  Chapter  III  to  the  detection 
of  known  signals  in  additive  atmospheric  noise.  In  particular,  the  prob¬ 
lem  of  deciding  between  the  two  _a  priori  equally  probable  hypotheses 


h^1^;  x(t)  =  m(t)  +  y(t) 
x(t)  =  y(t) 


0  <  t  <  T 


(3.2) 


with  the  smallest  probability  of  error  was  considered  in  detail.  The  addi¬ 
tive  atmospheric  noise  was  represented  by  the  generalized  "t"  model  and 
two  cases  were  identified: 

1 .  Case  I:  Short-Duration  Signals 

This  case  is  characterized  by  the  fact  that  the  slowly  varying 
modulating  process  a(t)  in  the  generalized  "t"  model  can  be  assumed 
constant  for  the  duration  of  the  transmitted  signal.  The  familiar  corre¬ 
lator  receiver  described  by  Eq,  (3.23)  was  found  to  be  optimum  in  this 
case,  and  the  probability  of  error  achieved  by  this  receiver  was  calcu¬ 
lated  in  Eq.  (3.44).  The  probability  of  error  result  for  high  SNR  is  the 
result  of  most  interest.  It  is  given  by 
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where  is  proportional  to  the  overage  received  whitened  noise  energy 

and  E^  is  the  energy  of  the  known  signal  m(t)  at  the  output  of  the 
whitening  filter.  The  slow  decrease  with  increasing  signal  energy  of  this 
result,  relative  to  the  exponential  decrease  achieved  by  the  correlator 
receiver  in  the  presence  of  additive  gaussian  noise  (see  Fig,  22),  leads 
to  the  conclusion  that  short-duration  signals  are  not  an  optimal  choice 
in  the  presence  of  additive  atmospheric  noise.  This  is,  of  course,  phys¬ 
ically  reasonable  since  the  distinctive  impulsive  nature  of  atmospheric 
noise  suggests  the  use  of  signals  of  long  duration  relative  to  that  of  the 
typical  noise  pulse. 


2.  Case  II:  Long-Duration  Signals,  The  General  Case 

This  case  is  characterized  by  the  fact  that  the  slowly  varying 
process  a(t)  is  now  allowed  to  vary  during  the  signaling  interval.  In 
particular,  Case  Ila  generalized  "t"  noise  was  considered  in  detail. 

It  is  specified  to  be  that  Case  II  noise  whose  complex  envelope  can  be 
described  by  Eq,  (a.76)  in  terms  of  2TB  statistically  independent  samples. 
This  is  proposed  as  the  simplest  representation  of  Case  II  signals  of  prac¬ 
tical  interest.  The  optimal  receiver  in  this  case  was  found  to  be  a 
logarithmic-correlator  receiver  (see  Fig.  23)  that  implements  the  rule: 

Choose  h^  iff 


T  T 

j  in[|x(t)|2  +  y2  Jdt  >  J  in[|x(t)  -  p(t)|2  +  y2]dt  ,  (3.97) 

0  0 

where  | x (t ) j  is  the  envelope  of  x(t),  | x ( t )  -  p(t)|  is  the  envelope  of 
[x(t)  -  m(t)3,  and  y  is  proportional  to  the  average  received  noise 
power.  An  exponentially  correct  estimate  of  the  probability  of  error 
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achieved  by  this  receiver  when  the  transmitted  signal  energy  is  distributed 
uniformly  in  the  signaling  interval  is  given  by 


v0Blm  M,  - ‘W  H;  ^  •  (3-138) 

— —  00 

Q 

where  f  ( 1  /^TB )  is  a  polynomial  in  (2TB)  E  =  QT  is  the  energy  of 

the  known  signal  m(t),  and  YQ  =  2y /B  is  proportional  to  the  average 
received  noiso  energy.  The  exponential  behavior  of  this  result  is 

strikingly  superior  to  that  given  by  Eq.  (3.45)  for  short-duration  signals. 
Furthermore,  it  was  shown  (see  Fig.  24)  that  this  error  performance  pre¬ 
dicted  for  the  log-correlator  receiver  is  significantly  better  than  either 
that  achievable  by  a  matched-filter  receiver  or  that  actually  achieved  in 
practice  by  a  typical  FSK  vlf  communication  system. 

Also,  to  demonstrate  that  the  asymptotic  performance  predicted  by 
Eq.  (3.138)  for  the  log-correlator  receiver  can  be  closely  approached  in 
practice,  we  note  the  special  case  Y^b/Q  =  3.  It  was  shown  in  Chapter  III 
til,? :  this  special  case  is  reasonably  achieved  in  practice,  and  that  it 
gives  the  result: 


(3.141) 


With  the  above  probability  of  error  results  in  mind,  it  is  inter¬ 
esting  to  note  that  there  is  an  alternate  realization  of  the  log-correlator 
receiver  that  is  approximately  equivalent  to  Eq.  (3.97)  when  the  condition 
(Y  B)/(4Q)  »  1  is  satisfied.  This  realization  takes  the  form: 

Choose  iff 


Re 


(|*(t)|2.  72] 


>  j  >(<>|2 


I|*(0|2  +  7*1 


dt 


(3.98) 
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This  receiver  rule  is  of  interest,  since  it  may  be  easier  to  implement 
than  Eq,  (3.97),  and  since  it  has  a  realization  closely  related  to  the 
matched-filter  realization  of  the  correlator  receiver  rule. 

The  above  summary  of  results  on  the  general  case  of  long-duration 
signals  leads  to  the  conclusion  that  the  error  performance  predicted  by 
Eq.  (3.138)  for  the  log-correlator  receiver  in  the  presence  of  additive 
atmospheric  noise  is  a  result  of  particular  interest.  This  follows  because 
Eq.  (3.138)  gives  the  only  theoretical  estimate  known  to  the  author  of 
the  error  performance  achievable  in  the  presence  of  additive  atmospheric 
noise.  Furthermore,  the  reasonable  simplicity  of  the  log-correlator  re¬ 
ceiver  indicates  that  this  level  of  performance  can  be  approached  by  prac-. 
tically  realizable  receivers. 

C.  RECOMMENDATIONS  FOR  FUTURE  WORK 

In  Chapter  IV,  the  operation  of  the  log-correlator  receiver  was  com¬ 
pared  with  that  of  a  nonlinear  receiver  form  commonly  used  in  practice  in 
the  presence  of  additive  atmospheric  noise.  A  qualitative  discussion  was 
presented  there  to  indicate  why  this  commonly  used  receiver  form  must 
perform  in  a  suboptimal  fashion.  However,  the  analytical  difficulties 
involved  in  making  this  comparison  quantitative,  plus  the  need  for  experi¬ 
mental  verification  of  the  performance  predicted  by  Eq,  (3.138)  for  the 
log-correlator  receiver,  lead  to  the  obvious  recommendation  that  the  log- 
correlator  receiver  be  built  and  compared  experimentally  with  receiver 
forms  commonly  used  in  the  presence  of  additive  atmospheric  noise. 

Also  in  Chapter  IV,  the  statistical  estimation  problem  in  the  presence 
of  atmospheric  noise  was  briefly  considered  as  it  pertains  to  the  problem 
of  obtaining  the  "best”  (in  some  sen^e)  estimate  of  a  whistler  mode  signal 
that  has  been  contaminated  by  additive  atmospheric  noise.  As  was  pointed 
out  there,  the  scientific  importance  of  the  whistler  estimation  problem 
indicates  that  the  generalized  "t"  model  can  perhaps  be  fruitfully  ap¬ 
plied  to  the  statistical  estimation  problem  in  the  presence  of  additive 
atmospheric  noise.  Related  to  this  problem,  it  is  suggested  that  another 
problem  of  interest  is  the  two-hypotheses  communications  problem,  in  which 
the  signal  whose  presence  is  to  be  tested  is  taken  to  be  random,  i.e.,  is 
known  only  statistically.  Furthermore,  it  is  proposed  that  a  procedure 
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that  may  give  results  applicable  to  the  whistler  estimation  problem  is  to 
consider  the  two-hypotheses  testing  problem  in  which  the  random  signal  is 
taken  to  be  a  gaussian  process. 

Finally,  returning  to  the  generalized  "t"  model  itself,  it  is  pro¬ 
posed  that  an  interesting  problem  consists  of  choosing  the  covariance 
matrices  and  introduced  in  Chapter  III  in  a  manner  which  repro¬ 

duces  the  probability  distribution  of  the  inter-level-crossing  interval 
observed  in  practice  for  received  atmospheric  noise.  In  this  regard  it 
is  pointed  out  that  the  calculations  of  the  inter-level-crossing  interval 
distribution  in  Chapter  II  indicate  the  existence  of  the  desired  co- 
variance  matrices. 
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APPENDIX  A.  EVALUATION  OF  INTEGRAL 

1.  An  integral  that  appears  in  consideration  of  the  modified  gener¬ 
alized  Mt"  model  has  the  form 

CXJ 

I  =  j  x°/2  exp[ax]dx,  n  odd  >  1,  a  <  0,  b  >  0  .  (A.l) 

b 

This  integral  can  be  evaluated  by  parts,  giving  the  result 
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2.  An  integral  that  appears  in  the  Case  Ila  probability-of -error 
calculations  has  the  form 


00 

11  *  I  — 5 - 2 — - 3 - 5-  **,  W  0 

j  («2t.2)[(x-b)2..2] 


•  Jtst 


ia)(x  -  ia)(x  -  b  +  ia)(x  -  b  -  ia. 


dx  .  (A.5j 


This  integral  is  easily  evaluated  as  a  contour  integral  as  described  in 
Fig.  27. 


X  -  plane 


X  *  ilal 


X  *  b  +  i  lol 
--X 


x  *  t>-  ilol 


FIG.  27.  DESCRIPTION  OF  CONTOUR  OF  INTEGRATION  USED 
TO  EVALUATE  EQ.  (A.5). 

Note:  l)  The  contour  of  integration  C  is  closed  at 
|x|  =00. 

2)  The  poles  of  the  integrand  plotted  above 
are  for  the  special  case  b  >  0. 
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APPENDIX  B.  PROOF  THAT  EQ.  (3.126),  n'(s)  =  0,  IS  UNIQUELY  SATISFIED 

BY  sol. 

4 

In  the  calculation  of  the  Chcrnoff  bound  on  P  in  Case  Ila,  it  is 

e 

necessary  to  find  the  value  of  s (finite)  >  0  which  satisfies  the  equation 


2 

i=l 


r  r  t  2  2  2 .  f  2  2. 2  T 

f  d*  fdy  is.  +?  .tJJ - -in  -  *j^-4 - - 

J  J  ((*  -  ♦  y2  ♦  s2]  L<*  -  t*,)2  +  y  ♦  S2  J 


/ dx  / 


dy 


/  2  2  2 . 
(*  +  y  +  S  > 


s-2 


=  0  . 


J  *  r  f  \2  2  2  1 

-oo  -oo  [(x  -  +  y  +  §  J 


(B.l) 


Noting  that  the  denominator  of  the  summand  is  positive  and  finite  for  all 
i  and  all  finite  s,  it  is  first  of  all  concluded  that  each  term  of  the 
sum  must  be  zero  if  the  sum  itself  is  to  be  zero,  since  for  any  given  s 
the  sign  must  be  the  same  for  each  term  in  the  sum.  Thus  it  is  concluded 
that  we  desire  that  value  of  s  >  0  which  satisfies  the  equation 


s-2 


I(.)  £  f  d*  fdy  *r  *  5  ) — ; 

-«<  t(x  -  n1)2  +  y2  +  £2] 


in 


r  2  2  .2  1 

_ x  +  y— u - -  o 

,  \2  2.2 
(x  -  Ri)  +  y  +  I  J 


Vi  =  1,  ....  N  .  (B.2) 


It  will  now  be  demonstrated  that  (B.2)  is  uniquely  satisfied  by  s  =  1. 
Setting  s  =  1  in  (B.2),  it  is  desired  to  evaluate 
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(8.3) 


Thus  it  follows  directly  that  l(s)  |  =0,  since  the  integrand  above  is 

s=l 

seen  to  be  odd  in  x. 

Now,  to  show  that  this  solution  is  unique,  consider  the  following: 


00  00 


fs  I(s)  =  /  d*  / 


dy 


-00  -00 


f  I- 


/  2  2 

(x  +  y  +  t  ) 


in 


2  2  t2 

x  t  y  l_£ 


,  \2  2  2 

(x  ~  p.± )  +  y  +1 


_d_ 

ds 


L\(x  - 


2  .2  \ 

y  +'A _ \ 

\2  2  ,2  / 
+  y  +  ft  / 


00  00 


J  dx  J*  dy 

-OJ  -0C 


/  2  2  2  x 

+  y  +  jji 


s-2 


r,  n2  2  2  -I 

[(x  -  ji1)  +  y  +  6  J 
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in 


2  2  2 
x  ±  y .  +  ft. 


I ,  \2  2  2 

Lu  -  ni)  +  y  +1 


(B.4) 
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|  t~ 


i  ; 


i 

M, 
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Thus,  d/ds  l(s)  >  OVs,  such  that  s  =  1  is  shown  to  be  the  unique 
solution  to  ^(s)  =  °»  and  the  proof  is  thus  complete. 
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APPENDIX  C.  PROOF  THAT  INFINITE  RECEIVER  BANDWIDTH  (2B  =  «)  MINIMIZES 

THE  CASE  II  a  P  RESULT 
e 

In  the  P^  calculations  of  Case  Ha,  an  exponentially  correct  esti¬ 
mate  of  the  probability  of  error  optimized  with  respect  to  signal  shape 
is  given  by  expression  (3.135): 


(3.135) 


Noting  that  YQ  and  Q  are  not  functions  of  the  receiver  bandwidth  2B, 
the  optimal  bandwidth  can  be  determined  by  optimizing  with  respect  to  k, 
where  k  =  Y  B/Q,  Thus,  in  terms  of  this  new  parameter 

v 


(C.l) 


vhere  E  =  QT  is  the  energy  of  the  known  signal.  Now,  neglecting  the 
f (l/(/2TB)  factor,  since  it  is  of  negligible  importance  in  comparison  to* 
the  factor  that,  varies  exponentially  with  2TB,  the  optimization  problem 


reduces  to  that  of  finding  the  value  of 
function 


k  >  k  . 

—  nun 


which  minimizes  the 
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8<k)  *  [rr~  2fi(1-  5!  I1  tti)]  • 


(C.2) 


where  k  .  =  Yrt/(2E).* 

min  0'  ' 

Perhaps  the  simplest  way  to  minimize  g(k)  is  to  plot  it  vs  k, 
since  it  is  seen  in  this  way  that  in  fact  g(k)  is  a  monotonically  de¬ 
creasing  function  of  k  >  0.  From  this  it  follows  directly  that  the  opti¬ 
mal  choice  for  the  receiver  bandwidth  is  B  =  oo,  and  the  proof  is  thus 
complete.  Such  a  plot  of  g(k)  is  not  presented  here,  but  a  formula 
from  which  it  can  be  easily  prepared  is  presented  as  follows:  Using 
Erdelyi  [Ref.  49],  items  2.8(l8)  and  2.8(l6),  it  is  seen  that 


F  U  L  5.  1  \  _  F  /I  ,  .  3.  1  \ 
2  1  \  *  2’  2‘  1  +  k/  “  2  1\2*  2’  1  +  k/ 


,1  +  k; 
2 


'l  +  k)1//2  +  1 


(1  +  k)l/2  -  1 


(C.3) 


Thus: 


g(k)  = 


- 1 

2(l  + 


Hfi  +  ^  ±  ill1 

.(1  ♦  - 1 


(C.4) 


and,  for  particular  values  of  k  of  interest,  it  is  seen  that  (c.4)  gives: 

1.  lim  g(k )  =  1; 
k  -*  0 

2.  g(o. 125 )«  0.821; 

3.  g(o.5)  «  0.684; 

4.  g(5 )  «  0.544; 

5.  lim  g(k)  0.5129. 
k  -»  00 


This  minimum  valu:  results  from  the  constraint  that  the  receiver  band¬ 
width  must  be  at  least  as  large  ns  the  signal  bandwidth;  i.e.,  2B  >  l/T, 

where,  as  mentioned  in  the  te>~,  the  choice  of  T  that  minimizes  P  is 

e 

clearly  T  =  co. 
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IS.  ABSTRACT 
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/  This  work  is  concerned  with  the  development  and  application  of  an  analytical 
model  for  atmospheric  noise,  which  is  radio  noise  originating  in  lightning 
discharges.  The  generalized  "t"  mode}  resulting  from  our  approach  is  in  good 
agreement  with  measured  data  and  describes  the  received  atmospheric  noise  y(t) 
as 


y(t)  -  a(t)  n(t) > 

where  n(t)  }j  (  narrowband  gaussian  process  and  a(t)  is  a  slowly  varying  random  ' 
process,  inajpendent  of  n(t),  that ."modulates"  n(t), 

\  ' 

The  detection  of  known  signals  in  &  Iditive  generalized  "t"  noise  is  considered, - 
and  the  receiver  that  minimizes  the  probability  of  error  found  to  be  a- nonlinear 
"logarithmic-correlator"  receiver.  Analytical  and  experimental  results  indicate 
that  the  probability  of  error  achieved  by  this  receiver  in  the  presence  of 
additive  atmospheric  noise. 
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